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PREFACE 


This  final  report  containing  research  supported  by  Grant  No.  830256 
from  the  Air  Force  Office  of  Scientific  Research  (AFOSR),  Bolling  AFB, 
Washington,  0.  C. ,  consists  of  four  parts:  Part  1:  Modelling  and 
Testing  of  Idealized  Rock  Joints;  Part  2:  Further  Development  of 
Hierarchical  Plasticity  Model  for  Joints;  Part  3:  Cyclic  Multi-Degree- 
of-Freedom  Shear  Device  with  Pore  Water  Pressure;  and  Part  4:  Development 
of  Simulated  Rock  Like  Material  and  Testing.  Previous  Progress  and 
Interim  Reports  were  submitted  in  April  1984,  September  1985, 

January  1986,  October  1986  and  February  1987.  Dr.  C.  S.  Desai  acted  as 
the  principal  investigator  at  the  University  of  Arizona,  Tucson,  and  the 
project  was  monitored  by  Dr.  Spencer  Wu  and  Col.  L.  D.  Hokanson  at  AFOSR. 

A  number  of  individuals  contributed  towards  the  research  results  and 
preparation  of  these  reports:  Dr.  C.  S.  Desai,  Dr.  K.  L.  Fishman, 

Mr.  Y.  Ma,  Mr.  D.  Rigby,  Dr.  T.  Kundu  and  Mr.  S.  Wang.  Part  1  is  based 
on  Dr.  Fishman's  doctoral  dissertation;  Part  2  is  based  on  the  ongoing 
doctoral  research  of  Mr.  Ma;  Part  3  is  based  on  the  ongoing  Master's  and 
doctoral  research  of  Mr.  Rigby,  which  is  also  partly  supported  from  a 
grant  from  the  NSF;  and  Part  4  is  based  on  the  Master's  research  of 
Mr.  Wang,  which  is  guided  by  Dr.  Desai  and  Dr.  T.  Kundu.  - - - 
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ABSTRACT 


A  rather  novel  plasticity  based  constitutive  model  to  describe 
the  response  of  simulated  (rock)  joints  under  cyclic,  quasi-static  and 
static  shear  is  developed.  Development  of  the  constitutive  model  includes 
both  mathematical  formalization  based  on  the  hierarchical  approach  and 
laboratory  testing.  The  mathematical  formulation  is  such  that  the  model  k 
is  basic  and  general  and  is  capable  of  predicting  observed  behavior  of 
joints.  Laboratory  test  results  are  used ^f or  the  determination  of  parame¬ 
ters  for  the  model/  and, for  comparison  with  model  predictions. 

The  constitutive  model  is  based  on  the  theory  of  incremental  plas¬ 
ticity.  A  generalized  three-dimensional  plasticity  model  capable  of 
predicting  the  behavior  of  geologic  solid  material  such  as  soil  and  rock 
is  specialized  to  describe  the  behavior  of  individual  rock  joints.  At' 
this  time,  the  model  allows  for  effects  of  initial  normal  stress,  states 
of  shear  and  normal  stress,  plastic  hardening,  nonassociativeness ,  volume 
changes  at  joints,  and  cycles  of  loading,  unloading  and  reverse  loading. 

The  test  program  was  conducted  on  simulated  joints.  The  simulated 
specimens  were  cast  in  concrete  with  a  variety  of  surface  geometries 
(angles  of  asperities).  Specimens  were  subjected  to  a  series  of  quasi¬ 
static  and  fast  cyclic  direct  shear  tests.  Tests  were  performed  with  a 
special  device  known  as  the  Cyclic  Multi-Degree-of-Freedom  (CYMDOF)  shear 
device;  minor  modifications  of  the  device  were  necessary  for  the  testing 
of  joints.  Quasi-static  tests  included  shear  loading,  unloading  and 
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reverse  loading,  and  fast  cyclic  tests  repeated  cycles  of  shear  loading 
at  the  frequency  of  1.0  Hz.  Tests  were  conducted  under  different  levels 
of  normal  stress  and  amplitudes  of  cyclic  displacement. 

The  constants  for  elastic  and  inelastic  responses  were  found 
from  the  laboratory  test  data.  Then  typical  observed  results  were  pre¬ 
dicted  by  integrating  the  incremental  plasticity  equations  which  were 
expressed  in  terms  of  the  constants.  The  predictions,  in  general,  were 
found  to  provide  satisfactory  correlation  with  the  observations. 

The  results  of  this  research  have  demonstrated  that  the  model 
described  herein  is  capable  of  capturing  many  aspects  of  rock  joint  be¬ 
havior  during  quasi-static  and  cyclic  shear  loading.  The  model  is 
sufficiently  simplified  so  that  it  can  be  easily  implemented  in  numeri¬ 
cal  techniques  such  as  the  finite  element  method.  Such  computational 
procedures  can  be  used  to  solve  practical  boundary  value  problems  in 
rock  mechanics  involving  static  and  dynamic  loads. 


CHAPTER  1 


INTRODUCTION 

This  study  is  concerned  with  the  response  of  idealized  joints  in 
rock  masses  subjected  to  quasi-static  or  cyclic  loadings.  The  cyclic 
loading  may  be  due  to  the  propagation  of  stress  waves  during  an  earth¬ 
quake  or  due  to  blast  loads  which  can  occur  either  above  or  below 
ground.  Many  times  rock  masses  are  considered  to  be  composed  of  con¬ 
tinuous  solid  material,  i.e.  discontinuities  in  the  form  of  cracks 
(joints)  are  not  considered  to  be  present.  This  idealization  may  not  be 
valid.  Rock  masses  are  seldom  continuous,  and  discontinuities  (joints) 
of  the  order  of  several  feet  to  hundreds  of  feet  long  may  be  present 
throughout  the  rock  mass,  and  may  have  a  particular  pattern  or  orienta¬ 
tion.  These  joints  may  have  been  created  in  the  rock  mass  due  to  a 
variety  of  tectonic  processes  including  previous  shear  along  faults 
during  an  earthquake  which  may  create  a  state  of  stress  such  that  ten¬ 
sile  fracture  of  the  rock  mass  along  particular  planes  occurs.  Also, 
joints  may  become  present  during  the  formation  of  the  rock  mass  due  ro 
temperature  gradients,  etc. 

Figure  1-1  presents  some  typical  problems  encountered  in  rock 
mechanics  where  the  presence  of  joints  in  the  rock  mass  may  be  signifi¬ 
cant.  Slope  stability  analysis  in  rock  masses  involves  the  determination 
of  critical  planes  of  weakness  (joints)  along  which  sliding  may  occur, 
and  then  the  determination  of  the  state  of  stress  along  the  joint  to 


ascertain  if  failure  is  imminent.  Structures  may  be  founded  on  rock 
masses  containing  joints,  and  shafts  or  caverns  may  be  constructed 
within  jointed  rock.  The  response  of  these  systems  to  various  excita¬ 
tion  is  an  interaction  type  of  problem  in  which  the  response  of  the 
structure  is  affected  by  the  response  of  the  medium  and  vice  versa. 

Analysis  involving  jointed  rock  masses  may  be  sensitive  to  the 
properties  of  the  individual  rock  joints  (O'Conner  et  al. ,  1986).  This 
is  particularly  true  during  cyclic  loading  conditions  where  relative 
motion  along  the  joint  plane  may  dissipate  energy.  An  individual  rock 
joint  has  several  degrees-of-f reedom,  as  shown  in  Figure  1-2.  The  joint 
is  depicted  as  a  zone  of  material  with  thickness,  h,  (l-2a).  The  joint 
may  undergo  compression  or  dilation,  V  ,  due  to  changes  in  normal  stress 
across  the  joint  plane,  (l-2b).  The  joint  may  undergo  relative  tan¬ 
gential  (shear)  displacement,  U^,  due  to  the  application  of  shear  stress 
across  the  joint  plane,  (l-2c)  ,  often  accompanied  by  dilation  or  com¬ 
pression.  Also,  due  to  eccentricities,  the  joint  may  undergo  rocking 
and  experience  a  nonuniform  distribution  of  normal  stress,  (l-2d).  At 
some  points  the  normal  stress  may  become  zero  and  separation  of  the 
joint  plane  may  occur.  Separation  of  the  joint  plane  is  known  as  de¬ 
bonding,  and  if  the  joint  plane  comes  back  into  contact,  rebonding 
occurs.  Joint  planes  may  undergo  cycles  of  rebonding  and  debonding  as 
described  by  Zaman  (1982),  Nagaraj  (1986).  In  order  to  describe  the  be¬ 
havior  of  rock  joints  properly,  methods  of  analysis  must  consider  as 
many  of  these  modes  of  deformation  as  possible. 
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The  finite  element  method  is  a  very  versatile  method  which  makes 
it  possible  to  incorporate  many  of  these  factors  in  the  analysis  (Desai 
at  al. ,  1982;  Nagaraj ,  1986).  It  is  possible  to  utilize  this  method, 
and  describe  the  behavior  of  the  individual  rock  joints  through  the  use 
of  special  joint  elements.  The  joint  elements  can  be  based  on  the  idea 
of  the  thin- layer  element  proposed  by  Desai  et  al.  (1981,  1983a, 

1984b).  The  formulation  of  these  elements  is  similar  to  that  of  the 
solid  elements.  The  joint's  are  considered  to  consist  of  a  thin  zone  of 
material,  the  properties  of  which  are  determined  from  appropriate 
laboratory  tests.  The  joint  properties  are  then  described  through  a 
constitutive  matrix  which  relates  increments  of  stress  to  increments  of 
strain.  The  relevance  of  the  joint  constitutive  matrix  within  the 
framework  of  the  finite  element  method  will  be  described  in  the 
following  section. 

Relevance  of  Joint  Constitutive  Matrix 
The  governing  matrix  element  equations  for  general  dynamic  prob¬ 
lems  are  expressed  as  (Zienkiewicz ,  1971;  Desai  and  Abel,  1972;  Bathe, 
1982) 

[m]  {q}  +  [c]  {q}  +  [k]  {q}  =  {Q  (t)}  (1-1) 

where  [m]  is  the  global  mass  matrix,  [c]  is  the  global  damping  matrix, 
[k]  is  the  global  stiffness  matrix,  {q }  is  the  vector  of  nodal  displace¬ 
ments,  and  {Q  (t)}  is  the  assemblage  nodal  load  vector  which  is  a 
function  of  time.  The  overdot  denotes  derivative  with  respect  to  time. 
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Terms  in  the  global  stiffness  matrix,  [k],  involve  contributions  from 
both  solid  and  joint  elements.  Joint  element  stiffness  matrices  are 
determined  from  the  following  relationship: 

[k],  =  /  [B]T  [C]  [B]  dV  (1-2) 


where  [k]_.  is  the  joint  element  stiffness  matrix,  and  [B]  is  the  strain- 
displacement  transformation  matrix  in  nonlinear  incremental  analysis. 


[C]^  is  the  joint  constitutive  matrix  relating  increments  of  stress  to 
increments  of  strain  as  follows: 


{da}  =  [C].  {de}  (1-3) 

where  {da}  is  the  vector  of  components  of  stress  increments,  and  {de}  is 
the  vector  of  components  of  strain  increments.  The  development  of  the 
joint  constitutive  matrix  is  the  focus  of  this  research. 

Procedure 

The  constitutive  matrix,  [ C ] ^  ,  is  obtained  on  the  basis  of 
appropriate  constitutive  laws.  The  main  steps  toward  the  development  of 
appropriate  constitutive  relationships  f  a  material  are  as  follows 
(Desai  and  Siriwardane,  1984): 

1.  Mathematical  formulation 

2.  Identification  of  significant  parameters 

3.  Determination  of  parameters  from  laboratory 
tests  and  verification,  which  can  involve  the 


following  two  additional  steps: 
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4.  Successful  prediction  of  a  majority  of 
observed  data  from  which  the  parameters 
were  determined,  and  of  other  test  data 
under  different  stress  paths 

5.  Satisfactory  comparisons  between  predic¬ 
tions  from  a  solution  scheme  in  which  the 
constitutive  law  is  introduced,  and  ob¬ 
servations  or  closed-form  solutions  for 
relevant  practical  boundary  value  problems. 

In  this  investigation  the  procedure  will  be  carried  through  the 
fourth  step.  Implementation  of  the  model  in  a  numerical  (finite  element) 
program  will  be  done  at  a  later  time. 

Scope  of  Research 

The  scope  of  the  research  performed  in  this  investigation  is  as 

follows : 

1.  Perform  a  comprehensive  series  of  tests  on  individual 
idealized  rock  joints.  These  tests  will  employ  a 
special  shear  device  known  as  CYMDOF  (Cyclic  Multi 
Degree-of-Freedom) ,  (Desai,  1980a)  which  has  undergone 
minor  modification  for  rock  joint  testing. 

2.  Develop  a  constitutive  model  for  rock  joints  during 
shear  based  on  observations.  The  model  must  have 
physical  meaning  and  obey  physical  laws. 


8 


3.  Determine  parameters  for  the  model  from  the 
test  data. 

4.  Verify  the  model  by  performing  back  predictions 
of  the  test  data. 

Chapter  2  of  this  report  will  provide  a  review  of  the  literature 
in  the  area  of  rock  joints.  Chapters  3,  4,  5  and  6  provide  details  of 
the  test  program  including  descriptions  of  the  device,  instrumentation, 
calibration  of  instruments,  preparation  of  samples,  and  the  test  series. 
Chapter  7  describes  the  proposed  constitutive  model,  and  Chapter  8  de¬ 
scribes  determination  of  parameters  and  back  prediction  of  test  data. 
Finally,  Chapter  9  presents  conclusions  of  the  investigation. 


ft 
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CHAPTER  2 


LITERATURE  REVIEW 

In  this  chapter  a  review  will  be  made  of  existing  literature 
pertaining  to  rock  joints.  The  survey  will  be  divided  into  three  sec¬ 
tions.  The  first  section  will  review  testing  of  rock  joints.  The 
remaining  two  sections  will  deal  with  modeling  of  rock  joint  behavior. 
Earlier  work  in  this  area  dealt  mainly  with  predicting  slip  across  the 
joint  plane  with  little  attention  paid  to  deformation  of  the  joint 
plane;  this  research  considers  both  in  a  single  framework.  Failure 
models  described  in  section  two  are  all  that  is  necessary  to  perform 
analyses  such  as  those  based  on  limiting  equilibrium  where  the  main  in¬ 
tention  is  to  design  against  failure.  The  third  and  final  section 
describes  models  developed  to  describe  deformation  of  rock  joints 
during  shear.  Development  of  deformation  models  is  necessary  to  perform 
deformation  analysis  with  such  advanced  techniques  as  the  finite  element 
method .  ► 


Rock  Joint  Testing 

A  review  of  various  test  methods  for  the  study  of  rock  joints, 
including  laboratory  and  in  situ  techniques  is  available  in  Goodman 
(1974).  The  strength  of  rock  joints  is  usually  obtained  in  the  labora¬ 
tory  on  the  basis  of  direct  shear  tests  conducted  under  dry  state  of 
joints  or  drained  conditions.  Here  a  joint  is  subjected  to  a  normal 


stress  and  a  shear  stress  is  applied  parallel  to  the  plane  of  the 
joint.  Some  experimenters  have  employed  the  triaxial  cell  for  the 
study  of  rock  joint  behavior.  The  following  paragraphs  discuss  pre¬ 
vious  investigations  and  provide  highlights  of  results  obtained.  The 
relevance  of  this  discussion  lies  in  the  fact  that  the  development  of  a 
constitutive  model  should  be  made  on  the  basis  of  observed  behavior. 
Although  the  model  may  not  include  all  of  these  features,  one  should  be 
made  aware  and  be  conscious  of  shortcomings  inherent  in  a  model. 

Jaeger  (1959)  studied  the  frictional  properties  of  rock  joints 
by  cutting  rock  cylinders  at  various  angles  to  their  axes  and  testing 
the  cylinders  in  a  triaxial  apparatus.  Different  types  of  joints  were 
created  by  filling  the  joints  with  commercial  patching  plaster,  grinding 
the  joint  surfaces  and  placing  them  together,  or  creating  natural  shear 
surfaces.  Jaeger  found  that  the  relationship  between  the  shear  stress 
at  failure  and  normal  stress  may  be  nonlinear.  Also,  the  coefficient  of 
friction  seems  to  vary  with  inclination  of  the  joint. 

Horn  and  Deere  (1962)  found  that  the  presence  of  water  increased 
the  frictional  resistance  of  smooth  mineral  surfaces  having  massive 
crystal  structures ,  but  that  of  minerals  having  layer  lattice  structures 
was  reduced.  However,  the  antilubricating  effect  on  the  crystal  struc¬ 
tured  minerals  is  reduced  as  the  surface  roughness  increases. 

Byerlee  (1967a)  found  that  the  coefficient  of  friction  for  rock 
joints  varied  with  displacement  and  normal  stress.  Also,  the  effective 
stress  principle  was  observed  to  hold  for  rock  joints.  In  another  in¬ 
vestigation,  Byerlee  (1967b) proposed  a  theory  which  assumes  that  sliding 
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is  initiated  between  surfaces  as  a  result  of  brittle  fracture  of 
asperities. 

Krsmanovic  (1967)  ran  a  series  of  direct  shear  tests  on 
samples  of  sandstone,  conglomerate  and  limestone  in  order  to  study  the 
initial  and  residual  strength  of  hard  rock.  Initially,  the  samples 
contained  no  joints  or  discontinuities,  and  it  was  found  that  the  ini¬ 
tial  shear  strength  of  the  specimen  was  reached  after  very  small 
deformations.  At  this  point  a  sharp  decrease  in  shear  strength 
occurred  as  deformation  continued  along  a  developed  shear  plane.  After 
large  deformations,  a  residual  shear  strength  was  attained.  It  was 
found  that  the  locus  of  points  defining  the  residual  shear  strength  of 
rock  joints  is  linear  with  respect  to  normal  stress  and  passes  through 
the  origin. 

Locher  (1968)  developed  an  apparatus  to  determine  the  friction 
angle  of  rock  joints  subjected  to  a  quasi-static  shear  load.  Displace¬ 
ment  normal  to  the  joint  plane  was  constrained.  The  device  displaced 
the  specimens  in  both  forward  and  reverse  directions,  and  the  normal 
stress  was  allowed  to  vary  throughout  the  test  due  to  the  tendency  of 
the  joint  to  compress  or  dilate.  Locher  found  that  the  angle  of  fric¬ 
tion  for  rock  joints  is  in  the  range  of  that  for  parent  media  such  as 
soil  and  rock,  the  difference  in  strength  coming  from  a  "cohesion"  term 
in  the  failure  criteria. 

Locher  and  Rieder  (1970)  ran  a  series  of  quasi-static  direct 
shear  tests  on  discontinuities  within  Jurassic  Limestone.  Laboratory 
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tests  with  and  without  a  filler  material  were  performed,  as  well  as 
in  situ  tests.  The  filler  material  decreased  the  friction  angle  by 
about  three  degrees,  and  the  laboratory  friction  angles  were  about 
five  degrees  less  than  those  determined  in  situ. 

Coulson  (1972)  ran  a  series  of  quasi-static  direct  shear  tests 
on  flat  rock  surfaces  having  different  "micro- roughnesses"  due  to 
grinding,  polishing  or  sandblasting.  Different  types  of  surface  damage 
were  studied  and  correlated  with  observed  variations  of  the  coefficient 
of  friction  with  displacement.  The  effects  of  water  on  the  surfaces 
were  also  investigated.  The  initial  coefficient  of  friction  was  ob¬ 
served  to  increase  with  an  increase  in  "micro- roughness"  and  to  decrease 
at  higher  levels  of  normal  stress.  At  low  levels  of  normal  stress  the 
residual  shear  strengths  of  rougher  surfaces  were  higher,  but  as  the 
normal  stress  was  increased,  the  residual  strength  was  essentially  con¬ 
stant  with  respect  to  "mic ro- roughness . " 

Christensen  et  al.  (1974)  studied  the  dependency  of  the  co¬ 
efficient  of  friction  for  rock  joints  on  the  state  of  stress.  A  rotary 
shear  apparatus  was  employed  which  could  apply  normal  stress  as  well  as 
confining  pressure  around  the  joint.  For  the  same  ratio  of  confining 
stress  to  normal  stress,  the  coefficient  of  friction  decreased,  and  the 
joint  stiffness  increased  with  increasing  normal  stress.  By  comparing 
the  states  of  stress  in  joints  during  shear  to  the  failure  locus  for 
intact  material,  it  was  postulated  that  the  closer  the  state  of  stress 
on  the  rock  joint  during  shear  to  the  failure  locus  of  the  intact  mate¬ 
rial,  the  higher  the  coefficient  of  friction  for  the  joint. 


Kutter  (1974)  developed  a  rotary  shear  apparatus  in  order  to 
study  the  effect  of  large  displacements  on  the  shear  strength  of  rock 
joints.  A  schematic  of  the  device  is  given  in  Fig.  2-1.  The  device  is 
displacement  controlled,  and  the  rotation  is  applied  to  the  upper 
sample  through  a  chain  and  sprocket  assembly.  The  lower  sample  is  pre¬ 
vented  from  rotating  but  is  allowed  to  pitch  and  roll  during  shear, 
allowing  full  contact  of  the  shearing  surface,  and  eliminating  ex¬ 
cessive  dilation  due  to  an  inclined  discontinuity.  The  torque  at  the 
interface' is  measured  with  strain  gages  mounted  to  the  reaction  arm 
that  prevents  the  rotation  of  the  lower  part  of  the  sample.  Kutter’ s 
preliminary  study  showed  that,  in  general,  rock  joints  do  not  seem  to 
require  large  displacements  to  reach  residual  strength. 

Martin  and  Miller  (1974)  used  direct  shear  equipment  to  test 
weathered  rock  or  rock  joints,  and  reported  results  from  strength  tests 
of  joints  in  sandstone  with  different  degrees  of  weathering.  A  sche¬ 
matic  of  the  device  is  shown  in  Fig.  2-2.  The  results  showed  decrease 
in  strength  with  increase  in  degree  of  weathering.  The  study  also 
suggested  a  procedure  for  isolating  the  dilatancy  in  the  rock  joints. 
They  adopted  the  model  proposed  by  Ladanyi  and  Archambault  (1970)  to 
describe  the  dilatant  component  of  the  joint  strength. 

Mogilevskaya  (1974)  found  that  the  effect  of  wetting  the  shear 
plane  on  the  strength  or  rock  joints  depended  on  the  mineralogy.  For 
crystalline  rock,  there  seemed  to  be  no  change,  but  for  other  types 
such  as  sedimentary  or  intrusive  rock,  there  was  a  decrease  in  the 
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Diagram  of  Direct  Shear 
[Martin  and  Miller  (1974)] 
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initial  friction  angle.  If  a  filler  material  exceeded  the  height  of 
the  macro-asperities,  then  the  shear  strength  of  a  rock  joint  depended 
on  the  properties  of  the  filler,  and  not  on  the  joint  morphology. 

Nikitin  et  al.  (1974)  studied  the  effects  of  impulse  loading  on 
rock  joints,  both  with  and  without  a  filler  material.  Impulse  loads 
were  induced  in  a  direct  shear  apparatus  by  applying  a  shear  stress 
which  was  near  failure  for  a  joint  under  a  given  normal  stress,  and 
then  reducing  the  normal  stress  abruptly.  This  caused  the  upper  sample, 
which  was  allowed  to  translate,  to  accelerate.  The  results  indicated 
that  for  joints  both  with  and  without  a  filler,  initially  under  impulse 
loading,  the  she^r  strength  was  higher  than  the  peak  "static"  shear 
strength.  £  relative  displacements  occurred,  the  shear  strength  under 
impulse  loading  decreased,  and  became  less  than  peak  static  strength. 

Pratt  et  al.  (1974)  performed  both  laboratory  direct  shear  and 
in  situ  shear  tests  on  samples  of  jointed  quartz  diorite.  The  effects 
of  specimen  size,  loading  path,  and  the  orientation  of  the  joint  on  the 
shear  strength  were  studied.  The  specimens  with  larger  joint  areas  had 
lower  shear  strengths.  This  could  be  due  to  differences  between  the 
total  area  of  the  joint  and  the  'apparent'  contact  area.  As  the  angle 
of  the  joint  with  respect  to  the  applied  load  increased,  the  shear 
strength  increased.  The  shear  strength  seemed  to  be  affected  by  load 
path.  Angles  of  friction  determined  from  direct  shear  tests  with  con¬ 
stant  normal  stress  were  different  than  those  obtained  from  proportional 
loading  tests. 
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Based  on  friction  tests  on  plaster  casts  of  rough  fractured 
surfaces  of  different  rock  types  (granite,  sandstone  and  limestone), 
Schneider  (1974)  found  that  the  geometry  of  the  joint  surface  has  sig¬ 
nificant  influence  on  the  stress-strain  response,  and  the  residual 
strength  is  dependent  upon  the  surface  geometry  and  material  parameters. 

Ross-Brown  and  Walton  (1975)  describe  a  portable  direct  shear 
box  which  is  schematically  shown.  Fig.  2-3.  The  device  was  designed 
such  that  as  much  data  as  possible  may  be  gathered  in  the  field.  How¬ 
ever,  the  device  is  not  designed  for  in  situ  testing.  It  can  be  deduced 
from  the  figure  that  the  shear  and  normal  stresses  are  applied  by  hy¬ 
draulic  rams  which  react  against  cables  attached  to  the  shear  box. 

Dieterich  (1978a,  1978b,  1979a,  1979b,  1981a,  1981b)  and 
Dieterich  and  Okubo  (1981)  presented  models  and  test  results  for  de¬ 
scribing  fault  behavior  including  stick-slip  phenomenon  for  predicting 
earthquakes.  Although  the  basic  mechanics  aspects  of  these  studies  are 
relevant  to  the  mechanisms  of  joints  or  interfaces,  the  present  study  is 
more  concerned  with  effect  of  blast  or  earthquake-type  loadings  on  joints, 
resulting  in  various  modes  of  deformation.  Hence,  detailed  reviews  of 
these  studies  are  not  included. 

Kane  (1981)  presented  a  review  of  static  strength  behavior  of 
rock  joints  and  reported  direct  shear  test  results  for  joints  with 
different  roughnesses  in  a  granite  gneiss.  Tests  were  performed  using 
natural  joints  and  joints  created  by  fracture  along  existing  cleavage 


planes . 


19 


The  foregoing  models  involved  tests  for  dry  joints  without 
pore  fluid  (water)  'ressures.  Goodman  and  Ohnishi  (1973)  developed  a 
static  direct  shear  device,  Fig.  2-4,  that  permitted  undrained  shear 
testing  of  joints  with  measurement  of  pore  water  pressures.  They  found 
that  the  behavior  of  joints  was  affected  by  compression  and  closure  of 
joint,  dilation  and  the  state  of  stress. 

Most  of  the  previous  studies  are  concerned  with  static  behavior 
of  dry  (or  wet)  joints  with  major  attention  to  strength  behavior  of 
rock  joints.  Recently,  there  are  a  few  studies  that  consider  cyclic 
behavior  of  joints  and  effect  of  pore  water  pressure.  Crawford  and 
Curran  (1981)  developed  a  direct-shear  device  to  examine  the  effect  of 
shear  rate  on  the  shear  resistance  of  smooth  artificial  joints  in  dif¬ 
ferent  rocks.  Gould  (1982)  and  Gillette  (1983)  describe  a  cyclic 
direct-shear  device.  Fig.  2-5,  and  its  use  for  studying  undrained  be¬ 
havior  of  rock  joints. 


Failure  Models 

Reviews  on  shear  strength  of  rock  joints  are  presented  by 
various  investigators,  Jaeger  (1959,  1971);  Barton  (1973,  1974); 
Goodman  (1974,  1976);  Jaeger  and  Cook  (1976);  Kane  (1981);  the  presen¬ 
tation  herein  has  consulted  these  publications. 

Amonton's  law  (1699)  for  metals  provides  a  basis  for  defining 
shear  force,  F,  parallel  to  the  contact  area  required  to  initiate 
sliding: 


F  =  UN 


(2-1) 
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where  N  =*  normal  force  and  p  ■  coefficient  of  friction,  which  is  inde¬ 
pendent  of  N  and  the  contact  area  A.  The  work  of  Bowden  and  Taber 
(1950)  for  metals  has  also  provided  a  basis  for  defining  rock  friction. 
In  these  theories  for  metals,  the  shearing  force  is  derived  from  force 
required  for  plowing  (force  necessary  to  plow  in  grooves  at  contact) 
and  that  required  to  shear  asperities.  Since  in  the  case  of  rocks  and 
minerals  failure  is  often  caused  by  brittle  fracture  with  plastic  flow 
and  'welding'  of  asperities,  the  above  mechanism  for  metals  may  not  be 
strictly  valid  for  rock  joints,  Jaeger  (1959). 

In  order  to  allow  for  elastic  behavior,  Archaud  (1958)  proposed 
the  following  expression  for  F: 

F  ■  UN01  (2-2) 

where  m  =  exponent  whose  value  varies  between  2/3  and  1;  for  purely 
elastic  behavior,  m  =  2/3  and  for  ductile  contact,  m  =  1.  For  the 
latter,  Eq.  (2-2)  reduces  to  Eq.  (2-1). 

Jaeger  and  Cook  (1976)  presented  the  following  Mohr-Coulomb 
criterion  for  shear  strength: 

t  »  s  +  pa  (2-3) 

o  n 

where  sq  =  shear  strength  (adhesion)  of  contact  surface  analogous  to 
cohesion  c  in  solids.  Lajtai  (1970)  suggested  that  the  linear  model  in 
Eq.  (2-3)  may  not  be  sufficient,  and  he  proposed  failure  mechanisms  in¬ 
volving  stress  dependency  and  components  of  stress  tensor. 
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Most  rocks  exhibit  strain  softening  type  behavior.  Fig.  2-6, 
where  a  residual  strength  is  reached  after  the  peak  strength  (Goodman, 
1974,  1976;  Jaeger,  1971;  Hoek  and  Bray,  1974).  The  schematic  of  peak 
and  residual  envelopes  is  shown  in  Fig.  2-7.  The  residual  strength  is 
considered  to  represent  the  strength  of  an  uncemented  joint.  The  ratio 
of  peak  to  residual  strength  is  defined  as  stress  ratio,  Sr;  for  a 
filled  joint,  the  ratio  can  be  as  high  as  4,  Barton  (1973). 

A  question  arises  as  to  the  definition  of  "failure"  because  its 
adoption  may  depend  on  various  factors  such  as  type  of  loading  and 
structure  to  be  designed  and  nature  of  the  (joint)  material.  Giuseppe 
(1970)  proposed  various  possible  definitions:  (a)  maximum  or  peak 
stress,  (b)  "knee"  of  the  T  vs.  uy  curve  (where  slope  change  can  occur), 
(c)  "inversion"  of  vertical  deformation  (prior  to  failure)  or  (d)  cur¬ 
rent  deformation  greater  than  previous  ones. 

Since  the  strength  behavior  of  joints  can  be  nonlinear, 

Jaeger  (1959,  1971,  1976)  proposed  the  following  relation.  Fig.  2-8: 

-ba 

x  *  c.  (1  -  e  n)  +  a  tan  $  (2-4a) 

p  j  n  r 

For  very  large  values  of  a  ,  Eq.  (2-4a)  becomes 

x  =  c.  +  o  tan  (2-4b) 

p  j  n  r 

and  for  a  =  0,  it  gives 
n 

dx 

_£  = 

da 

n 


c,b  +  tan  <J) 

j  r 


(2-4c) 


Figure  2-8.  Empirical  Shear  Strength  Curve 
[Jaeger  (1971)] 
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Until  now  the  discussion  has  been  limited  to  planar  contact 
surfaces  with  shear  stress  parallel  to  it.  Now  consider  a  contact 
surface  such  that  it  is  not  parallel  to  the  shear  stress,  but  inclined 
at  an  angle  i.  Fig.  2-9.  Then  the  expression  for  shear  stress  becomes 
(without  sq),  Patton  (1966) 

x  =*  a  tan  (<£  +  i)  (a  <  a  )  (2-5a) 

n  rU  v  n  n 

x  =  c  +  a  tan  <})  (2-5b) 

n 

* 

where  =  critical  or  yield  stress,  and  cj>^  =*  mineral  friction  angle 
between  smooth  and  even  surfaces;  4>  =*  4>r  3  residual  friction  angle;  at 
larger  shear  strains,  <J>^  a  <j)^.  It  was  noted  by  Patton  that  an  apparent 
cohesion,  c,  can  exist  due  to  interlocking  of  asperities.  Figure  2-10 
shows  the  bilinear  failure  envelope  proposed  by  Patton.  One  of  the 
characteristics  of  this  model  is  that  it  predicts  normal  (relative)  dis¬ 
placement,  v^,  for  any  horizontal  (relative)  displacement,  u^;  in  other 
words,  during  ur>  the  volume  of  the  specimen  increases  or  it  dilates. 

The  angle  i,  which  denotes  the  slope  (dv  /dup  of  v^  vs.  u^  plot,  can  be 
used  to  define  the  measure  of  the  dilation.  The  tendency  to  dilate  de¬ 
creases  with  increasing  normal  stress.  The  dilation  approaches  zero  as 
the  normal  stress  approaches  the  unconfined  compressive  strength  of  the 
asperities.  At  this  stage,  the  asperities  shear  off  rather  than  be 
ridden  over  (Barton,  1973;  Goodman,  1974,  1976;  Martin  and  Miller,  1974). 

Barton  (1973)  and  Barton  and  Chouby  (1977)  proposed  the  following 
model  for  strength  of  rock  joints: 
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x  =  an  tan  (JRC  log  (-£-  +  4>b>) 

n 


(2-6) 


where  JRC  =  joint  roughness  coefficient,  which  is  an  empirical  value 
between  0  and  20  (for  smooth  and  rough  joint),  JCS  *  joint  wall  (uncon¬ 
fined)  compressive  strength  for  low  a  and  confined  for  high  CT  ,  <j>,  =* 

n  n  b 

residual  shear  strength  of  flat,  nondilatant  rock  (dry  or  wet)  surfaces. 
A  plot  of  Barton's  model  is  shown  in  Fig.  2-11.  In  the  context  of 
plasticity,  the  nonassociative  rule  at  yield  (which  represents  dilation 
rate  at  peak)  from  Barton's  model  can  be  expressed  as 


n  „  .JRC  .  ,JCS,, 

IT  =  tan  (~  108  (~}) 

s  n 


(2-7) 


where  de^,  d£g  =  increments  of  normal  and  shear  strains,  respectively. 
The  dilatant  behavior  of  a  joint  is  affected  by  factors  such  as  width  of 
asperities,  roughness  and  angle  of  asperities. 

Ladanyi  and  Archambault  (1970)  developed  the  following  model  for 
rock  joint  behavior  and  yield  surface,  Fig.  2-12: 


a  (1  -  a  )  (tan  i  +  tan  $  )  +  a  S 
n  _ s  _ p  s  r 

1  -  (1  -  a  )  tan  i  tan  d>, 
s  f 


(2-8a) 


where  ag  ■  relative  length  of  shear  through  asperities,  =  statistical 
average  of  friction  angle  with  varying  i  (cf)^  s  <j>^  +  i  ave.),  and  Sr  = 
shear  strength  of  asperities.  The  flow  rule  with  this  model  is  given  by 


WWl  ill  >K«|)l«fN  !m 


Figure  2-11.  Barton  Model  (1973) 
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dc  P 

— —  =  tan  i  for  a  <  a  (2-8b) 

,  p  n  nt 

de  y 
s 


where  a  =  stress  at  which  the  asperities  are  sheared  off  and  i  s  0, 
nt 

and  the  superscript  p  denotes  plastic  strain.  When  =*  a  ,  the  joint 

can  be  considered  to  have  reached  a  critical  state  and  the  behavior  can 

be  considered  to  transit  from  that  of  joint  to  that  of  intact  rock. 

Equation  (2-8b)  defines  a  set  of  plastic  potentials  and  since  they  are 

not  equal  to  the  function  describing  the  yield  surface,  the  flow  will  be 

nonassociative;  hence,  tan  i  j4  . 

n 

Models  for  Deformation  of  Rock  Joints 
Goodman  et  al.  (1968)  proposed  a  model  for  rock  joint  in  which 
the  constitutive  or  stiffness  relation  for  the  joint  element  is  ex¬ 
pressed  as 


(2-9) 


where  k  ,  k  =  stiffness  for  coupling  between  shear  and  normal  re¬ 
sponses.  In  many  applications,  it  is  assumed  k^  »  k  =  0.  Schematic 
force  displacement  curves  for  the  stiffnesses  are  shown  in  Fig.  2-13, 
where  w^  (vf)  and  w^  (u^.)  are  (relative)  normal  and  shear  displacements. 

Ghaboussi  and  Wilson  (1973)  utilized  a  constitutive  model  which  in¬ 
volved  normal  and  shear  stiffnesses  as  above  which  can  be  linear  or  can 
vary  nonlinearly  (with  stress).  They  also  proposed  that  dilatancy  in 
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Figure  2-13-  Data  From  Hypothetical  Direct 
Shear  Test  on  a  Rock  Joint 
(Goodman  et  al  (1966)] 
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rock  joints  can  be  incorporated  by  using  the  idea  of  the  Cap  model 
(DiMaggio  and  Sandler,  1971),  (Fig.  2-14).  The  fixed  "failure"  en¬ 
velope,  f £ ,  and  the  yield  caps,  are  expressed  as 

ff  (V  °s)  =  0  (2- 10a) 

f  (a  ,  a  ,  k)  =»  0  (2-10b) 

c  n  s 

where  <  =  strain  hardening  parameter.  Although  the  above  plasticity 
model  was  implemented  for  boundary  value  problems ,  the  question  of  de¬ 
termination  of  required  parameters  remained. 

Zienkiewicz  et  al.  (1977)  proposed  an  elastic-viscoplastic 
model  for  joint  behavior.  A  Mohr-Coulomb  yield  criteria  was  employed 
with  both  associative  and  nonassociative  flow  rules.  The  nonassociative 
flow  rule  adopted  a  potential  function,  Q,  of  the  form 

Q  =  d  +  d  tan  ip  -  <  =  0  (2-11) 

pq  q 

where  d  ■  shear  stress,  3  =  normal  stress,  i V  (between  0  and  tf>)  the 

pq  q 

internal  angle  of  friction,  and  k  =  parameter. 

Roberds  and  Einstein  (1978)  presented  a  model  that  can  include 
elastic,  plastic  and  viscous  effects  in  rock  discontinuities.  The  con¬ 
cept  is  based  on  using  two  or  multisurface  plasticity  models  such  as 
critical  state  or  cap.  The  model  appears  to  be  essentially  in  a  quali¬ 
tative  form,  and  includes  qualitative  comparisons  with  other  available 


models  for  joints. 
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Figure  2-14.  Capped  Plasticity  Model  for  Joints 
[DiMaggio  and  Sandler  (1971), 
Ghaboussi  et  al  (1973)] 
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Hsu-Jun  (1979)  proposed  a  plasticity  based  model  for  joints  in 
which  the  Mohr-Coulomb  criterion  is  used  as  the  yield  function. 

Pande  (1979)  presented  a  model  for  rock  mass  behavior  derived 
from  individual  representation  of  joint  behavior  by  using,  for  instance. 
Barton's  equation  (Eq.  2-6),  and  "solid"  rock  behavior  by  using  a  yield 
function  expressed  in  terms  of  stress  invariants. 

Heuze  and  Barbour  (1982)  presented  a  review  of  some  of  the 
available  models  for  joints. 

Plesha  (1986)  describes  an  incremental  plasticity  model  which 
may  be  used  to  model  joints  with  dilatancy  and  contact  surface  degrada¬ 
tion  (damage) .  The  model  distinguishes  between  macroscopic  and 
microscopic  features  of  a  discontinuity.  Through  macroscopic  considera¬ 
tions,  an  incremental  constitutive  law  is  derived  which  is  applicable  to 
a  large  class  of  contact  friction  problems.  By  idealizing  the  micro¬ 
structure  to  consist  of  interlocking  asperity  surfaces,  the  constitutive 
equations  are  specialized  for  the  description  of  rock  joints  including 
effects  such  as  dilatancy,  asperity  surface  degradation  and  bulking. 
Asperity  degradation  is  described  using  a  simple  tribological  relation¬ 
ship  which  assumes  that  degradation  is  a  function  of  the  plastic 
tangential  work,  WtP. 

Barton  and  Bandis  (1987)  modified  Barton’s  model  further  to  con¬ 
sider  the  variation  in  shear  strength  and  dilation  angle  during  the 
course  of  shearing.  These  effects  are  incorporated  on  the  basis  of  a 
"mobilized  roughness"  concept.  Roughness  is  mobilized  at  the  onset  of 
dilation,  has  its  maximum  contribution  at  peak  shear  displacement,  and 
then  declines  due  to  surface  mismatch  and  wear. 
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Drumm  and  Kane  (1987)  describe  the  nonlinear  inelastic  response 
of  rock  joints  utilizing  an  incremental  plasticity  based  model,  similar 
to  the  Cap  model  of  DiMaggio  and  Sandler  (1971),  but  formulated  in  terms 
of  natural  joint  coordinates.  Elliptical  caps  or  moving  yield  surfaces 
are  determined  from  experimentally  observed  dilation  rates.  The  cus¬ 
tomary  fixed  yield  surface  is  replaced  by  a  pair  of  slip  surfaces  which 
describe  peak  and  residual  strengths  of  the  joint  plane. 

Divakar  et  al.  (1987)  model  the  shear  behavior  of  rough  inter¬ 
faces  using  a  crack  stiffness  matrix  formulation.  The  behavior  is 
considered  to  be  nonlinear  elastic,  and  expressions  are  formulated  for 
shear  stress  vs.  slip  and  slip  vs.  dilation  at  constant  normal  stress. 
Coupling  is  incorporated  by  relating  normal  stresses  to  peak  shear 
stresses  by  means  of  a  "failure  surface. "  Introducing  this  coupling 
effect  provides  off-diagonal  terras  in  the  elastic  stiffness  matrix. 

From  the  above  review,  it  appears  that  major  attention  in  the 
past  has  been  given  to  characterizing  strength  behavior  of  joints. 
Limited  attention  has  been  given  to  define  deformation  behavior  by  using 
linear  or  piecewise  linear  idealization,  and  this  is  limited  to  cases  of 
monotonic  loading.  There  are  hardly  any  models  that  include  the  entire 
deformation,  failure  and  softening  in  a  single  framework.  A  major  ob¬ 
jective  of  this  research  is  to  begin  to  develop  such  a  general  concept 
on  the  basis  of  laboratory  tests,  including  quasi-static  and  cyclic 
loadings  on  simulated  joints  in  concrete. 


CHAPTER  3 


DESCRIPTION  OF  TEST  DEVICE 

Large-scale  fast  and  slow  or  quasi-static  direct  shear  tests 
were  performed  on  (idealized)  rock  joint  surfaces  using  the  translation 
portion  of  the  cyclic  multi-degree-of-f reedom  shear  device  (CYMDOF) . 

The  CYMDOF  device  is  described  by  Desai  (1980a,  1981),  Drumm  and  Desai 
(1983),  and  Desai  et  al.  (1983b).  Minor  modification  of  the  test  box  was 
undertaken  such  that  vertical  displacement  and  rotation  of  the  upper 
sample  about  an  axis  in  the  plane  of  shear  perpendicular  to  the  direc¬ 
tion  of  the  shear  were  permitted. 

The  device,  as  shown  in  Fig.  3-1,  consists  of  an  upper  and  lower 
sample  holder,  a  load  frame,  a  system  for  guiding  the  motion  of  the 
upper  sample  holder,  a  hydraulic  system  for  applying  loads,  a  control 
system,  a  measurement  system  for  measuring  loads  and  displacements,  and 
a  data  acquisition  system. 

Upper  and  Lower  Sample  Holders 

The  upper  and  lower  sample  holders  are  shown  in  Figs.  3-2  (a) 
and  (b) .  The  upper  sample  holder  is  smaller  than  the  lower  sample 
holder  so  that  during  shear  testing  no  loss  of  contact  area  of  the  upper 
sample  shear  surface  will  occur.  Also  depicted  are  the  plates  required 
for  the  rail  system  which  guides  the  upper  sample  motion  with  respect  to 
the  lower  sample  while  allowing  for  appropriate  degrees-of-f reedom. 


The  upper  sample  holder  has  a  pair  of  plates  protruding  near  the 
bottom,  centered  on  the  front  face,  for  attachment  of  the  horizontal  ac¬ 
tuator  which  exerts  the  shear  force  on  the  sample.  These  plates  are  as 
close  to  the  bottom  of  the  sample  holder  as  possible  in  order  to  reduce 
the  moment  that  the  shear  force  exerts  on  the  upper  sample.  Holes  at 
the  front  of  the  sample  holder  allow  instrumentation  wiring  to  pass 
through. 

The  lower  sample  holder  has  a  cutout  portion  on  one  side  and  in 
the  front,  although  only  the  front  cutout  is  shown  in  Fig.  3-2  (b) . 

These  cutouts  allow  for  observation  of  the  joint  during  shear  testing, 
and  also  provide  access  to  instrumentation  to  be  discussed  later.  A 
leveling  plate  was  placed  on  the  inside  bottom  of  the  holder  such  that 
the  bottom  sample  could  be  leveled  prior  to  testing.  The  bottom  sample 
holder  is  bolted  at  the  bottom  to  the  test  frame  so  that  it  is  held 
fixed  during  testing. 

Securing  Samples  to  Sample  Holders 

Samples  are  secured  to  sample  holders  with  a  commercially  avail¬ 
able,  quick-setting,  nonshrinking  cement  known  as  "Burke  Stone," 
manufactured  by  the  Burke  Company,  San  Mateo,  California.  Samples  are 
cast  with  slightly  undersized  dimensions  with  respect  to  the  sample 
holders.  Upper  and  lower  sample  holders  have  inside  dimensions  (L  x  W) 
of  12  in.  x  12  in.  (30.5  cm.  by  30.5  cm.)  and  16  in.  x  16  in. 

(40.64  cm.  x  40.64  cm.),  respectively,  while  upper  and  lower  samples 
have  inside  dimensions  of  11.5  in.  x  11.5  in.  (29.2  cm.  x  29.2  cm.)  and 
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15.5  in.  x  15.5  in.,  respectively.  This  allows  the  samples  to  be 
leveled  and  centered  within  the  sample  holders,  then  cemented  into 
place  by  filling  the  gap  between  the  sample  and  sample  holder  with  Burke 
Stone.  Since  removal  of  the  samples  from  the  sample  holders  after 
testing  would  be  very  difficult  if  the  entire  gap  is  filled  with  Burke 
Stone,  only  intermittent  lengths  of  the  gap  are  filled,  forming  cement 
plugs  around  the  periphery  of  the  specimen.  The  cement  plugs,  depicted 
in  Fig.  3-3,  are  1.0  to  1.5  in.  deep  and  1.0  to  1.5  in.  (2.54  to  3.81  cm) 
long.  The  spaces  between  the  cement  plugs  are  used  as  an  aid  in  re¬ 
moving  the  plugs  after  testing. 

The  top  sample  is  placed  in  the  upper  sample  holder  in  such  a 
manner  that  both  the  sample  and  the  holder  are  level  prior  to  placing 
the  cement  plugs.  Initially  the  top  sample  is  set  on  a  level  plate  with 
its  structured  surface  facing  down,  and  the  upper  sample  holder  is  sus¬ 
pended  above  the  lower  sample  by  means  of  a  winch.  Fig.  3-4.  The  upper 
sample  holder  is  lowered  around  the  upper  sample  until  contacting  two 
protruding  metal  plates  at  the  bottom  side  of  the  sample.  Set  screws 
are  placed  around  the  periphery  for  leveling  the  upper  sample  holder. 
After  centering  and  leveling  the  upper  sample  holder,  cement  plugs  are 
placed  around  the  upper  periphery  as  depicted  in  Fig.  3-3.  After 
allowing  the  cement  to  set,  the  upper  sample  and  holder  are  overturned 
and  cement  plugs  are  placed  around  the  bottom  periphery  as  well. 

The  bottom  sample  is  placed  in  the  bottom  sample  holder  in  such 
a  manner  that  it  is  level  with  respect  to  the  horizontal, 


Upper  Sample 


Lower  Sample 


Figure  3-4.  Setting  of  Upper  Sample  Holder  to  Upper  Sample 


and  the  asperities  on  the  surfaces  of  the  upper  and  lower  samples  will 
mate  when  the  device  is  assembled.  The  bottom  sample  is  lowered  into 


the  bottom  sample  holder  by  means  of  a  winch  (Fig.  3-5).  A  level  plate 
at  the  base  of  the  bottom  sample  holder  is  used  to  level  the  sample. 

Tne  position  of  the  bottom  sample  is  adjusted  such  that  asperi¬ 
ties  on  the  surfaces  of  the  top  and  bottom  samples  run  in  the  same 
direction  after  assembling  the  device.  This  is  necessary  for  proper 
mating  of  the  joint  at  the  onset  of  testing.  The  rear  edge  of  the  upper 
sample  holder  is  used  as  a  reference  and  tiie  angle  of  an  asperity  trough 
on  the  surface  of  the  upper  sample  with  respect  to  this  edge  is  deter¬ 
mined  as  shown  in  Fig.  3-6.  Then  the  bottom  sample  is  positioned  in  the 
bottom  sample  holder  such  that  its  asperity  peaks  make  the  same  angle 
with  respect  to  a  reference  line  on  the  bottom  sample  holder.  This 
reference  line  is  such  that  it  will  be  parallel  to  the  rear  edge  of  the 
upper  sample  holder  after  assembling  the  device.  After  positioning  and 
leveling  the  lower  sample,  cement  plugs  are  placed  around  the  periphery. 

Rail  System 

The  device  is  designed  to  shear  an  interface  or  rock  joint  by 
applying  shear  and  normal  forces  to  a  top  sample.  With  respect  to  the 
bottom  sample,  which  is  held  fixed,  the  top  sample  is  allowed  to  displace 
horizontally,  vertically  and  rotate  about  an  axis  in  the  plane  of  shear 
perpendicular  to  the  direction  of  shear.  Top  and  bottom  samples  are 
secured  to  sample  holders  which  are  then  placed  in  a  load  frame  containing 
actuators  for  applying  vertical  and  horizontal  forces  as  shown  in 
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Figure  3-6.  Determination  of  the  Angle  of  an  Asperity 
Trough  with  Respect  to  the  Rear  Edge  of 
the  Upper  Sample  Holder 
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Fig.  3-7.  Motion  of  the  upper  sample  is  guided  by  means  of  a  rail 
system  depicted  in  Fig.  3-8.  The  rail  system  consists  of  two  1  in. 

(2.54  cm.)  hardened  steel  rods,  one  on  either  side  of  the  sample 
holders,  which  are  constrained  by  four  slotted  bearings,  two  on  each 
side  of  the  bottom  sample  holder.  The  slotted  bearings  are  made  from 
oil  impregnated  bronze.  Details  of  the  slotted  bearings  can  be  seen, 
Fig.  3-8.  These  bearings  allow  translation,  vertical  (upward  or  down¬ 
ward)  displacement  up  to  0.25  in.  (0.635  cm.),  and  rotation  about  an 
axis  in  the  plane  of  shear,  perpendicular  to  the  direction  of  shear,  of 
the  top  sample  with  respect  to  the  bottom  sample.  Lateral  motion  of  the 
upper  sample  or  rotations  other  than  those  considered  above  are 
constrained. 


Load  Frame 

The  load  frame  is  depicted  in  Fig.  3-7.  Wide  flange  beams 
(W  10  x  45)  were  used  in  the  construction  of  the  load  frame.  The  frame 
was  designed  to  withstand  a  vertical  or  horizontal  load  up  to  30  tons 
and  a  frequency  of  about  2  Hz. 

The  load  frame  is  self-supporting  and  no  attachment  to  the  floor 
is  necessary.  Provision  is  made  in  the  frame  for  securing  the  lower 
sample  holder  and  for  the  vertical  and  horizontal  hydraulic  actuators 
used  for  application  of  normal  and  shear  loads.  A  winch  and  trolley 
system  is  installed  on  the  load  frame  to  aid  in  the  assembly  of  the 
specimen  in  the  shear  device. 
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Detail  of  Rail  System  for 
Guiding  Motion  of  Upper 
Sample 
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The  loading  system  consists  of  the  hydraulic  system  used  to 
apply  loads  including  the  vertical  and  horizontal  hydraulic  actuators 
and  all  attachments  necessary  to  secure  the  actuators  to  the  load  cells 
and  the  upper  sample  holder.  The  system  is  shown  schematically  in 
Fig.  3-9. 

At  the  heart  of  the  load  system  is  a  hydraulic  pump  capable  of 
providing  a  fluid  pressure  up  to  3000  psi  (21  MPa)  at  a  flow  rate  of 
25  gallons /min.  This  pump  services  a  servo-controlled  electro-hydraulic 
system,  to  be  discussed  subsequently,  which  controls  the  flow  of  hy¬ 
draulic  fluid  to  two  hydraulic  actuators.  One  actuator  is  mounted  to 
the  load  frame  such  that  it  reacts  against  a  vertical  member  of  the 
frame  and  provides  the  shear  force  transmitted  to  the  upper  sample. 

The  other  actuator  mounts  to  the  top  of  the  load  frame  in  the  center 
and  provides  the  normal  force  transmitted  to  the  upper  sample.  Both 
actuators  have  a  stroke  of  8  in.  (20  cm.).  The  horizontal  actuator  is 
capable  of  applying  a  force  of  12  kips  (53  KN)  and  the  vertical  actuator 
is  capable  of  applying  a  30  kip  (133  KN)  force. 

The  actuators  are  connected  to  the  load  frame  and  load  cells 
using  connections  and  fittings  as  shown  in  Fig.  3-9.  The  connections 
are  designed  to  provide  a  tight  fit  while  not  allowing  the  load  system 
to  inhibit  the  motion  of  the  upper  sample  box.  Both  the  horizontal  and 
vertical  actuators  are  pinned  to  the  load  frame  using  clevises  which 
fix  the  end  of  the  actuator  with  respect  to  translation  but  allow  rota¬ 
tion  about  any  axis.  Thus  the  ends  of  the  actuators  which  contact  the 


upper  sample  are  allowed  to  move  with  the  upper  sample  as  they  rotate 
about  fixed  ends.  The  horizontal  actuator  piston  threads  directly  into 
a  load  cell  which  is  mounted  to  a  circular  plate.  Another  clevis 
assembly  joins  the  circular  plate  to  the  upper  sample  holder  at  the 
center  of  the  front  face,  2  in.  (5.08  cm.)  up  from  the  bottom.  The 
vertical  actuator  piston  feeds  directly  into  a  clevis  assembly  that  is 
threaded  to  a  vertical  load  cell.  Both  clevis  assemblies  at  the  load 
cell  end  of  the  horizontal  and  vertical  actuators  contain  ball  bearings 
which  allow  the  ends  of  the  actuators  to  rotate  about  any  axis. 

Control  System 

The  hydraulic  actuators  are  controlled  with  a  closed- loop 
servo-controlled  electro-hydraulic  system.  This  system  consists  of  an 
MTS  model  #436  control  unit,  an  MTS  service  manifold,  two  MTS  model  # 406 
controllers  (one  for  each  actuator),  a  frequency  generator,  and  a  Moog 
model  #A076  servovalve  in  line  with  each  actuator.  The  system  is  shown 
schematically  in  Fig.  3-10. 

The  hydraulic  actuators  may  be  controlled  by  the  closed-loop 
electro-hydraulic  system  in  either  a  displacement-controlled  mode  or  a 
load-controlled  mode.  The  system  compares  a  command  signal  generated  by 
the  frequency  generator  with  a  feedback  signal  from  instruments  measuring 
loads  or  displacements  realized  by  the  test  arrangement.  The  command 
signal  describes  the  desired  motion  or  force  of  the  actuator  while  the 
feedback  signal  indicates  the  actual  motion  or  force  of  the  actuator.  A 
command  signal,  in  terms  of  amplitude  and  frequency  of  the  controlled 
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variable,  either  load  or  displacement,  is  sent  from  the  electronic  con¬ 
troller  to  the  servovalve  mounted  on  the  hydraulic  actuator.  A  feedback 
signal  indicating  the  value  of  the  controlled  variable  is  returned  to 
the  controller  and  compared  with  the  command  signal.  A  deviation  be¬ 
tween  the  command  signal  and  the  feedback  signal  results  in  an  error 
signal  that  is  returned  as  a  command  signal  to  reduce  the  error  to  zero, 
or  close  the  loop. 

Measurement  of  Test  Variables 

Measurements  made  during  testing  consist  of  the  vertical  and 
horizontal  forces  transmitted  to  the  upper  sample  via  the  hydraulic  ac¬ 
tuators,  as  well  as  the  horizontal  motion,  the  vertical  displacement, 
and  rotation  about  an  axis  in  the  plane  of  shear  perpendicular  to  the 
direction  of  shear  of  the  upper  sample.  Figure  3-11  indicates  the 
placement  of  instrumentation,  and  the  variables  measured. 

Vertical  and  Horizontal  Loads 

Vertical  loads  are  measured  with  a  Strainsert  load  cell, 

Model  //FFL  (45/±30)  U(C)-SPKT  which  has  a  capacity  of  30,000  lbs. 

(133  kN)  in  both  tension  and  compression.  Horizontal  loads  are  measured 
with  a  Strainsert  load  cell  Model  //FFL  (18/+12)  U-(3/±2)  which  has  a 
capacity  of  12,000  lbs.  (53  kN)  in  both  tension  and  compression.  Both 
load  cells  are  bolted  to  load  plates  and  measure  loads  through  the  use 


of  strain  rosettes. 


Left  Side  View 
(inductance  coils  not  shown) 


Right  Side  View 


1.  Vertical  load  cell 

2.  Horizontal  load  cell 

3.  Horizontal  displacement  LVDT 

4.  Inductance  coils 


5.  Inside  front  displacement  sensor 

6.  Inside  back  displacement  sensor 

7.  Outside  front  vertical  LVDT 

8.  Outside  back  vertical  LVDT 


Figure  3-11.  Placement  of  Instruments 


Horizontal  Displacement 

Horizontal  displacement  of  the  upper  sample  is  measured  using  a 
Scheavitz  linear  variable  differential  transformer  (LVDT) 

Model  #PCA-117-1000.  This  LVDT  is  spring  loaded  and  has  a  stroke  of 
plus  or  minus  1  in.  (2.54  cm.).  Attachment  of  the  LVDT  is  made  at  the 
side  of  the  bottom  sample  holder,  and  the  gage  head  of  the  LVDT  contacts 
a  plate  which  is  suspended  from  the  rail  that  guides  the  upper  sample. 
Locating  the  LVDT  in  this  manner  bypasses  measurement  of  displacements 
which  are  extraneous  to  the  motion  of  the  upper  sample.  These  extraneous 
displacements  include  some  small  amount  of  slip  occurring  in  the  clevises 
used  in  the  actuator  connections ,  and  any  motion  of  the  lower  sample 
holder. 

Vertical  Displacements 

Vertical  displacements  of  the  upper  sample  holder  are  made  at 
five  locations  as  shown  in  Fig.  3-11.  These  include  two  LVDT's  at  the 
far  corners  of  the  bottom  sample  which  are  attached  to  the  lower  sample 
holder  and  contact  the  rail  guiding  the  upper  sample  (outside  front  and 
back  vertical  LVDT's),  two  displacement  sensors  which  attach  to  plates 
protruding  from  the  side  of  the  front  sample  near  the  bottom  and  contact 
a  smooth  aluminum  plate  attached  to  the  cutout  portion  on  the  side  of  the 
lower  sample  holder  (inside  front  and  back  displacement  sensors) ,  and  a 
pair  of  inductance  coils  which  are  placed  in  cavities  formed  in  the 
center  of  the  top  and  bottom  samples.  Measuring  the  vertical  displace¬ 
ment  at  five  locations  allows  measurements  of  vertical  displacement  as 
well  as  rotations,  if  any,  of  the  upper  sample. 
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Outside  Front  and  Back  LVDT's.  The  outside  front  and  back 
LVDT's  are  Sheavitz  Model  #PCA-1  16-300.  This  model  is  externally  spring 
loaded  and  has  a  stroke  of  plus  or  minus  0.30  in.  (76.2  mm). 

Inside  Front  and  Back  Displacement  Sensors.  The  inside  front 
and  back  displacement  sensors  are  Model  #9610  Spring  Return  Linear  Posi¬ 
tion  Modules  purchased  from  Duncan  Electronics,  Costa  Mesa,  California. 
These  displacement  sensors  have  a  stroke  of  plus  or  minus  0.5  in. 

(1.27  cm.)  and  are  DC-operated  instruments  which  utilize  a  Wheatstone 
Bridge  arrangement.  Displacements  are  measured  via  a  change  in  resis¬ 
tance  induced  as  the  piston  of  the  instrument  extends  or  retracts. 

Inductance  Coils.  The  inductance  coils  used  to  measure  vertical 
displacements  near  the  center  of  the  joint  by  placing  them  in  cavities 
formed  at  the  center  of  the  upper  and  lower  structured  surfaces  were 
purchased  from  Bison  Instruments,  Minneapolis,  Minnesota.  The  coils 
come  as  part  of  a  system  which  includes  the  coils  (sensors)  and  an  ex¬ 
ternal  instrument  package  (Bison  Model  #41Q1A  Soil  Strain  Gage  Indicator). 
The  system  is  similar  in  concept  to  one  described  by  Selig  (1964)  ,  and  a 
full  description  can  be  found  in  the  Bison  Gage  Instruction  Manual. 

Strain  is  measured  as  a  change  in  spacing  between  two  reference  points 
divided  by  the  original  spacing  (gage  length) .  The  system  is  depicted 
schematically  in  Fig.  3-12  (a)  and  (b) . 

The  sensors,  each  of  which  is  a  disk-shaped  coil,  are  placed  in 
near  parallel  and  coaxial  alignment.  One  sensor  is  placed  in  the  cavity 
of  the  upper  sample,  and  the  other  is  placed  in  the  cavity  of  the  bottom 
sample.  They  are  separated  a  distance  over  which  the  strain  is  to  be 
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Figure  3.12  (a).  Pictorial  Diagram  of  Soil  Strain 

Gage  (Truesdale  and  Anderson,  1964) 
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Figure  3.12  (b) .  Block  Diagram  of  Soil  Strain 

Gage  (Truesdale  and  Anderson,  1964) 


averaged.  The  sensors  are  2.0  in.  (5.08  cm.)  in  diameter,  which  is  a 
standard  size  currently  manufactured.  One  inch  (2.54  cm.)  diameter 
coils  are  also  available.  The  separation  of  the  sensors  is  related  to 
the  electro-magnetic  coupling  between  the  two.  By  means  of  an  inductance 
bridge.  Fig.  3-12  (b) ,  an  output  voltage  as  a  function  of  strain  may  be 
obtained  since  a  change  in  spacing  from  the  initial  spacing  produces  an 
imbalance.  The  sensors  operate  at  any  spacing  between  one  and  four  times 
the  nominal  sensor  diameter.  The  effects  of  rotational  or  transverse 
movement,  to  be  discussed  later,  and  moisture  and  temperature  are  usually 
negligible.  The  proximity  of  metals  may  affect  the  output.  This  will  be 
discussed  in  Chapter  4,  Calibration  of  Instruments. 

The  Model  //4101A  Soil  Strain  Indicator  to  which  the  coils  are 
connected  contains  all  necessary  driving,  amplification,  balancing, 
read-out  and  calibration  controls  and  a  self-contained  power  supply. 

Fig.  3-13.  The  bridge  balance  is  accomplished  by  means  of  the  phase  and 
amplitude  controls  using  the  meter  to  indicate  null.  The  amplitude  dial 
reading  corresponds  to  the  sensor  spacing.  Changes  in  spacing  may  be 
determined  by  meter  deflection  from  zero  or  by  voltage  output  recorded  by 
a  data  acquisition  system  connected  to  the  rear  panel.  The  calibration 
control  is  used  to  set  up  the  output  sensitivity  so  that  it  corresponds 
to  the  desired  amount  of  strain,  about  0.4  percent  in  this  research. 

This  particular  instrument  was  designed  for  dynamic  strain  measurement 
and  it  has  the  ability  to  sense  changes  as  fast  as  0.25  microseconds. 


Data  Acquisition  System 


A  Digital  Electronic  Corporation  M1NC-11  (Micro  Instrumentation 
Computer)  was  used  for  data  acquisition  for  tests  with  the  CYMDOF  device. 
The  MINC-11  consists  of  a  64  K  RAM  microcomputer  system  designed  for 
laboratory  data  acquisition  and  data  reduction.  The  MINC  used  in  con¬ 
junction  with  the  CYMDOF  device  is  equipped  with  two  preamplifiers  and 
an  analog  to  digital  (A/D)  converter.  The  preamplifiers  step  up  or  down 
all  voltage  signals  so  that  they  are  between  I  and  5  volts  and  the  A/D 
converter  allows  for  analog  data  such  as  the  voltage  output  from  load 
cells,  LVDT's,  displacement  sensors,  and  inductance  coils  to  be  con¬ 
verted  into  digital  form  for  storage  and  manipulation  by  the  computer. 

The  system  may  accommodate  up  to  eight  preamplified  channels,  and  up  to 
64  channels  with  no  amplification.  The  accuracy  of  readings  using  the 
preamplifier  is  plus  or  minus  0.005  volts.  The  system  is  capable  of 
reading  a  channel  every  microsecond,  and  the  A/D  converter  takes  about 
40  microseconds  to  process  the  signal. 

The  MINC  utilizes  8.0  in.  (20.3  cm.)  floppy  magnetic  disks  for 
data  and  program  storage,  and  both  BASIC  and  FORTRAN  programming 
languages  can  be  used. 

The  operation  of  the  MINC  and  the  data  acquisition  program  are 
fully  described  elsewhere  (Drumm  and  Desai,  1983).  The  major  features 
of  the  data  acquisition  and  data  reduction  programs  will  be  briefly  de¬ 


scribed  here. 
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Acquisition  of  data  is  triggered  by  an  internal  time  clock, 
which  is  also  a  part  of  the  MINC  system.  During  a  slow  (quasi-static) 
cyclic  test  data  from  all  the  instruments  are  sampled  every  three 
seconds.  During  a  fast  cyclic  test  readings  of  all  instruments  are 
taken  continuously  for  cycles  1  through  10,  then  continuously  for  two 
cycles  after  cycles  number  20,  50,  100,  150,  200,  250,  300,  350,  400, 

450  and  500.  Data  from  all  the  instruments  may  be  sampled  up  to  sixteen 
times  during  a  given  fast  cycle.  The  voltage  signal  from  each  instrument 
contains  some  "noise"  or  high-frequency  variation  in  voltage.  Since 
these  deviations  in  the  voltage  signal  may  be  recorded  by  the  MINC,  two 
voltage  readings  are  taken  for  each  sampling  and  averaged  during  the  data 
reduction  process.  Thus,  during  a  fast  cyclic  test,  for  each  load  cycle 
for  which  data  is  recorded,  thirty-two  actual  voltage  readings  per  cycle 
are  taken  and  stored  in  the  memory  of  the  microcomputer  for  the  duration 
of  the  test.  Upon  completion  of  the  test,  the  computer  writes  all  the 
test  data  onto  a  floppy  disk  for  permanent  storage.  The  voltages  may 
then  be  printed  out  to  obtain  a  hard  copy  of  the  raw  data,  or  the  data 
may  be  submitted  to  a  data  reduction  program. 

The  data  reduction  program  computes  an  arithmetic  average  of  the 
two  data  readings  taken  at  each  sampling  to  yield  one  voltage  sample. 

From  the  average  voltage  values,  the  shear  stress,  relative  tangential 
displacement,  normal  stress,  and  vertical  relative  displacement  are  de¬ 
termined.  Stresses  are  calculated  by  dividing  loads  by  the  area  of  the 


rock  joint  (idealized)  surface,  which  gives  an  average  value  for  stress. 
The  maximum  and  minimum  shear  stress  and  the  maximum  and  minimum  rela¬ 
tive  displacement  within  a  given  cycle  are  also  determined. 

The  results  of  the  data  reduction  program  may  be  printed  out  to 
yield  a  listing  of  all  the  computed  stresses  and  relative  displacements 
for  each  cycle  for  which  readings  were  taken.  The  maximum  and  minimum 
shear  stresses  and  relative  tangential  displacements  are  also  output. 

Assembly  of  the  Device 

After  securing  the  upper  and  lower  samples  to  their  respective 
holders,  the  structured  surfaces  of  the  samples  are  cleaned  and  the  top 
sample  placed  on  top  of  the  bottom  sample  in  the  load  frame  with  a  winch 
and  trolley  system.  By  suspending  the  upper  sample  slightly  above  the 
lower  sample,  the  rail  system  can  be  installed  and  the  horizontal  actua¬ 
tor  connected  to  the  upper  sample  holder  without  disturbing  the 
structured  surface.  Before  lowering  the  upper  sample  completely,  the 
centers  of  the  top  and  bottom  samples  are  aligned  with  the  horizontal 
actuator  using  center  marks  placed  on  the  samples  before  they  are  secured 
to  their  sample  holders.  Aligning  the  samples  in  this  manner  assures  as 
close  a  mating  as  possible  between  the  upper  and  lower  structured 
surfaces.  Also,  the  Bison  gages  placed  in  the  center  of  the  top  and 
bottom  samples  need  to  be  coaxially  aligned  as  close  as  possible  prior 
to  testing.  After  lowering  the  upper  sample  to  its  proper  position,  the 


horizontal  actuator  is  leveled. 
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The  vertical  load  system  is  installed  such  that  the  load  plate 
is  placed  level  on  top  of  the  upper  sample  and  the  vertical  load  actua¬ 
tor  is  directly  above  the  center  of  the  samples  prior  to  testing.  In 
order  to  accomplish  this,  a  leveling  plate  is  placed  atop  the  upper 
sample,  and  leveled.  The  vertical  load  plate  with  the  load  cell  is 
placed  on  this  level  plate.  The  vertical  actuator  is  then  trolleyed 
into  position  above  the  center  of  the  samples  and  secured.  Vertical 
loads  can  now  be  applied  to  the  sample  by  lowering  the  piston  of  the  ac¬ 
tuator  into  a  yoke  which  is  connected  to  the  load  cell. 

Once  the  device  is  assembled,  electrical  connections  are  made 
between  the  measurement  instruments,  the  control  units,  readout  devices, 
and  any  excitation  or  signal  conditioners  required  by  the  instruments. 

The  vertical  and  horizontal  load  cells,  the  horizontal  LVDT ,  and  the  two 
(front  and  back)  outside  LVDT's  are  connected  to  MTS  406  controllers. 
These  controllers  provide  excitation  and  signal  conditioning  required  by 
the  instruments  and  can  monitor  device  outputs  to  be  used  for  motion  or 
force  control.  The  control  units  can  also  direct  the  output  signals  to 
the  MINC-1I  data  acquisition  system. 

The  Bison  gages  used  to  measure  vertical  displacements  at  the 
center  of  the  joint  are  connected  to  a  Bison  Model  4101A  Soil  Strain 
Indicator.  This  instrument  provides  signal  conditioning  and  excitation 
to  the  coils  and  directs  the  output  signals  to  the  MINC-1I  data  acqui¬ 


sition  system. 
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The  two  (front  and  back)  linear  motion  sensors  placed  on  the 
protruding  plates  of  the  upper  sample  receive  excitation  of  10  volts 
from  a  DC  power  supply.  Sensor  outputs  are  wired  directly  to  the 
MINC-11  data  acquisition  system. 


CHAPTER  4 


CALIBRATION  OF  INSTRUMENTS 

The  instruments  employed  including  load  cells,  LVDTs,  displace¬ 
ment  sensors,  and  inductance  coils  respond  to  a  change  in  the  quantity 
to  be  measured,  be  it  load  or  displacement,  with  a  change  in  voltage 
output.  Calibration  of  instruments  involves  the  establishment  of  a 
definite  relationship  between  the  voltage  output  of  the  instrument  and 
the  quantity  measured.  The  recorded  voltage  output  of  the  instruments 
may  also  be  affected  by  factors  other  than  the  quantity  of  interest. 
These  factors  include  changes  in  temperature,  the  proximity  of  metals  or 
magnetic  fields,  capacitance  of  cables  used  in  making  connections,  the 
excitation  provided  to  the  instruments,  signal  conditioning,  and  charac¬ 
teristics  of  the  data  acquisition  system.  These  effects  must  either  be 
shown  to  be  negligible,  compensated  for,  or  controlled  in  some  manner 
such  that  the  change  in  voltage  due  to  the  measured  quantity  alone  is 
essentially  realized.  The  degree  to  which  this  may  be  accomplished,  and 
how  well  a  definite  relationship  between  the  measured  quantity  and  the 
voltage  output  is  established  will  determine  the  accuracy  and  the  re¬ 
peatability  of  the  instrument. 

Knowledge  of  the  accuracy  and  repeatability  of  instruments  is 
especially  important  when  studying  the  shear  behavior  of  rock  joints, 
particularly  with  regard  to  displacements.  The  tangential  displacement 
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which  occurs  prior  to  reaching  the  proportional  limit,  and  the  amount 
of  compression  which  occurs  prior  to  dilation  during  direct  shear 
testing  of  rock  joints  are  extremely  small  quantities.  The  constitutive 
model  requires  determination  of  certain  parameters  such  as  the  initial 
shear  stiffness,  and  parameters  which  are  a  part  of  the  yield  and  po¬ 
tential  functions  employed  in  the  plasticity  formulation,  which  are 
determined  using  these  small  quantities.  How  accurately  these  parame¬ 
ters  are  evaluated  is  directly  related  to  the  accuracy  and  repeatability 
of  the  instruments  used  to  make  measurements  during  testing. 

Details  of  instrument  calibrations  are  included  in  part  II  of 
the  final  report  for  this  project  submitted  to  the  Air  Force  Office  of 
Scientific  Research.  This  report  is  entitled  "Constitutive  Modeling  of 
Idealized  Rock  Joints  Under  Cyclic  Loading;  Final  Report,  Part  II." 

Only  calibration  of  inductance  coils  will  be  discussed  in  this  chapter. 

As  justification  for  including  this  chapter,  results  from  in¬ 
ductance  coil  measurements  are  shown  to  be  reliable.  To  the  author's 
knowledge,  this  is  the  first  time  these  instruments  have  been  used  for 
the  measurement  of  vertical  displacement  during  direct  shear  testing. 
Proper  calibration  procedures  must  be  established  as  well  as  the  pre¬ 
cision,  accuracy  and  repeatability  of  the  instrument. 

Inductance  Coils 

The  Bison  Instrument  Inductance  Coils  (Bison  Gages)  were  cali¬ 
brated  using  the  calibration  fixture  depicted  in  Fig.  4-1.  The 


Figure  4-1.  Calibration  Instrument 
for  Inductance  Coils 
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calibration  fixture  is  made  of  wood  such  that  the  fixture  does  not 
affect  the  magnetic  field  surrounding  the  coils.  Coils  are  mounted  to 
the  calibration  fixture  such  that  they  are  coaxially  aligned.  One  coil  is 
held  fixed  and  the  other  coil,  attached  to  a  sliding  wooden  rod,  may  be 
translated  by  turning  a  micrometer  at  the  end  of  the  fixture.  The 
micrometer  is  manufactured  by  Schaevitz  and  has  a  total  throw  of  1.0  in. 
(2.54  cm.)  and  an  accuracy  of  0.0001  in.  (0.0025  mm.).  During  calibra¬ 
tion  the  coils  are  connected  to  the  Bison  Model  #4101A  Soil  Strain 
Indicator  which  supplies  excitation  and  signal  conditioning  required  by 
the  coils,  and  directs  output  to  the  MINC-11  data  acquisition  system. 

The  calibration  of  the  Bison  gages  involves  two  phases  of  cali¬ 
bration,  spacing  calibration  and  meter  calibration.  Both  are  fully 
described  in  the  "Bison  Instruments  Soil  Strain  Gage  Model  4101A  Instruc¬ 
tion  Manual."  The  spacing  calibration  consists  of  mounting  the  coils  on 
the  calibration  fixture  and  determining  settings  (amplitude  dial)  on  the 
indicator  unit  required  to  balance  the  instrument  (inductance  bridge)  at 
different  coil  spacings.  The  meter  calibration  is  intended  to  determine 
the  displacement  that  will  cause  the  meter  to  deflect  100%  (5.0  volts) 
from  an  initially  balanced  condition  which  is  termed  as  100%  full  scale 
displacement.  The  resulting  calibration  allows  determination  of  the 
initial  spacing  of  the  coils  placed  in  the  test  specimens  (idealized 
rock  joints) ,  and  accumulated  deformation  from  the  voltage  output  (meter 
deflection)  at  a  given  time.  The  displacement  is  calculated  by  multi¬ 
plying  the  per  cent  meter  deflection  by  the  100%  full  scale  deflections 
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determined  from  the  calibration.  Typical  calibration  plots  for  both 
the  spacing  and  meter  calibrations  are  presented  in  Fig.  4-2  (a)  and  (b) . 

The  procedure  for  using  these  plots  for  measurements  during 
testing  is  as  follows: 

-  After  assembly  of  the  device,  the  strain  indicator  is  balanced. 
The  initial  coil  spacing  is  determined  from  Fig.  4-2  (a)  corresponding 

to  the  amplitude  dial  setting  shown  on  the  indicator  unit.  For  example, 
if  the  amplitude  reading  is  900,  the  initial  spacing  is  3.35  in. 

(8.51  cm.). 

-  When  the  coils  move  during  testing,  the  indicator  shows  a 
positive  (extension)  or  negative  (compression)  reading.  The  ordinate 
in  Fig.  4-2  (b)  shows  100%  full  scale  displacement.  Then,  depending 
upon  the  meter  deflection  (some  percentage  of  5.0  volts),  displacements 
are  calculated  as  a  percentage  of  the  full  scale  displacement  for  the 
given  initial  spacing.  For  example,  if  the  meter  reading  is  x,  the 
measured  displacement  for  a  spacing  of  3.35  in.  would  be  (0.013)  (x)/100, 
where  0.013  corresponds  to  an  amplitude  dial  setting  of  900. 

Linearity  of  Meter  Deflection 

The  calibration  assumes  the  deflection  to  be  a  linear  function 
of  coil  displacement.  This  was  verified  by  setting  a  pair  of  coils  on 
the  calibration  fixture,  and  by  comparing  voltage  outputs  with  applied 
displacements.  This  verification  was  found  to  be  satisfactory,  particu¬ 
larly  for  displacements  above  0.0008  in.  (0.02  mm.).  For  displacements 
as  small  as  0.0003  in.  (.0076  mm.),  the  error  is  as  high  as  40%. 


Table  4-1  shows  details  of  this  verification. 


SETS 


78 


Precision 

The  precision  of  the  Bison  gages  was  checked  by  repeatedly 
mounting  the  same  pair  of  coils  on  the  calibration  fixture  and  comparing 
results.  Table  4-2  summarizes  these  results.  The  precision  of  the  in¬ 
strument  is  within  0.1%,  based  on  a  full  scale  output  of  5  volts. 

Effects  on  Calibration 

The  voltagf'  output  of  the  Bison  gages  may  be  affected  by  factors 
such  as  misalignment  of  the  coils,  rotations  of  the  coils  during  deforma¬ 
tion,  and  nearby  metals  or  magnetic  fields. 

Misalignment.  Misalignment  of  the  Bison  gages,  as  depicted  in 
Fig.  4-3,  is  a  shift  in  the  relative  position  of  the  coils  such  that  they 
are  no  longer  coaxial.  The  effect  of  coil  misalignment  on  the  calibra¬ 
tion  of  the  Bison  gages  is  of  interest  for  two  reasons.  First,  the  coils 
may  not  be  perfectly  coaxial  when  placed  within  the  specimens  due  to 
errors  in  measurement.  Secondly,  as  the  upper  sample  is  displaced  during 
direct  shear,  the  coils  displace  with  respect  to  one  another  and  become 
misaligned.  The  calibration  assumes  that  the  coils  are  initially  co¬ 
axially  aligned.  The  purpose  for  studying  the  effect  of  coil 
misalignment  is  twofold: 

1.  To  determine  the  accuracy  of  the  Bison  gage  readings,  given 
the  uncertainty  involved  in  their  initial  position;  i.e.,  misalignment 
due  to  errors  in  measurement  during  placement  within  shear  test  specimens. 

2.  To  determine  a  correction  factor  to  be  applied  to  readings 
after  a  known  amount  of  misalignment  has  taken,  place;  i.e.,  misalignment 
due  to  tangential  displacement  imposed  during  direct  shear  testing. 


Table  4-2.  Verification  of  Precision  of  Bison  Soil  Strain  Gage 

Amplitude  Setting  =  900 

Reading  Initial  Spacing  Strain  Z  F.S.  Strain  %  F.  S. 


# 

(in.) 

Compression 

Extension 

1 

4.504 

-5.45 

7.04 

2 

4.504 

-5.56 

6.98 

3 

4.507 

-5.49 

7.03 

4 

4.502 

-5.55 

7.02 

5 

4.513 

-5.51 

7.04 

F.S.  =  Full  Scale 
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The  effect  of  coil  misalignment  on  the  calibration  was  studied 
using  a  special  vertical  offset  caliper,  fastened  to  the  calibration 
fixture  as  shown  in  Fig  'i-4.  One  Bison  gage  is  secured  to  the  calibra¬ 
tion  fixture  as  before.  The  other  is  fastened  to  the  vertical  offset 
caliper.  The  vertical  offset  caliper  enables  one  Bison  gage  to  move  up 
or  down  vertically  in  measured  increments  with  respect  to  the  other, 
initially  coaxially  aligned,  coil.  The  Bison  gage  calibration  was  per¬ 
formed  with  the  gages  coaxially  aligned,  misaligned  by  0.1  in.  (2.54  mm.), 
and  misaligned  by  0.5  in.  (1.27  cm.). 

The  meter  calibration  at  different  offsets  was  performed  con¬ 
sidering  two  different  cases.  In  one  case,  the  Bison  gages  were  balanced 
in  the  coaxially  aligned  configuration,  and  then  offset.  This  case 
corresponds  to  direct  shear  test  conditions  where  the  coils  may  be  ini¬ 
tially  coaxially  aligned  when  the  instrument  is  balanced,  and  then 
misaligned  due  to  tangential  displacement  during  shear.  For  the  second 
case,  the  Bison  gages  were  balanced  in  the  offset  position,  and  this  case 
was  used  to  ascertain  any  errors  which  may  be  present  due  to  uncertainties 
involved  in  the  initial  alignment  of  the  coils. 

Figures  4-5  (a)  and  (b)  present  results  from  the  meter  and 
spacing  calibrations  at  the  various  degrees  of  alignment.  Figure  4-5  (a) 
indicates  that  the  spacing  calibration  is  not  significantly  affected  by 
offsets  as  high  as  0.5  in.  (1.27  cm.).  The  difference  between  the 
spacing  calibration  for  the  coaxially  aligned  case  and  the  case  of  0.1  in. 
(1.27  cm.)  offset  differs  from  the  others  by  1%  at  a  spacing  of  around 
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3.8  in.  (9.65  cm.). 


Figure  4-4.  Vertical  Offset  Caliper 


Figure  4-5(a).  Effect  of  Coil  Misalignment  on  Spacing  Calibration 
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Figure  4-5 (b)  indicates  that  when  the  Bison  gages  are  balanced 
in  the  initial  coaxially  aligned  position  and  then  offset  by  0.10  in. 
(2.54  cm.),  the  calibration  is  shifted  by  0.005  in.  (0.127  mm.)  from 
the  coaxially  aligned  case.  The  shift  is  such  that  during  extension 
(positive  displacement) ,  the  100%  full  scale  displacement  at  a  given 
amplitude  dial  setting  is  decreased  by  0.005  in.  (0.127  mm.)  above  the 
100%  full  scale  displacement  at  the  same  amplitude  dial  setting  for  the 
case  when  the  coils  remain  in  the  coaxial  alignment.  During  compression 
(negative  displacement)  the  absolute  value  of  the  100%  full  scale  dis¬ 
placement  at  a  given  amplitude  dial  setting  when  the  coils  are  offset 
0.10  in.  (2.54  mm.)  is  increased  0.005  in.  (0.127  cm.)  from  the  coaxially 
aligned  case.  The  shift  in  the  calibration  curve  for  an  offset  of 
0.5  in.  (1.27  cm.)  is  not  presented  since  this  shift  was  so  large  that 
the  needle  of  the  voltage  meter  displaced  the  full  amount  due  to  the 
offset  alone. 

Using  this  data,  an  adjustment  to  the  meter  calibration  was  de¬ 
termined  which  could  account  for  the  effects  of  coil  misalignment, 
given  the  misalignment  occurred  with  respect  to  an  assumed  initial  co¬ 
axially  aligned  configuration.  The  adjustment  assumes  that  the  shift  in 
the  calibration  curve  between  the  case  when  the  coils  remain  in  the  co¬ 
axial  alignment  and  when  they  are  misaligned  is  linear  with  respect  to 
misalignment  for  misalignments  between  0.0  and  0.1  in.  (0.0  to  2.54  cm.). 
The  adjustment  was  incorporated  into  the  calibration  as  follows: 
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1.  After  assembly  of  the  shear  device,  as  discussed  in 
Chapter  3,  the  Bison  gages  are  initially  balanced. 

2.  Using  the  amplitude  dial  reading,  and  assuming  the  coils  are 
initially  in  coaxial  alignment,  the  100  per  cent  full-scale  displacement 
is  determined  from  the  meter  calibration. 

3.  During  a  direct  shear  test,  the  tangential  displacement  and 
the  voltage  output  (meter  deflection)  of  the  Bison  indicator  unit  are 
determined. 

4.  The  100  per  cent  full-scale  displacement  determined  in 
Step  2  is  multiplied  by  the  per  cent  meter  deflection  determined  in 
Step  3  to  yield  the  unadjusted  vertical  displacement. 

5.  An  adjustment  is  calculated  based  on  the  tangential  displace¬ 
ment  (coil  misalignment)  which  corrects  the  vertical  displacement 
determined  in  Step  4.  The  adjustment  is  equal  to  the  ratio  of  the  tan¬ 
gential  displacement,  determined  in  Step  2,  to  0.1  in.  (2.54  mm.), 
multiplied  by  0.005  in.  (0.127  ram.).  This  adjustment  is  only  valid  for 
offsets  between  0.0  and  0.1  in.  (0.00  to  2.54  mm.). 

Table  3  presents  results  from  the  meter  calibration  of  the  Bison 
gages  when  misaligned  by  0.10  in.  (2.54  mm.),  and  balanced  in  this  off¬ 
set  configuration,  and  compares  these  results  with  the  meter  calibration 
of  the  coils  when  coaxially  aligned.  The  table  indicates  that  when  the 
Bison  gages  are  initially  balanced  in  the  offset  configuration,  the 
difference  in  100  per  cent  full-scale  displacement  for  a  given  amplitude 
dial  setting  from  the  coaxially  aligned  case  is  within  0.001  in. 
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Table  4-3.  Comparison  of  Meter  Calibration 

With  Coils  Aligned  and  Misaligned 


Full  Scale  Displacement 


Offset 


Amplitude 

Setting 

Exten. 

XU. 

Comp. 

.  L 

Exten. 

XU  • 

Comp. 

• 

800 

0.009 

0.010 

0.010 

0.009 

850 

0.010 

0.011 

0.011 

0.010 

900 

0.014 

0.013 

0.014 

0.013 

• 

950 

0.018 

0.014 

0.018 

0.015 

1000 

0.019 

0.015 

0.020 

0.020 

i 


(.025  mm.).  Therefore,  the  uncertainty  involved  with  the  initial  posi¬ 
tion  of  the  coils  relative  to  one  another  is  not  considered  significant 
for  offsets  less  than  0.10  in.  (2.54  mm.). 

Rotations.  During  direct  shear  testing,  not  only  does  the  upper 
sample  translate  with  respect  to  the  lower  sample,  but  a  (small) 
amount  of  rotation  of  the  upper  sample  about  an  axis  in  the  plane  of 
shear  perpendicular  to  the  direction  of  shear  may  also  be  present. 
Effects  of  coil  rotations  on  the  Bison  gage  calibration  are  harder  to 
evaluate  than  the  effects  of  coil  misalignment  because  of  the  difficulty 
in  producing  rotation  without  axial  movement.  Truesdale  and  Schwab 
(1965)  found  that  for  angular  rotation  less  than  plus  or  minus  10  de¬ 
grees,  the  effect  was  not  significant.  Angular  rotation  in  excess  of 
this  amount  is  not  anticipated  during  direct  shear  testing. 

Surrounding  Metal.  With  the  Bison  gages  placed  within  the 
samples,  and  the  shear  device  assembled  as  described  in  Chapter  3,  the 
coils  are  surrounded  by  the  metal  sample  holders  as  shown  in  Fig.  4-6. 
The  effect  of  surrounding  metal  on  the  calibration  of  the  Bison  gages 
was  checked  by  placing  the  gages,  mounted  to  the  calibration  fixture, 
within  the  upper  sample  holder  as  depicted  in  Fig.  4-7.  Since  the 
upper  sample  holder  has  inside  dimensions  of  12  in.  by  12  in.  (30.5  cm. 
by  30.5  cm.),  the  boundaries  of  the  metal  box  were  6  in.  (15.25  cm.) 
from  the  center  of  the  inductance  coils,  or  about  three  coil  diameters 
away.  With  the  coils  in  this  configuration,  a  spacing  and  meter  cali¬ 
bration  was  performed  with  the  instrument.  Table  4-4  presents  results 


Bison  Gages  Su 
Sample  Holder 


of  the  calibration  performed  in  the  presence  of  the  surrounding  metal, 
and  compares  them  to  the  results  of  the  calibration  performed  in  the  ab¬ 
sence  of  the  surrounding  metal. 

The  coil  spacing  at  a  given  amplitude  dial  setting  in  the 
presence  of  surrounding  metal  differed  by  less  than  5  per  cent  from  the 
coil  spacing  at  the  same  amplitude  dial  setting,  but  in  the  absence  of 
surrounding  metal.  Results  from  the  meter  calibrations  indicate  that 
the  100  per  cent  full-scale  displacements  at  a  given  amplitude  dial 
setting  in  the  presence  of  surrounding  metal  differed  between  6  per  cent 
to  15  per  cent  from  the  100  per  cent  full-scale  displacements  at  the 
same  amplitude  dial  setting,  but  in  the  absence  of  surrounding  metal. 
This  difference  is  less  at  higher  amplitude  dial  settings.  The  errors 
in  the  calibration  of  the  Bison  gages  due  to  the  presence  of  surrounding 
metal  appear  to  be  acceptable,  and  the  calibrations  performed  in  the  ab¬ 
sence  of  surrounding  metal  were  used. 


CHAPTER  5 


DESCRIPTION  OF  SAMPLES 

Ideally,  laboratory  testing  intended  to  study  rock  joint  behavior 
should  be  performed  on  rock  joint  samples  obtained  from  the  field.  Due 
to  the  size,  number,  and  consistency  of  samples  required  for  this  study, 
obtaining  samples  from  the  field  was  neither  practical  nor  economical. 
Artificial  rock  joint  surfaces  were  prepared  by  casting  concrete  against 
a  structured  surface.  The  structured  surface  had  a  "sawtooth"  shaped 
pattern  intended  to  simulate  asperities  on  the  rock  joint  surface. 

Figure  5-1  depicts  the  "sawtooth"  shaped  asperity  pattern.  Idealized 
rock  joints  were  constructed  by  casting  top  and  bottom  samples  such  that 
when  placed  together,  the  asperities  from  the  top  and  bottom  produce  a 
closely  mated  joint,  as  shown  in  Fig.  5-1.  Different  rock  joint  rough¬ 
nesses  were  simulated  by  varying  the  angle  of  the  "sawtooth"  asperity 
side  with  respect  to  the  horizontal.  Samples  were  cast  having  asperity 
angles  of  0,  5,  7,  9  and  15  degrees;  the  first  denotes  "smooth"  surface. 

Forms 

Forms  for  casting  concrete  samples  were  constructed,  as  shown  in 
Fig.  5-2.  Bottom  sample  forms  had  inside  dimensions  of  (L  x  W  x  H) 

15.5  in.  x  15.5  in.  x  5.5  in.  (39  cm.  x  39  cm.  x  14  cm.).  Top  sample 
forms  had  inside  dimensions  of  11.5  in.  x  11.5  in.  x  5.5  in.  (29  cm.  x 
29  cm.  x  14  cm.).  Chamfer  strips  were  placed  around  the  perimeter  of  the 
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Figure  5-2.  Forms  for  Casting  Concrete  Samples 
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bottom  of  each  form  to  reduce  stress  concentrations  on  the  leading 
edges  of  the  samples  during  shear.  The  form  bottom  for  each  sample  con¬ 
sisted  of  the  structured  surface. 


!• 


4 


t 


Provisions  were  made  in  the  formwork  such  that  instrumentation 
could  be  placed  on  the  specimens.  Cylindrical  cavities  were  formed  in 
the  center  of  the  top  and  bottom  samples  so  that  inductance  coils  could 
be  placed  for  measurement  of  normal  displacement  close  to  the  joint 
surface.  These  cavities  were  formed  so  that  one-half  inch  (1.27  cm.)  of 
concrete  covered  the  bottom  of  the  cavity.  Also,  metal  plates  were 
attached  front  and  back  on  the  bottom  of  a  side  form  for  the  top  sample. 
When  partially  embedded  in  concrete,  these  plates  allowed  attachment  of 
linear  voltage  sensors  to  the  top  sample.  Figure  5-3 (a)  and  (b)  is  a 
sketch  of  the  sample  geometry. 

Method  of  Casting 

Using  six  sets  of  forms,  as  described  in  the  preceding  section, 
six  specimens  were  cast  at  a  time.  Initially,  six  bottom  specimens  were 
cast,  as  shown  in  Fig.  5-4.  After  one  day  of  curing,  samples  were  re¬ 
moved  from  the  forms  and  overturned.  A  piece  of  corrugated  sheet  metal, 
previously  used  as  the  structured  surface  at  the  form  bottom,  was  then 
placed  on  the  top  of  each  sample,  covering  the  newly  exposed  surfaces, 
which  exhibited  "sawtooth"  shaped  asperities.  The  top  sample  was  then 
cast  against  the  structured  surface  of  the  bottom  sample.  Fig.  5-5.  In 
this  manner,  a  close  fit  between  the  asperities  of  the  top  and  bottom 


I 


samples  was  assured. 


Top  View 


Front  View 


1" 


2.54 


Figure  5-3(b).  Sketch  of  Bottom  Sample  Geometry 


Figure  5-4.  Casting  of  Bottom  Samples 
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Figure  5-5.  Casting  of  Top  Samples 
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A  concrete  and  a  grout  mix  were  employed.  Details  of  the  mixes 
are  provided  in  the  next  section.  The  grout  mix  was  placed  for  the 
first  one-half  inch  (1.27  cm.)  near  the  joint  surface  for  both  the  top 
and  bottom  samples.  This  was  to  ensure  that  none  of  the  larger  aggre¬ 
gates  used  in  the  concrete  mix  were  near  the  joint  surface  since  these 
might  dislodge  during  testing  and  affect  the  results. 

Compaction  of  grout  and  concrete  during  casting  was  accomplished 
with  a  vibration  table  and  a  hand  held  concrete  vibrator.  Grout  placed 
near  the  structured  surface  of  the  bottom  samples  was  compacted  with  the 
vibration  table.  Fig.  5-6.  The  concrete  used  in  the  rest  of  the  sample 
was  compacted  with  the  hand  held  concrete  vibrator.  Compaction  of  grout 
placed  near  the  structured  surface  of  the  top  sample  was  accomplished  by 
touching  the  hand  held  concrete  vibrator  to  the  corrugated  sheet  metal. 
The  concrete  used  in  the  rest  of  the  sample  was  compacted  with  the  hand 
held  concrete  vibrator  as  before. 

Samples  were  cured  for  one  day  after  casting  the  top  samples, 
and  forms  were  removed.  Samples  were  then  placed  in  a  moist  room  to 
cure  for  28  days.  Upon  removal  from  the  curing  room,  the  joint  struc¬ 
tured  surfaces  were  cleaned  with  a  sponge  and  a  wire  brush  to  remove 
debris  or  residue  from  the  surfaces.  Figure  5-7  is  a  picture  of  a 
typical  sample  after  cleaning. 

Concrete  and  Grout  Mixes 

Table  5-1  shows  the  proportions  of  materials  used  in  the  concrete 
and  grout  mixes.  Figure  5-8  shows  gradation  curves  for  the  sand  and 


Table  5-1.  Concrete  and  Grout  Mixes 


(a)  Cement  Mix 

cement:  sand:  coarse  agg. :  water 
1:  3.26:  2.9:  0.75 

(b)  Grout  Mix 


cement:  sand:  water 


Gradation  Curve  for  Sand  and  Coarse  Aggregat 
Used  in  Concrete  Mix  Desisn 
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aggregates.  The  concrete  mix  used  was  the  same  as  that  used  by 
Salami  (1985) ,  so  that  the  model  he  developed  for  solid  concrete  could 
be  used  with  the  model  developed  here  for  joints  to  characterize  an 
idealized  jointed  rock  mass  system  at  a  later  time. 

Both  the  concrete  and  grout  mixes  were  designed  for  a  28-day 
unconfined  compressive  strength  of  around  2500  psi  (17.25  MPa). 

Cylinders  for  testing  the  28-day  unconfined  compressive  strength  of  both 
concrete  and  grout  mixes  were  prepared  during  sample  casting.  Each 
strength  determination  involved  testing  three  such  samples.  The  results 
of  the  three  tests  were  averaged  to  yield  the  28-day  unconfined  com¬ 
pressive  strength  for  a  given  mix.  Displacement  readings  were  taken 
during  testing  of  one  of  the  three  cylinders  so  that  a  determination  of 
Young's  modulus  after  28  days  for  a  given  mix  can  be  made.  Table  5-2 (a) 
and  (b)  presents  28-day  unconfined  compression  test  results  for  grout  and 
concrete  mixes.  Figures  9(a)  and  (b)  present  stress-strain  data.  An 
average  28-day  unconfined  compressive  strength  of  2500  psi  (17.25  MPa) 
prevailed  for  both  concrete  and  grout  mixes.  Young's  modulus  for  the 
concrete  and  grout  was  determined  as  3,500,000  psi  and  2,000,000  psi 
(24,000  and  13,780  MPa),  respectively. 
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Table  5-2 (a). 


Twenty-Eight  Day  Unconfined  Compressive 
Strength  Test  Results  for  Concrete  Mixes 


Date  of  Cast 

Cylinder 

f'c  (psi) 

f'c  (ave  for  cast) 

4/26/84 

1 

2759 

2 

2936 

3 

3148 

2948 

5/10/84 

1 

2511 

2 

2405 

3 

2830 

2582 

6/19/84 

1 

2582 

2 

2405 

3 

2830 

2511 

10/25/85 

1 

1875 

2 

1981 

3 

1946 

1934 

11/8/85 

1 

1592 

2 

1592 

3 

1592 

1592 

2/6/86 

1 

3325 

2 

3272 

3 

1592 

3278 

6/9/86 

1 

1627 

2 

2476 

3 

1839 

1981 

7/2/86 

1 

3537 

2 

3537 

3 

3077 

3384 

1  psi  -  6.89  KPa 


average  =  2526  psi 
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Table  5-2 (b) .  Twenty-Eight  Day  Unconfined  Compressive 


• 

Strength  Test  Results  for  Grout  Mixes 

Date  of  Cast 

Cylinder  f/c  (psi)  f'c  (ave.  for  cast) 

2/11/86 

1 

1704 

2 

1627 

3 

1450 

1627 

6/13/86 

1 

2724 

2 

1733 

• 

3 

1839 

2100 

7/9/86 

1 

3962 

2 

3750 

3 

3891 

3868 

• 

1  psi  =  6.89  KPa 

average 

=•  2531  psi 

% 


% 


% 
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Concrete  Cylinder 
Cast  2/6/86 
Test  3/10/86 


Grout  Cylinder 
Cast  2/11/86 
Test  3/11/86 


CHAPTER  6 


DESCRIPTION  OF  TEST  SERIES 

A  test  program  was  undertaken  to  investigate  the  behavior  of 
joints  during  quasi-static  and  cyclic  shear  loading  by  varying  factors 
such  as  initial  normal  stress,  roughness  (asperity  angles),  rate  of 
loading,  and  amplitude  of  displacements.  This  chapter  describes  the  pur¬ 
pose  of  the  investigation,  and  details  of  the  tests. 

Purpose  of  the  Test  Series 

The  purpose  of  the  tests  was  to  investigate  the  effects  of  dif¬ 
ferent  factors  on  the  shear  behavior  of  rock  joints  and  to  use  this  data 
base  for  defining  the  parameters  in  the  plasticity  model  described  in 
Chapter  7.  Also,  test  results  not  used  to  find  parameters  for  the  plas¬ 
ticity  model  were  provided,  and  back  predictions  of  these  test  results 
using  the  plasticity  model  serve  as  a  necessary,  though  not  sufficient, 
condition  for  verification  of  the  plasticity  model. 

Factors  studied  include  the  rate  of  loading,  roughness  of  the 
joint,  relative  displacement  during  shear,  the  number  of  cycles  during 
cyclic  loading,  and  state  of  stress. 

Details  of  the  Test  Series 

The  test  series  investigated  both  normal  and  shear  behavior  of 
the  joints.  Normal  behavior  of  the  joints  was  investigated  by  performing 
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normal  load  tests,  described  in  the  next  section.  Shear  behavior  was 
studied  by  performing  large-scale  direct  shear  tests.  Fast  cyclic  and 
slow  cyclic  or  quasi-static  tests  were  performed  in  order  to  ascertain 
the  effects  of  rate  of  loading  as  well  as  number  of  cycles  on  the  shear 
behavior.  Most  of  the  direct  shear  tests  were  displacement-controlled. 
Some  force-controlled  testing  was  performed  for  validation  purposes. 

Normal  Stress  Tests 

The  objective  of  the  normal  stress  tests  was  to  study  the  re¬ 
sponse  of  the  joints  to  normal  stress.  This  entailed  applying  increments 
or  decrements  of  normal  stress  and  measuring  the  vertical  compression  or 
expansion  of  the  joint. 

Method .  Normal  load  tests  were  performed  by  assembling  the  de¬ 
vice  so  that  the  joint  was  mated.  Initial  readings  of  the  vertical 
displacement  measuring  devices  were  taken.  At  this  point,  a  small  ini¬ 
tial  normal  load  was  present  due  to  the  weight  of  the  upper  sample,  and 
the  vertical  load  plate.  Increments  of  normal  load  were  applied  to  the 
upper  sample  in  the  absence  of  shear  forces  in  the  plane  of  the  joint. 
Average  normal  stresses  were  calculated  by  dividing  the  normal  load  by 
the  cross-sectional  area.  As  each  increment  of  normal  stress  was 
applied,  vertical  displacement  readings  were  taken. 

Load  Increments.  In  order  to  determine  the  elastic  normal  stiff¬ 
ness  and  plastic  deformation  of  the  joint ,  cycles  of  loading,  unloading 
and  reloading  were  applied.  A  plot  of  normal  stress  versus  vertical 
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displacement  for  a  typical  normal  stress  test  (9-degree  surface)  is 
presented  in  Fig.  6-1.  Starting  from  the  initial  state  of  stress, 
normal  stress  increments  of  5  psi  (34  KPa)  were  applied  until  a  normal 
stress  of  20  psi  (138  KPa)  was  reached  (0-1).  Then  the  sample  was  un¬ 
loaded  in  decrements  of  5  psi  (138  KPa)  until  the  initial  state  of  stress 
was  reached  (1-2) .  Another  cycle  of  loading  was  applied  in  increments  of 
10  psi  (68.9  KPa)  until  a  stress  level  of  40  psi  (276  KPa)  was  attained 
(2-3-4) ,  and  once  again  the  sample  was  unloaded  to  the  initial  state  of 
stress,  but  this  time  in  decrements  of  10  psi  (68.9  KPa)  (4-5).  Finally, 
the  sample  was  reloaded  in  increments  of  10  psi  (68.9  KPa)  until  a  normal 
stress  of  50  psi  (345  KPa)  was  reached  (5-6). 

Factors.  Normal  stress  tests  were  performed  on  samples  with  dif¬ 
ferent  asperity  side  angles,  including  5,  7,  9  and  15  degrees.  In  this 
manner,  the  effect  of  joint  roughness  on  normal  behavior  was  investigated. 

Quasi-Static  Shear  Tests 

The  purpose  of  the  quasi-static  shear  tests  was  to  subject  the 
joint  to  one  slow  cycle  of  displacement  in  direct  shear.  The  tests  were 
displacement-controlled  with  a  rate  of  displacement  of  0.05  in. /min. 

(1 . 27  mm/min.) . 

Cycle  of  Displacement.  The  tests  were  displacement-controlled 


via  the  control  unit  which  received  a  triangular  displacement  function 
from  the  function  generator.  The  phase  and  amplitude  of  the  function 
were  adjusted  to  obtain  the  waveform  depicted  in  Fig.  6-2(a).  The 
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frequency  was  set  so  that  the  rate  of  displacement  was  .05  in. /min. 

(1.27  mm. /min.)  given  the  maximum  amplitude  of  displacement,  0.25  in. 

(6. 35  mm. ) . 

Figure  6-2 (b)  shows  the  motion  of  the  upper  sample  due  to  the 
input  function  of  Fig.  6-2(a).  Initially,  the  upper  sample  was  in  a 
position  such  that  the  joint  was  fully  mated,  shown  as  datum  in 
Fig.  6-2(b).  The  upper  sample  was  then  displaced  in  the  horizontal  di¬ 
rection  (forward  pass),  about  0.25  in.  (6.35  mm.),  to  a  point  where  the 
peaks  of  the  asperities  from  the  top  and  bottom  surfaces  were  almost  in 
contact.  During  the  forward  pass,  the  upper  sample  traveled  up  the  as¬ 
perity  sides  of  the  bottom  sample.  The  upper  sample  was  then  returned  to 
its  initial  position  (reverse  pass).  During  the  reverse  pass,  the  upper 
sample  traveled  down  the  asperities  of  the  bottom  sample.  The  sample 
was  initially  loaded  in  shear,  unloaded  after  the  upper  sample  was  dis¬ 
placed  to  its  forward  position,  and  then  load  reversal  in  direct  shear 
occurred  as  the  sample  was  displaced  back  to  the  initial  position. 

Method.  To  perform  quasi-static  shear  tests,  the  device  was 
assembled  and  initial  readings  of  the  horizontal  and  vertical  displace¬ 
ment  measuring  devices  were  taken.  A  normal  stress  was  applied,  and 
vertical  displacement  due  to  the  application  of  the  normal  stress  was 
recorded.  A  cycle  of  displacement  was  applied  via  the  input  function 
received  from  the  function  generator.  Readings  of  horizontal  load  and 
the  vertical  and  horizontal  displacement  were  taken  during  the  test. 
Average  shear  stress  was  computed  by  dividing  the  horizontal  load  by  the 


(i)  Datum 


H  b- 


0.25  (in.) 


(ii)  Forward  Pass 


Unloading  in  Shear  Shear  Load  Reversal 

(iii)  Reverse  Pass 


Figure  6-2 (b) .  Motion  of  Upper  Sample  Due  to 
Input  Function  of  Fig.  6-2 (a) 
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projected  cross-sectional  area  of  the  joint.  Figures  6-3(a)  and  (b) 
are  typical  plots  from  a  quasi-static  direct  shear  test  showing  the 
variation  of  shear  stress  (t) ,  and  normal  displacement  (v  ) ,  respec¬ 
tively,  with  horizontal  displacement  (ur) . 

Factors .  Quasi-static  shear  tests  were  performed  on  surfaces 
with  different  roughnesses,  including  cast-smooth,  ground- smooth,  and 
asperity  side  angles  of  5,  7,  9  and  15  degrees.  For  each  roughness, 
tests  were  performed  using  normal  stresses  of  5,  10,  20,  and  50  psi 
(34.5,  69,  138,  and  345  KPa) .  A  test  with  a  normal  stress  of  100  psi 
(690  KPa)  was  also  performed  on  the  cast-smooth  joint. 

Damage.  During  shear  the  surfaces  of  rock  joints  may  become 
damaged  or  worn.  This  damage  may  be  visible  on  the  joint  surfaces  in 
the  form  of  striations ,  broken  asperities,  and/or  gouge.  The  effects  of 
shear  on  the  joint  surface  were  studied  by  comparing  measurements  made 
on  the  joint  surface  before  shearing  to  measurements  made  on  the  joint 
surface  after  shear.  These  measurements  include  profilometer  measure¬ 
ments,  similar  to  those  described  by  Farrington  (1983),  photographs  of 
the  joint  surfaces,  and  measurements  of  the  weight  of  the  gouge  material 
after  shear. 

Multiple  Tests  With  One  Sample.  Between  tests  with  a  different 
normal  stress,  the  device  was  disassembled,  and  the  samples  were  cleaned 
of  debris.  In  this  manner,  the  same  sample  was  used  for  all  tests  having 
a  different  normal  stress.  This  assumes  that  previous  testing  did  not 
cause  significant  damage  to  the  joint  surface. 


This  assumption  was  verified  by  performing  a  test  with  a 
normal  stress  of  20  psi  (138  KPa)  on  a  fresh,  cast-smooth  sample,  and 
comparing  these  results  with  results  of  a  similar  test  performed  on  a 
sample  which  had  been  subjected  to  500  cycles  of  displacement  with  an 
amplitude  of  0.8  in.  (2.03  cm.)  under  a  normal  stress  of  50  psi 
(345  KPa).  Figure  6-4  shows  no  significant  difference  in  the  results 
of  the  two  tests. 

As  a  further  check.,  a  quasi-static  shear  test  was  repeated  on 
the  same  sample  to  see  if  the  first  cycle  of  displacement  affected 
the  results  of  the  second  cycle.  The  device  was  disassembled  between 
the  tests  and  any  gouge  present  removed.  This  investigation  was  done 
on  a  sample  with  an  asperity  side  angle  of  5  degrees  under  a  normal 
stress  of  50  psi  (345  KPa).  The  comparison  of  results  can  be  seen  in 
Fig.  6-5.  These  results  indicate  that  if  the  surfaces  are  cleaned  of 
debris  between  the  tests,  the  results  of  the  second  test  are  not  sig¬ 
nificantly  different  from  the  results  of  the  first  test. 

Fast  Cyclic  Tests 

The  purpose  of  the  fast  cyclic  tests  was  to  study  the  effects 
of  rate  of  displacement  and  number  of  cycles  on  the  shear  behavior  of 
joints.  Joints  were  subjected  to  up  to  500  cycles  of  displacement  at  a 
frequency  of  1  Hz.  Results  from  fast  cyclic  tests  can  be  compared  to 
results  from  quasi-static  tests  to  investigate  the  effect  of  rate  of 
displacement.  Measurements  taken  from  different  cycles  of  displacement 
during  the  test  can  be  compared  to  study  the  effect  of  the  number  of 
cycles  on  the  shear  behavior. 
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r 

Figure  6-4.  Comparison  of  Results  from  Quasi-static  Shear  Testing 
On  a  Fresh  Sample  and  a  Sample  Previously  Subjected 
to  500  Cycles  of  Displacement 
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Figure  6-5.  Comparison  of  Results  from  Two  Quasi-Static  Shear 
Tests  On  the  Same  Sample  With  Debris  Cleaned 
From  the  Shear  Surface  Between  Tests 
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Displacement  Cycle.  The  tests  were  displacement-controlled 
via  the  control  unit  which  received  a  sinusoidal  displacement  function 
from  the  function  generator.  The  displacement  function  had  a  frequency 
of  1  Hz.  The  amplitude  and  phase  of  the  function  were  adjusted  to  ob¬ 
tain  the  symmetric  waveform  depicted  in  Fig.  6-6  (a). 

Figure  6-6 (b)  shows  the  motion  of  the  upper  sample  due  to  the 
input  function  of  Fig.  6-6 (a).  Initially,  the  upper  sample  was  in  a 
position  such  that  the  joint  was  fully  mated  (datum).  The  displacement 
during  one  cycle  can  be  considered  to  be  composed  of  four  phases: 

Phase  1  -  The  upper  sample  was  displaced  in  the  positive  hori¬ 
zontal  direction  until  the  maximum  positive  amplitude  of  cyclic 
displacement,  U^,  was  reached.  During  this  part  of  the  displacement 
cycle,  the  top  sample  rode  up  the  right  asperities  of  the  bottom  sample. 

Phase  2  -  The  displacement  was  reversed,  corresponding  to  un¬ 
loading  and  load  reversal.  During  this  part  of  the  displacement  cycle, 
the  upper  sample  rode  back  down  the  asperities  of  the  bottom  sample 
until  returning  to  the  datum. 

Phase  3  -  Displacement  continued  in  the  negative  direction  as  in 
Phase  2.  After  passing  through  the  datum,  the  upper  sample  traveled  up 
the  left  asperities  in  the  negative  direction  until  reaching  the  nega¬ 
tive  extreme  of  the  maximum  amplitude  of  cyclic  displacement. 

Phase  4  -  Once  again  the  displacement  was  reversed  corresponding 
to  unloading  and  load  reversal.  The  upper  sample  traveled  down  the 
asperities  of  the  bottom  sample  in  the  positive  direction  until 
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Figure  6-6 (a). 
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Figure  6-6 (b).  Motion  of  Upper  Sample  Due  to  the 
Input  Function  of  Fig.  6-6 (a) 
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returning  to  the  datum.  Overtopping  of  asperities  did  not  take  place 
at  any  time  during  the  displacement  cycle. 

Method .  After  adjustment  of  the  displacement  function,  the  pro¬ 
cedure  for  running  a  fast  cyclic  test  was  similar  to  the  procedure  for 
running  a  quasi-static  test,  described  in  the  previous  section.  How¬ 
ever,  for  the  fast  cyclic  tests,  500  cycles  of  displacement  were 
applied.  Figures  6-7 (a)  and  (b)  are  typical  plots  from  a  fast  cyclic 
test  showing  the  variation  of  shear  stress  and  normal  displacement,  re¬ 
spectively,  with  horizontal  displacement  for  different  cycles. 

Factors.  Fast  cyclic  tests  were  performed  on  surfaces  with  dif¬ 
ferent  roughnesses  including  cast-smooth,  ground-smooth,  and  asperity 
side  angles  of  5,  7,  and  9  degrees.  Tests  were  performed  for  normal 
stresses  of  5,  10,  20  and  50  psi  (34.5,  69,  138,  and  345  KPa)  and  maxi¬ 
mum  amplitudes  of  cyclic  displacement  of  0.05,  0.10,  and  0.15  in.  (1.27, 
2.54,  and  3.81  mm.). 

Damage.  Damage  which  may  occur  to  the  joint  surface  due  to  the 
cyclic  shear  displacement  was  studied  in  a  manner  similar  to  that  de¬ 
scribed  for  the  slow  cyclic  tests.  Profilometer  measurements  and 
pictures  of  the  joint  surface  were  taken  before  and  after  shear.  Gouge 
material  present  on  the  joint  surface  after  testing  was  weighed.  These 
measurements  were  compared  for  different  tests. 

Effects  on  Cyclic  Shear  Behavior.  The  effects  of  rate  of  loading 


or  gouge  generated  on  the  joint  surface  due  to  cyclic  shear  displacement 


on  the  cyclic  shear  behavior  of  the  rock  joint  were  also  investigated. 
After  performing  a  fast  cyclic  test,  a  quasi-static  test  under  the  same 
normal  stress  was  performed  without  cleaning  debris  on  the  joint  surface 
generated  during  the  fast  cyclic  test.  Comparisons  were  made  between 
results  from  the  following  "cases:" 

"Case  1"  -  A  quasi-static  test  performed  prior  to  fast  cyclic 
testing  on  a  joint  with  the  same  roughness  and  under  the  same  normal 
stress  (debris  not  present) 

"Case  2"  -  The  first  cycle  of  the  fast  cyclic  test  (debris  not 

present) 

"Case  3"  -  The  quasi-static  test  performed  after  the  fast  cyclic 
test  (debris  present) 

"Case  4"  -  The  last  cycle  of  the  fast  cyclic  test  (debris 
present) . 

Figures  6-8 (a)  and  (b)  present  typical  results  from  testing  with 
a  9-degree  surface  with  normal  stress  20  psi  (138  KPa)  and  Urm  =  0.10  in. 
(2.54  mm.).  Figure  6-8 (a)  shows  that  there  is  no  significant  difference 
between  results  in  terms  of  strength  from  Case  1  and  Case  2. 

Figure  6-8 (b)  indicated  the  same  for  Cases  3  and  4.  Since  the  rate  of 
loading  in  Cases  1  and  2  or  Cases  3  and  4  was  different,  it  may  be  con¬ 
cluded  that  the  rate  of  loading  did  not  have  a  significant  effect  on  the 
shear  behavior  within  the  range  of  rates  investigated.  The  accumulation 
of  gouge  may  affect  the  peak  shear  stress  by  as  much  as  40  percent,  as 
evidenced  by  the  difference  in  results  between  the  tests  with  no  debris 
present  (Cases  1  and  2)  and  the  tests  with  debris  present  (Cases  3  and  4) . 


Ur  (in.) 

Figure  6-8(b) .  Comparison  of  Results  From  a  Quasi-static 
Test  Performed  After  a  Cyclic  Test 
(Surfaces  Not  Cleaned,  Case  3)  and  the 
Last  Cycle  of  a  Cyclic  Test  (Case  4) 
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Multiple  Tests  With  One  Sample.  On  the  basis  of  the  above 
comparison,  it  was  assumed  that  changes  in  the  shear  behavior  of  the 
rock  joints  tested  in  this  program  were  largely  due  to  the  formation  of 
gouge  material.  Therefore,  a  series  of  cyclic  tests  at  various  normal 
stresses  can  be  performed  on  a  single  sample  as  long  as  the  debris  is 
cleaned  from  the  joint  surfaces  between  tests  because  no  significant 
damage  to  the  structure  of  the  surfaces  occurs  as  a  result  of  gouge 
formation. 

Summary  of  Test  Program 

Table  6-1  summarizes  the  tests  performed  during  this  investiga¬ 
tion  including  quasi-static  and  cyclic  shear  tests.  Each  test  has  a 
designation  which  identifies  the  test  series,  the  roughness  of  the 
joint,  type  of  test  and  amplitude  of  cyclic  displacement  if  applicable, 
and  the  normal  stress.  Descriptions  of  these  designations  can  be  found 
at  the  bottom  of  Table  6-1.  In  addition  to  shear  testing,  normal  stress 
tests  were  performed  on  each  surface.  Results  from  this  testing  program 
will  be  discussed  in  Part  II  of  the  final  report  for  this  project  en¬ 
titled  "Constitutive  Modeling  of  Idealized  Rock  Joints  Under  Cyclic 
Loading;  Final  Report,  Part  II."  In  this  report,  test  results  will  only 
be  discussed  with  respect  to  determination  of  parameters  for  the  plas¬ 
ticity  model  and  back  prediction  of  results  with  the  model. 
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Table  4, 


Continued 


Test 


Roughness 


ur  (in) 


an  (psi) 


STATIC 


I-5-S-05  7°  Cast 

0.25 

5 

I-5-S-10 

0.25 

10 

• 

I-5-S-20 

0.25 

20 

I-5-S-50 

0.25 

50 

CYCLIC 

I-5-C2-IQ 

0.10 

10 

% 

I-5-C2-20 

0.10 

20 

STATIC 

I-6-S-05  9°  Cast 

0.25 

5 

I-6-S-10 

0.25 

10 

• 

I-6-S-20 

0.25 

20 

I-6-S-50 

0.25 

50 

CYCLIC 

I-6-C1-05 

0.05 

5 

• 

I-6-C1-10 

0.05 

10 

Static  Test 

After  I-6-C1-10 

_ 

I-6-C1-20 

.05 

20 

% 

I-6-C2-05 

0.10 

5 

I-6-C2-10 

0.10 

10 

I-6-C2-20 

0.10 

20 

Static  Test 

After  I-6-C1-20 

• 

STATIC 

I-2-S-05  15°  Cast 

0.5 

5 

I-2-S-20 

0.5 

0 

I-2-S-50 

0.5 

50 

• 

Table  4,  Continued 


Test  Series  Designation 


I-l-S-10  Series  I  -  Roughness  1  -  Static  -  -  10  psi 

i_2-C1-20  Series  I  -  Roughness  2  Cyclic  TesE  Type  1  - 
*  20  psi,  and  so  on. 

Roughness 


Smooth,  Cast  -  Cast  on  smooth  glass  -  i 
Ground,  Flat  -  Ground  smooth  using  grinding  machine  -  0 
5°  Cast  -  Asperity  angle  =»  5°  -  4 

7°  Cast  -  Asperity  angle  =  7°  -  5 

9°  Cast  -  Asperity  angle  =  9°  -  6 

15°  Cast  -  Asperity  angle  -  15°  -  2 


Table  6-1.  Summary  of  Test  Program 
Explanations  given  below. 
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Test 

Roughness 

ur  (in) 

an  (psi) 

STATIC 

I-l-S-10 

Smooth  Cast 

0.5 

10 

I-l-S-20 

0.5 

20 

I-l-S-50 

0.5 

50 

I-l-S-100 

0.5 

CYCLIC 

100 

I-l-Cl-10 

0.05 

10 

I-l-Cl-20 

0.05 

20 

I-1-C2-10 

0.10 

10 

I-1-C2-20 

0.10 

20 

I-1-C2-50 

0.10 

STATIC 

50 

I-0-S-05 

Ground  Flat 

0.25 

5 

I-0-S-20 

0.25 

20 

I-0-S-50 

0.25 

CYCLIC 

50 

I-0-C2-05 

0.10 

5 

I-0-C2-10 

0.10 

10 

I-0-C2-20 

0.10 

STATIC 

20 

I-4-S-05 

5°  Cast 

0.15 

5 

I-4-S-10 

0.25 

10 

I-4-S-20 

0.25 

20 

I-4-S-50 

0.15 

CYCLIC 

50 

I-4-C2-10 


0.10 


10 


CHAPTER  7 


DESCRIPTION  OF  PLASTICITY  MODEL 

In  order  to  describe  the  behavior  of  jointed  rock  masses 
properly,  it  is  necessary  to  describe  the  behavior  of  the  individual 
rock  joints. 

Rock  joints  in  nature  may  contain  filler  materials,  and  may  be 
saturated  with  water.  The  surface  features  of  a  rock  joint  can  be  highly 
irregular,  and  the  asperities  of  the  top  and  bottom  portions  of  the  joint 
may  be  mated,  or  there  may  be  some  initial  offset  in  the  mating  surfaces. 
The  state  of  stress  to  which  the  rock  joint  is  subjected  may  be  such  that 
shear  takes  place  in  any  arbitrary  direction  with  respect  to  the  asperi¬ 
ties,  and  the  response  of  the  joint  may  be  affected  by  initial  stress  and 
state  of  stress.  Cyclic  or  quasi-static  loading  may  be  applied.  During 
shear  the  joint  may  tend  to  harden  or  soften,  to  compress  or  dilate,  and 
rotation  may  occur  causing  contact  pressures  to  vary.  Development  of  a 
constitutive  model  which  considers  all  of  these  features  is  a  difficult 
task. 

Constitutive  models  based  on  the  theory  of  incremental  plasticity 
have  been  proposed  by  Ghaboussi  and  Wilson  (1973),  Roberds  and  Einstein 
(1978),  Hsu-Jung  (1979),  Plesha  (1986),  Drumm  and  Kane  (1987).  By  em¬ 
ploying  the  theory  of  plasticity,  normal  and  shearing  behavior  of  the 
joint  during  shear  are  coupled.  In  this  study,  use  of  a  generalized 
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yield  function  which  can  model  hardening  behavior  through  a  series  of 
continuous  yield  surfaces  which  expand  towards  an  ultimate  surface 
during  deformation  is  suggested.  To  simplify  the  formulation  in  this 
initial  study,  the  effects  of  filler  material,  pore  water  pressure, 
irregular  surface  features,  initial  offset,  and  softening  are  not  con¬ 
sidered.  Some  of  these  aspects  will  be  considered  in  subsequent  studies. 

The  model  allows  for  effects  of  initial  normal  stress,  states  of 
shear  and  normal  stress,  plastic  hardening,  nonassociativeness ,  volume 
changes  at  joints,  load  reversal  and  cyclic  loading. 

This  chapter  describes  the  specialization  of  a  three-dimensional 
generalized  plasticity  model  to  describe  shear  behavior  of  individual 
rock  joints.  The  specialization  is  described  in  terms  of  monotonic  shear 
loading,  and  then  adaptations  are  described  such  that  the  model  may  be 
used  for  unloading,  load  reversal,  and  cyclic  shear  loading.  The 
physical  significance  of  parameters  as  well  as  salient  details  of  parame¬ 
ter  determination  will  be  discussed. 

Specialization  Of  a  3-D  Generalized  Model 
The  following  generalized  yield  function  written  in  terms  of 
stress  invariants  was  proposed  by  Desai  (1980b) ,  Desai  and  Faruque  (1984), 
Desai  et  al.  (1986),  Desai  et  al.  (1987): 

F  -  J2D  -  +  «12)  (1-S  Sr>“  *  J2D  '  Fb  Fs  ’  0  W'1* 
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where  is  the  first  invariant  of  the  stress  tensor,  J ^  is  the  second 

invariant  of  the  deviatoric  stress  tensor,  S_^.  .  The  parameters  a,  n, 

y,  8  and  m  are  the  response  functions,  and  Sr  equals  stress  ratio  such 
1/3  1/2 

as  /J 2Q  and  Lode  angle,  where  is  the  third  invariant  of 

and  Fg  are  called  basic  and  shape  functions,  respectively. 

Presently  y,  8  and  m  are  treated  as  constants  associated  with 
the  ultimate  envelope,  which  is  defined  as  the  locus  of  points  corre¬ 
sponding  to  asymptotic  stress  to  stress-strain  curves  for  different 
tests.  The  parameter  a  is  chosen  as  a  growth  or  hardening  function  of 
specific  internal  variables  such  as  plastic  strain  trajectory 
E,  =*  /  where  equals  the  plastic  strain  tensor. 

The  above  function  describes  a  series  of  continuous  yield 
surfaces  which  expand  towards  an  ultimate  surface  during  deformations. 

It  has  been  shown  by  Desai  et  al.  (1986)  that  the  function  is  convex  when 
plotted  in  any  stress  space,  gives  orthogonal  intersection  with  the 
axis  when  plotted  in  versus  /j^  space,  and  is  equal  to  zero  on  the 
ultimate  surface  which  intersects  the  axis  at  infinity.  Details  of 
the  hierarchical  approach  are  given  by  Desai  et  al.  (1986)  and  Desai 
et  al.  (1987). 


Analogies 

In  this  investigation  the  shear  behavior  of  idealized  rock  joints 
is  studied.  The  form  of  Eq.  (7-1)  is  specialized  to  describe  yielding  of 
the  joint  plane.  Since  yield  is  to  occur  on  a  specified  plane,  the  func¬ 
tion  is  written  in  terms  of  the  normal,  a  ,  and  shear,  x,  components  of 

n 

stress  across  the  joint  plane,  instead  of  the  stress  invariants.  For 
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solids,  F  controls  the  shape  of  the  yield  function  plotted  in  principal 
s 

stress  space.  When  specialized  to  the  joint  plane,  the  function  is  only 

plotted  in  a  versus  t  space.  Therefore,  F  will  be  taken  as  unity, 
n  s 

The  parameter  n  will  control  the  shape  of  the  yield  function  plotted  in 
an  versus  x  space  and  act  as  the  phase  change  parameter  denoting  transi¬ 
tion  from  contractive  to  dilative  response. 

In  order  to  accomplish  the  specialization,  several  analogies  will 
be  drawn  between  quantities  that  are  defined  for  a  three-dimensional 
state  of  stress  and  for  the  state  of  stress  across  a  joint  plane. 

Table  7-1  summarizes  the  analogies  which  involve  quantities  appearing  in 
the  yield  and  growth  functions.  The  reasoning  behind  the  analogies  is 
discussed  in  the  following  sections. 

Analogies  for  Yield  Function.  Shear  stress,  x,  on  a  joint  plane 
is  considered  to  be  analogous  to  ,  and  normal  stress,  a^,  on  a  joint 

plane  is  considered  to  be  analogous  to  J^.  The  reason  for  this  analogy 
is  the  similar  manner  by  which  these  quantities  affect  the  failure  con¬ 
dition  for  different  states  of  stress.  When  a  three-dimensional  state  of 
stress  is  considered,  J  may  cause  distortion  of  materials  and  may 
add  confinement  which  increases  the  resistance  of  a  material  to  devia- 
toric  stress.  Similarly,  the  shear  stress,  x,  across  a  joint  plane  may 
induce  relative  tangential  displacement  of  the  joint  and  application  of 
normal  stress,  a^,  may  increase  resistance  to  shear  stress.  Using  this 
analogy  between  and  x,  and  and  a^,  Eq.  (7-1)  may  be  written  as 
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Table  7-1.  Analogous  Quantities 


Solid 


Joint 


2D 


£..  =  e 
ii  v 


kk  . 

eij  =  £ij  '  3  6ij 

5  -  /  (dejj 

«d  -  !  (d.J 


/  [(dvrP)2  +  (durP)2]1/2 
/  [(dvrp)2]1/2 
/  [  (du  p)2]1/2 


2D 

£ij 


'ii 


ij 

P 

ij 

P 

kk 

P 


n 


v 


first  invariant  of  stress  tensor 
second  invariant  of  deviatoric  stress  tensor 
strain  tensor 
volumetric  strain 

deviatoric  component  of  strain  tensor 
plastic  strain  tensor 
volumetric  plastic  strain 

deviatoric  component  of  plastic  strain  tensor 
normal  stress 
shear  stress 

relative  vertical  displacement 
relative  tangential  displacement 
plastic  relative  vertical  displacement 
plastic  relative  tangential  displacement 
plastic  strain  trajectory 
volumetric  component  of  £ 
deviatoric  component  of  £ 


* 
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F 


2  ,  n  2 

t  +  aa  -  y  a 
n  n 


0 


(7-2) 


Since  Eq.  (7-2)  is  of  the  same  form  as  Eq.  (7-1),  all  of  the 
properties  of  Eq.  (7-1)  are  retained  including  a  *  0  on  the  ultimate 
surface,  convexity,  and  orthogonality  of  the  intersection  with  the 
axis.  Figure  7-1  presents  a  plot  of  a  typical  function  in  T  versus  cr 
space  for  a  joint  to  be  described  later. 

Analogies  for  Growth  Function.  The  growth  function  needs  to  be 

expressed  in  terms  of  internal  state  variables  which  are  believed  to  be 

related  to  the  hardening  process.  In  order  to  accomplish  this,  two  more 

analogies  need  to  be  drawn  between  quantities  which  describe  deformation 

for  rock  joints  and  those  for  general  three-dimensional  states  of  stress. 

It  is  proposed  that  relative  vertical  displacement,  of  rock  joints  is 

analogous  to  volumetric  strain,  e^,  in  general  3-D  problems,  and  relative 

tangential  displacement,  u^,  of  rock  joints  is  analogous  to  deviatoric 

strain,  e,,  =  -  1/3  <5...  These  analogies  are  based  upon  similar 

ij  ij  kk  ij 

stress-deformation  relationships  when  the  stress  analogies  discussed  in 
the  preceding  section  are  considered. 

Typical  results  from  a  hydrostatic  compression  test  performed  on 
a  granular  material  and  a  normal  stress  test  on  a  rock  joint  are  compared 
in  Fig.  7-2.  The  shape  of  the  versus  z^  is  similar  to  that  of  the 
versus  v  curve. 


r 

Next  consider  the  analogous  stress  quantities  and  T.  An  ini¬ 
tially  isotropic  material  requires  the  application  of  deviatoric  stress 
for  deviatoric  strain  to  occur.  In  rock  joints  shear  stress  induces  rela¬ 
tive  horizontal  displacement. 


Comparison  of  Results  From  a  Hydrostatic  Compression  Test 
and  a  Normal  Stress  Test 


4. 


i 


# 


% 
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Mathematical  Specialization 

It  is  possible  to  derive  Eq.  (7-2)  without  the  analogy  of  the 

preceding  sections.  Consider  a  right  handed  coordinate  system,  as  shown 

in  Fig.  7-3,  where  z  is  the  direction  of  normal  load  application  on  the 

joint  plane,  x  is  the  direction  of  shear,  and  the  y  axis  is  in  the  joint 

plane  perpendicular  to  the  direction  of  shear.  For  the  case  of  direct 

shear,  t  =  t  *  0  since  it  is  assumed  that  there  is  no  relative  dis- 
xy  yz 

placement  of  the  joint  in  the  y  direction,  and  shear  strain  across  the 
joint  plane  is  defined  in  terms  of  relative  displacement.  Also,  it  is 
assumed  that  body  forces  are  negligible.  Thus,  the  following  equilibrium 
equations  for  stress  may  be  written: 


3a 


3x 


.  .  X  ,  xz 

(a)  ~5I  +  — 


3a 

(b)  -  0 


(7-3) 


3x 


xz 


If  the  joint  plane  is  considered  to  be  thin,  is  approxi¬ 

mately  zero  and  Eq.  (7-3)  becomes 


3a 

_ x 

3x 


=  0 


(7-4a) 


3a 


3y 


(7-4b) 


Assuming  for  the  case  of  direct  shear  that  a^  and  a^  are  zero  on 


the  boundaries ,  and  T 


xy 


x  =  a  »  a 

zy  x  y 


0,  and  a  =  a  and  t  =  t  ,  then 
z  n  xz 
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J,  =a  and 
1  n 


(7-5a) 


t  n  2 

J2D  =  T  +  T 


(7-5b) 


Substitution  of  Eq.  (7-5)  in  Eq.  (7-1)  with  F  »  1  yields 


-2—  +  t2  +  a  a11  -  ya2  31  0 

3  n  n 


(7-6) 


which  reduces  to 


t2  +  a an  -  y'  a2  =  0 

n  n 


(7-7) 


where  y'  =  y  -  — ,  which  is  the  same  as  Eq.  (7—2) 


Nonassociativeness 

It  has  been  shown  that  for  many  geologic  materials,  the  plastic 
strain  increment  is  not  normal  to  the  yield  surface  as  assumed  with  the 
associative  flow  rule.  This  can  lead  to  a  nonassociative  flow  rule  where 
the  plastic  strain  increment  is  assumed  to  be  normal  to  the  plastic  po¬ 
tential  function,  Q.  A  simple  approach  has  been  proposed  in  which  Q  is 
obtained  by  correcting  the  growth  function,  a,  in  the  yield  surface  F 
during  plastic  straining  (Desai  and  Siriwardane,  1980;  Frantziskonis  et 
al.  ,  1986).  In  this  case,  Q  is  written  as 

Q  =  t2  +  a.  an  -  ya  2  (7-8) 

Q  n  n 


where 
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a.  -  a  +  <  (aT  -  a)  (1  -  r  )  (7-9) 

Q  I  v 

<  is  a  material  constant  and  a  is  the  value  of  a  at  the  initiation  of 
nonassociativeness ,  which,  for  a  direct  shear  test,  is  considered  to  be 
just  after  the  normal  stress  has  been  applied  prior  to  the  initiation  of 
shear.  The  parameter  r^  is  equal  to 

Elastic-Plastic  Constitutive  Matrix 

In  an  elastic-plastic  constitutive  law,  incremental  stresses  and 
incremental  strains  are  related  through  an  elastic-plastic  constitutive 
matrix.  The  derivation  involves  decomposition  of  deformations  into 
elastic  and  plastic  parts,  and  a  stress-elastic  strain  relationship. 
These,  together  with  either  an  associative  or  nonassociative  flow  rule, 
which  relates  plastic  strains  to  the  state  of  stress,  the  consistency 
condition,  and  considering  a  to  be  a  function  of  the  history  of  plastic 
strain  trajectories,  yield  the  following  expression  for  the  elastic- 
plastic  constitutive  matrix: 


which  is  symmetric  for  the  associative  case  and  nonsymmetric  for  the  non¬ 
associative  case.  Details  of  this  derivation  can  be  found  in  Appendix  I. 
In  this  formulation,  relative  displacements  are  treated  as  strains.  For 
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input  in  a  finite  element  program,  strain  would  be  computed  from  rela¬ 
tive  displacements  by  considering  the  joint  to  be  of  finite  thickness 
and  dividing  relative  displacements  by  this  thickness  (Desai  et  al. , 
1982). 

Load  Reversal 

Load  reversal  causes  the  shear  direction  of  the  upper  sample 
with  respect  to  the  lower  sample  to  be  reversed.  The  response  of  the 
joint  after  load  reversal  depends  on  two  factors: 

1.  Displacement-induced  anisotropic  response  of  the  joint  with 
respect  to  shear  displacement. 

2.  The  accumulation  of  gouge  and  damage  which  occurred  as  a  re¬ 
sult  of  shear  displacement  prior  to  load  reversal. 

Displacement-Induced  Anisotropy.  The  displacement-induced  ani¬ 
sotropic  response  of  rock  joints  will  be  discussed  with  reference  to  the 
slow  cyclic  tests  described  in  Chapter  6.  Consider  the  cycle  of  dis¬ 
placement  depicted  in  Fig.  6-2 (b).  Initially,  the  upper  sample  is  in  a 
position  such  that  the  rock  joint  is  considered  fully  mated.  At  this  in¬ 
stant,  the  response  of  the  joint  during  shear  is  independent  of  the 
direction  of  shear,  and  the  joint  is  considered  to  be  isotropic  with  re¬ 
spect  to  the  direction  of  shear.  As  the  upper  sample  is  displaced 
during  the  forward  pass,  it  travels  up  the  asperity  sides  of  the  bottom 
sample.  When  the  direction  of  the  upper  sample  is  reversed,  the  vertical 
displacement  that  is  coupled  to  the  shear  displacement  is  also  reversed 
as  the  upper  sample  travels  back  down  the  asperities  of  the  bottom  sample 


returning  to  the  datum. 
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The  reversal  of  travel  up  or  down  an  asperity  side  not  only 
changes  the  deformation  characteristics  of  the  joint,  but  also  the  shear 
strength.  This  can  easily  be  realized  by  referring  to  Fig.  7-4,  which 
depicts  a  plane  inclined  to  the  direction  of  the  driving  force.  The 
failure  criteria  relates  the  shear  strength  of  the  joint  to  the  normal 
stress,  asperity  angle  and  direction  of  travel  (up  or  down  incline). 
Hence,  after  shifting  the  upper  sample  from  the  initial  mated  joint  con¬ 
dition,  the  shear  behavior  of  the  joint  depends  upon  the  direction  of 
shear,  and  the  joint  must  be  considered  to  be  anisotropic  with  respect 
to  shear. 

Displacement-induced  anisotropy  is  modeled  by  adaptation  of  the 

isotropic  hardening  model,  as  shown  in  Fig.  7-5.  After  the  joint  is 

displaced  from  the  mated  position  (0-1) ,  parameters  for  the  yield  and 

growth  functions  are  for  a  forward  or  reverse  direction.  Determination 

of  these  parameters  will  be  discussed  in  the  next  section.  Thus,  a 

different  function  denoted  F  ,  F. ,  F„  or  F_  is  used  to  define  yield,  de- 

o  1  Z  J 

pending  on  the  sign  of  the  shear  stress.  Note  that  there  is  no 
discontinuity  in  the  yield  between  unloading  and  load  reversal  (2-3)  be¬ 
cause  yield  functions  for  both  forward  and  reverse  passes  have 
orthogonal  intersections  with  the  abscissa. 

Prior  to  load  reversal,  the  joint  is  unloaded  in  shear  (1-2),  as 
depicted  in  Fig.  7-5.  Currently,  unloading  is  considered  to  be  elastic. 
However,  after  unloading,  the  yield  surface  collapses.  Thus,  when  the 
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shear  stress  is  reduced  to  zero,  the  joint  is  under  the  action  of  the 
applied  normal  stress,  and  just  prior  to  load  reversal  (2-3),  the  stress 
point  lies  on  a  current  yield  surface. 

Gouge  and  Damage.  Gouge  and  damage  generated  on  the  joint 
surfaces  due  to  shear  displacement  prior  to  load  reversal  tend  to  change 
the  state  of  the  joint,  affecting  the  shear  behavior.  In  order  to  de¬ 
scribe  the  shear  behavior  after  load  reversal  properly,  the  model  must 
consider  the  change  in  state  which  has  occurred  to  the  joint  surface  as 
a  result  of  the  forward  pass.  This  is  considered  since  the  evaluation  of 
parameters  for  the  reverse  pass  are  determined  from  test  data  obtained 
after  the  sample  has  been  subjected  to  shear  from  the  forward  pass. 

Application  of  the  Plasticity  Model  to  Cyclic  Loading 

In  this  section  development  of  the  aforementioned  plasticity 
model  to  the  case  of  cyclic  shear  loading  is  presented.  Each  cycle  of 
displacement  is  considered  to  consist  of  four  phases,  as  discussed  in 
Chapter  6.  These  phases  describe  the  direction  of  travel  of  the  top 
sample  with  respect  to  the  bottom  sample.  The  isotropic  hardening  model 
is  applied  to  each  phase  of  the  displacement  cycle.  Parameters  are  de¬ 
termined  for  each  phase  of  a  given  cycle,  and  then  parameters  for  a  given 
phase  may  be  expressed  as  functions  of  number  of  cycles. 

For  the  sake  of  completeness,  descriptions  of  each  phase  of  the 
displacement  cycle  for  a  cyclic  shear  test  are  given  here.  This  dis¬ 
cussion  is  with  reference  to  Fig.  6-6 (b). 


Phase  1:  Starting  from  an  initially  mated  condition  (datum), 
the  upper  sample  is  displaced  up  the  right  asperities  (positive  direc¬ 
tion)  an  amount  equal  to  the  maximum  amplitude  of  cyclic  displacement. 

Phase  2:  The  upper  sample  is  unloading  in  shear,  and  the  direc¬ 
tion  of  shear  is  reversed  (negative  direction).  During  this  phase, 
elastic  shear  displacement  occurs  first.  Then  elasto-plastic  deforma¬ 
tions  due  to  load  reversal  develop  and  continue  until  the  upper  sample 
returns  to  the  datum. 

Phase  3:  Displacement  continues  in  the  negative  direction  until 
the  maximum  amplitude  of  cyclic  displacement  is  reached.  During  this 
phase,  the  upper  sample  travels  up  the  left  asperity  sides  of  the  bottom 
sample  in  the  negative  direction.  With  the  exception  of  cycle  one,  to  be 
discussed  later,  this  phase  is  similar  to  Phase  1,  only  the  direction  of 
tangential  displacement  is  reversed. 

Phase  4:  Once  again,  the  direction  of  travel  is  reversed  (posi¬ 
tive  direction) ,  during  which  initial  elastic  unloading,  followed  by 
plastic  deformation,  occurs  until  the  upper  sample  returns  to  the  datum. 
This  phase  is  similar  to  Phase  2. 

As  previously  discussed,  after  unloading  the  yield  surface  is 
assumed  to  collapse  so  that  upon  load  reversal  the  stress  point  lies  on  a 
current  yield  surface.  Between  Phases  2  and  3  and  4  and  1,  the  parame¬ 
ters  change  to  reflect  the  transition  from  down  to  up  asperity  side 
travel.  No  unloading  occurs,  but  there  is  a  shift  in  the  yield  surface 
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depicted  in  Fig.  7-6.  This  shift  in  the  yield  surface  occurs  because 
prior  to  the  phase  transition  the  response  is  approaching  ultimate,  but 
the  instant  after  the  phase  transition  the  response  is  no  longer  near 
ultimate. 

It  is  important  to  point  out  that  in  the  present  testing,  at  no 
time  during  the  deformation  history  do  the  asperities  of  the  top  sample 
override  the  asperities  of  the  bottom  sample  as  depicted  in  Fig.  7-7. 

This  would  cause  an  instability  with  regard  to  the  state  of  stress,  which 
is  not  currently  built  into  the  model. 

Cyclic  Hardening.  Changes  in  the  behavior  of  the  joint  plane  as 
a  result  of  cyclic  shear  loading  are  assumed  to  be  described  primarily 
through  a  change  in  the  ultimate  response.  This  is  based  on  observations 
made  from  the  cyclic  test  series.  The  results  of  this  test  series  indi¬ 
cate  that  the  ultimate  stress  during  Phases  1  and  3  (up  asperity  travel) 
increases  with  number  of  cycles.  The  history  of  deformation  which  oc¬ 
curred  during  previous  cycles  is  considered  since  the  parameters 
determined  during  a  given  cycle  are  based  on  deformations  which  occur  to 
a  joint  surface  that  has  experienced  deformation  the  previous  cycles. 

Parameters 

The  elastic-plastic  model  discussed  herein  requires  constants  for 
the  elastic  response  of  the  joint  as  well  as  parameters  which  appear  in 
the  yield  and  potential  functions.  The  elastic  response  of  the  joint  is 
described  through  two  constants,  ks  and  k  ,  the  elastic  shear  and  normal 
stiffnesses,  respectively.  The  yield  and  potential  function  utilize 
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parameters  which  describe  the  ultimate  state,  the  growth  or  hardening 
behavior,  and  nonassociativeness.  The  physical  significance  of  the 
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various  parameters  and  salient  details  of  their  determination  are  dis¬ 
cussed  in  this  section.  Parameter  determination  from  test  data  will  be 
discussed  in  Chapter  8. 

Elastic  Constants 

The  elastic  shear  stiffness,  kg ,  is  the  slope  of  the  unloading 
part  of  the  shear  (x)  versus  relative  tangential  displacement  (u^)  curve 
obtained  from  a  direct  shear  test,  as  shown  in  Fig.  7-8.  This  value 
varies  depending  on  the  normal  stress  across  the  joint  plane  at  the  ini¬ 
tiation  of  shear. 

The  elastic  normal  stiffness,  k  ,  is  the  initial  slope  of  a 

n 

normal  stress  (a  )  versus  relative  vertical  displacement  (v  )  curve  ob- 
n  r 

tained  from  a  normal  stress  test,  as  depicted  in  Fig.  7-9. 

Parameters  for  the  Yield  Function 

The  following  sections  describe  procedures  for  finding  parameters 
for  the  yield  function  including  the  ultimate  parameter  y,  the  phase 
change  parameter  n,  and  the  growth  function. 

Ultimate  Parameters.  As  previously  discussed,  the  growth  func¬ 
tion  is  equal  to  zero  on  the  ultimate  surface.  Thus  Eq.  (7-2) 
specializes  to 
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where  x^  is  usually  ten  to  twenty  percent  greater  than  the  peak  shear 
stress,  as  shown  in  Fig.  7-10.  The  only  parameter  which  needs  to  be  de¬ 
termined  is  y.  It  is  found  by  squaring  the  slope  of  the  "best  fit"  line 
from  a  plot  of  t  versus  a  (Fig.  7-11). 

Parameter  n.  The  parameter  n  may  be  found  by  considering  the 
point  during  shear  where  the  change  in  vertical  displacement  is  zero,  as 
shown  in  Fig.  7-12.  The  shear  stress  at  this  point  is  T  .  If  the 

plastic  strain  increment  is  assumed  to  be  normal  to  the  yield  surface, 

3F 

then  at  this  point  =  0.  From  Eq.  (7-2) 

n 

II-  =  naa  n_1  -  3ya  =0  (7-12) 

3a  n  n 

n 


and 


a 

n 
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=  ( 


3y.n 

nor 


2 


(7-13) 


If  Mzs  is  defined  as  the  slope  of  the  line  yielding  the  locus  of 
points  on  a  plot  of  x  versus  (Fig.  7-13)  ,  then 


x 

=  M  2  -  -aa  n~2  +  y  =  y  (1--)  (7-14) 

2  zs  n  n  / 

a 

n 


The  value  of  n  may  be  determined  from  this  expression  once  M  is  known. 

Growth  Function.  In  general,  the  growth  function  can  be  con¬ 
sidered  a  function  of  the  history  of  plastic  strain,  plastic  work,  and 
internal  state  variables  (Desai  and  Faruque,  1984;  Desai  et  al. ,  1986; 


, 


< 


< 


< 


T 


T 


T 


T 
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and  Desai  et  al. ,  1987).  In  this  study,  relationships  between  a  and  the 
history  of  plastic  strain  trajectories,  and  relationships  between  a  and 
plastic  work  were  investigated.  It  was  found  that  plotting  a  versus  the 
total  trajectory  of  plastic  strains,  E,  =  /  (u^V  +  (v^V^.  yielded 
the  most  consistent  trends.  A  plot  of  a  versus  £  can  be  seen  in 
Fig.  7-14.  The  shape  of  this  plot  indicates  that  a  relationship  of  the 
form 

ct  =  a  £b  (7-15) 

where  a  and  b  are  material  constants  can  be  appropriate.  The  constant  b 
is  negative  such  that  a  becomes  smaller  as  £  increases. 

Nonassociative  Parameter 

The  nonassociative  parameter,  k,  may  be  determined  by  considering 
the  relation  between  the  vertical  and  horizontal  displacement;  as  a 
simplification,  the  slope  near  ultimate  conditions  can  be  used 
(Frantz iskonis  et  al. ,  1986).  This  may  be  shown  as  follows: 

Consider  Eq.  (7-8) 

aQ  =  a  +  <  (c^  -  a)  (1  -  r ^ 

at  ultimate  conditions  a  equals  zero,  hence 

aQ  =  koi  (1  -  rv) 
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so  that 


< 


Vai 


Now,  from  the  flow  rule 


dv  P 

( — ”)  Ult 
du  P 


(IQ-zlQ)  ult 


(7-16) 


(7-17) 


and  we  obtain 


dv  F  /  1  \ 

(7“V  Ult  =  (aQ  n  Qn  n' 

du 

r 


-  2y a  ) / 2t 
n 


ult 


Therefore, 


dv  .  n 

a  -  [(-JL-)  ult  2xult)  +  2yan)/nan^-1> 
^  du 


(7-18) 


Load  Reversal 

As  previously  discussed,  after  load  reversal,  parameters  for  the 
yield  and  growth  function  must  be  determined  for  the  reverse  pass.  The 
procedure  is  similar  to  that  for  the  forward  pass. 

Yield  and  Growth  Function.  The  ultimate  state  parameter  y,  the 
parameter  n,  and  the  constants  a  and  b  in  the  growth  function  must  be 
determined  for  the  reverse  pass.  The  value  of  y  for  the  reverse  pass  is 
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found  in  the  same  manner  as  y  for  the  forward  pass.  The  ultimate  shear 
stress,  t  ,  is  plotted  versus  for  the  reverse  pass,  and  y  is  taken 
as  the  square  of  the  slope  of  the  "best  fit"  line  passing  through  these 
points . 

The  same  value  for  the  parameter  n  which  was  used  during  the 
forward  pass  is  adopted  for  the  reverse  pass.  This  is  because  the  phase 
change  behavior  is  assumed  to  be  related  to  damage  occurring  to  smaller 
asperities  superimposed  on  the  sides  of  the  larger  sawtooth  asperities. 
These  smaller  asperities,  or  microasperities,  are  related  to  the  com¬ 
ponents  of  the  concrete  mix  as  discussed  by  Kulhawy  and  Peterson  (1979). 

Parameters  a  and  b  for  the  growth  function  during  the  reverse 
pass  are  found  after  y  is  determined  for  the  reverse  pass.  The  proced¬ 
ure  is  identical  to  that  for  the  forward  pass.  The  parameters  are 
determined  from  a  plot  of  a  versus  £  for  the  reverse  pass. 

Cyclic  Loading 

Parameters  to  be  determined  for  each  phase  of  displacement 
during  a  given  cycle  include  the  ultimate  parameter,  y,  the  phase  change 
parameter,  n,  and  parameters  for  the  growth  function.  Due  to  the  simi¬ 
larities  between  Phases  1  and  3,  and  Phases  2  and  4  (Fig.  6-6(b)),  only 
two  sets  of  ultimate  parameters  are  needed  to  describe  a  given  cycle  of 
displacement.  The  same  is  true  of  the  parameters  for  the  growth  func¬ 
tion  with  the  exception  of  Phase  1  during  the  first  cycle  of  displacement. 
The  same  value  for  n  is  used  for  all  phases  since  the  phase  change  be¬ 
havior  is  assumed  to  be  similar  for  all  phases  of  the  displacement  cycle. 
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Growth  function  parameters  are  determined  based  on  the  values 
of  the  ultimate  parameters  used  for  a  given  cycle  of  displacement  and 
may  be  functions  of  number  of  cycles  themselves. 

Summary 

This  section  has  described  the  specialization  of  a  3-D  General¬ 
ized  Plasticity  Model  to  describe  the  behavior  of  rock  joints  under  the 
action  of  shear  and  normal  stresses.  The  model  was  extended  beyond  the 
case  of  monotonic  shear  loading  to  include  shear  load  reversal  and 
cyclic  shear  loading.  Methods  for  the  determination  of  parameters  were 
detailed.  In  the  next  chapter,  parameters  will  be  determined  from  a  set 


of  test  data  and  back  predictions  performed  in  order  to  demonstrate 
application  of  the  model. 


CHAPTER  8 


BACK  PREDICTIONS 

Parameters  for  the  plasticity  model  discussed  in  Chapter  7  are 
determined  from  the  test  data.  These  parameters  characterize  the  be¬ 
havior  of  the  joint  plane  during  shear  as  generally  as  possible.  Using 
these  determined  parameters,  back  predictions  of  experimental  results 
were  performed  for  three  cases  of  shear:  quasi-static  loading,  quasi¬ 
static  loading  with  unloading  and  load  reversal,  and  cyclic  shear  loading, 
or  displacement  with  the  frequency  of  1  Hz,  as  described  in  Chapter  6. 

This  chapter  describes  the  determination  of  parameters  for  these  three 
cases  of  shear.  Back  predictions  of  test  results  from  the  three  cases 
are  presented  and  discussed.  Because  of  the  large  number  of  figures  de¬ 
picting  back  predictions,  only  results  from  the  9-degree  surface  will  be 
presented  in  this  chapter.  Other  back  predictions  are  presented  in 
Appendix  II. 


Quasi-Static  Loading 

The  case  of  quasi-static  loading  was  investigated  with  test  re¬ 
sults  for  four  different  surfaces.  These  include  the  ground  flat  surface, 
and  surfaces  with  asperity  side  angles  of  5 ,  7  and  9  degrees. 

Determination  of  Parameters 

Five  parameters  need  to  be  determined  in  order  to  describe  the 
case  of  quasi-static  loading.  These  include  the  ultimate  parameter,  y, 
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the  phase  change  parameter  n,  growth  function  parameters  a  and  b,  and, 
for  the  case  of  the  nonassociative  flow  rule,  the  parameter  K.  The  pro¬ 
cedures  for  determination  of  these  parameters  were  discussed  in 
Chapter  7.  In  this  section,  the  data  used  to  determine  the  parameters 
is  presented,  and  the  values  from  various  surfaces  are  compared  and  dis¬ 
cussed.  Table  8-1  presents  a  summary  of  parameters  obtained  from  the 
different  surfaces. 

Ultimate  Parameter  y.  Figure  8-1  is  a  plot  of  ultimate  shear 
stress  versus  normal  stress  obtained  from  quasi-static  shear  tests  per¬ 
formed  with  the  9-degree  surface.  A  "best  fit"  line  is  passed  through 
the  data  points  based  on  a  linear  regression  analysis.  This  line  de¬ 
fines  the  ultimate  surface,  and  the  square  of  its  slope  yields  the  value 
of  y.  Table  8-1  presents  the  values  of  y  obtained  for  the  different 
surfaces.  The  parameter  y  is  lowest  for  the  flat  surface  and  consis¬ 
tently  increases  with  roughness. 

Phase  Change  Parameter  n.  Figure  8-2  is  a  plot  of  the  shear 
stress  at  the  point  of  phase  change  versus  normal  stress  (Phase  Change 
Envelope)  for  the  9-degree  surface.  It  should  be  noted  that  the  point 
of  phase  change  during  testing  of  rock  joints  is  not  always  a  well 
defined  point.  Many  times,  only  an  initial  flat  portion  of  the  curve  is 
evident  prior  to  dilatation  rather  than  a  measured  compression.  This 
accounts  for  the  evident  scatter  in  the  data. 

Using  a  linear  regression  analysis,  a  "best  fit"  line  was  de¬ 
termined  from  the  data  points  for  each  surface.  Based  on  the  slope  of 
this  line  and  the  values  of  y  already  determined,  a  value  of  n  was 
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Table  8-1.  Summary  of  Parameters  Obtained  From 
Different  Surfaces  for  the  Case  of 
Quasi-Static  Loading 


Parameter 

Surface 

Flat 

5° 

7° 

9° 

Y 

0.36 

0.42 

0.78 

0.31 

n 

2.50 

2.50 

2.50 

2.50 

a 

0.02312 

0.01063 

0.03053 

0.04736 

b 

-0.116 

-0.293 

-0.223 

-0.162 

< 

0.70 

0.70 

0.70 

0.55 

Average  Values  of  Growth  Function  Parameters 
a  and  b,  and  Nonassociative  Parameter  < 

a  -  0.02791 

b  =  -0.1985 

<  =  0.66 


(psi) 


(psi) 


computed  for  each  surface.  No  measurable  dilation  or  compression  oc¬ 
curred  during  the  testing  of  the  ground  flat  surface,  and  this  indicates 
that  a  value  of  n  equal  to  2.0  should  be  used  for  this  case  since  M 

zs 

is  zero  [see  Eq.  (7-14)].  Table  8-2  presents  the  value  of  n  determined 
for  each  surface.  Due  to  the  scatter  in  the  data,  and  the  small  varia¬ 
tion  in  n  between  surfaces,  an  average  value  of  n  =  2.5  was  used  for  all 
surfaces.  By  using  the  same  value  of  n  for  all  surfaces,  it  is  assumed 
that  the  point  of  phase  change  during  shear  of  rock  joints  is  more  re¬ 
lated  to  the  microasperity  structure,  which  may  be  similar  between  these 
surfaces,  and  the  manner  in  which  the  macroasperities  of  the  top  and 
bottom  samples  "fit"  prior  to  the  initiation  of  shear,  than  to  varia¬ 
tions  of  the  sawtooth  (macroasperity)  structure  apparent  between  the 
surfaces  investigated  during  this  study. 

Growth  Function  Parameters.  Figure  8-3  is  a  plot  of  growth 
function  values,  a  versus  the  trajectory  of  the  total  plastic  displace¬ 
ments  E,,  plotted  from  data  for  the  9-degree  surface.  The  figure 
presents  data  from  tests  performed  under  different  levels  of  normal 
stress. 

Assuming  a  form  for  the  growth  function  a  ■  a^,  as  discussed 
in  Chapter  7,  the  parameters  a  and  b  were  determined  via  a  linear  re¬ 
gression  analysis  of  the  data  points  plotted  on  a  transformed 
logarithmic  plot.  The  parameters  a  and  b  are  considered  as  constants 
over  the  range  of  applied  normal  stresses  for  a  given  surface. 

Table  8-1  presents  the  values  of  the  growth  function  constants  a  and  b 
for  different  surfaces. 
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The  values  of  a  and  b  for  the  different  surfaces  presented  in 
Table  8-1  may  be  averaged  to  yield  growth  function  constants  to  be  used 
for  all  the  different  roughnesses.  These  values  are  shown  at  the  bottom 
of  Table  8-1.  By  using  the  same  growth  function  constants  for  all  the 
roughnesses ,  it  is  ass.jned  that  the  hardening  process  is  more  related  to 
the  microasperity  structure  than  to  the  difference  in  macroasperities 
between  the  surfaces.  Back  predictions  using  growth  function  constants 
determined  for  each  roughness  and  for  all  the  surfaces  will  be  compared 
in  the  next  section. 

Nonassociative  Parameter  K.  The  nonassociative  parameter  <  is 
computed  from  test  results  based  on  the  slope  of  the  vertical  relative 
displacement  versus  tangential  relative  displacement  near  ultimate  con¬ 
ditions,  as  outlined  in  Chapter  7.  These  computations  are  presented  in 
Appendix  III.  During  back  predictions,  it  was  found  that  for  best  re¬ 
sults,  these  computed  values  may  need  to  be  adjusted  slightly.  This  may 
be  due  to  the  fact  that  the  procedure  for  the  determination  of  k  assumes 
it  to  be  a  constant  where  in  reality  it  may  not  be.  Table  8-1  presents 
the  values  of  tc  used  in  the  back  prediction  of  test  results  for  the 
different  roughnesses. 

Back  Predictions  of  Quasi-Static  Loading 

Using  the  parameters  discussed  in  the  preceding  section,  back 
predictions  were  made  and  compared  with  experimental  results.  Both 
associative  and  nonassociative  flow  rules  were  employed.  Also,  back 
predictions  were  performed  with  averaged  growth  function  constants  and 
nonassociative  parameter  k. 
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Algorithm.  The  back  predictions  were  performed  by  numerical 
integration  of  the  elastic-plastic  constitutive,  Ke^,  matrix,  described 
in  Chapter  7.  Because  of  the  hierarchical  nature  of  the  plasticity 
model,  the  same  formulation  was  used  for  both  the  associative  and  non- 
associative  cases.  For  the  associative  case,  <  =  0.  The  integration 
is  performed  using  an  incremental  technique  similar  to  one  described  by 
Desai  and  Faruque  (1984).  Small  increments  of  relative  tangential  dis¬ 
placements  are  input  and  the  corresponding  increments  of  shear  stress 
and  relative  vertical  displacement  are  computed  under  the  condition  of 
constant  normal  stress.  As  with  all  numerical  techniques,  the  solution 
obtained  in  chis  manner  is  an  approximate  one,  the  accuracy  of  which  in¬ 
creases  with  a  decrease  in  the  increment  size.  In  this  study,  relative 
tangential  displacement  increments  of  0.0001  in.  (0.0025  mm.)  were 
found  to  yield  acceptable  results.  The  algorithm  employed  for  these 
back  predictions  is  as  follows: 

1.  Input  the  normal  stress,  increments  of  relative  tangential 
displacement,  parameters  for  the  joint  plane  including  y,  n,  a,  b,  <, 
and  the  elastic  shear  and  normal  stiffnesses. 

2.  Starting  from  a  state  of  zero  displacement  and  stress,  apply 
the  normal  stress  and  compute  the  value  of  the  growth  function,  a,  based 
on  the  state  of  stress. 

3.  Set  the  value  of  a  at  the  initiation  of  nonassociative  be¬ 
havior,  a equal  to  the  value  of  a  just  after  the  application  of  normal 
stress,  but  prior  to  the  initiation  of  shear.  At  this  point,  the 
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vertical  component  of  the  trajectory  of  the  history  of  plastic  strain, 

E  »  is  equal  to  the  total  trajectory  of  the  history  of  plastic  strain, 

E  (computed  from  a),  and  the  ratio  of  Ev/E»  ,  is  equal  to  unity. 

4.  Using  values  of  a,  ,  r^  and  the  inputted  parameters,  com¬ 
pute  the  modified  growth  function,  ,  used  in  the  potential  function. 
Compute  the  derivatives  of  the  yield  and  potential  surfaces,  and  form 
the  KeP  matrix. 

5.  Solve  the  system  of  equations  and  determine  the  increment  of 
shear  stress  and  relative  vertical  displacement. 

6.  Update  the  total  stress  and  relative  displacement  vectors. 

7.  Compute  increments  of  plastic  relative  displacement  using 
the  flow  rule,  associative  or  nonassociative. 

8.  Obtain  updated  values  for  E»  Ev,  and  a. 

9.  If  this  is  the  first  increment  of  shear  displacement,  check 
the  value  of  the  yield  function,  F.  If  F  is  negative,  adjust  the  shear 
stress  to  make  F  =  0.  This  is  done  by  changing  the  elastic  shear  stiff¬ 
ness  on  the  first  increment  such  that  F  is  equal  to  zero  at  the  end  of 
the  increment . 

10.  If  the  total  relative  tangential  displacement  is  not 
reached,  return  to  Step  4.  If  the  total  tangential  displacement  is 
reached,  stop. 

Results .  Figures  8-4  to  8-8  present  back  predictions  of  quasi¬ 
static  shear  tests.  Both  test  data  and  back  predicted  results  are 
shown.  Parameters  utilized  in  this  analysis  are  presented  in  Table  8-1, 
where  the  growth  function  constants  and  k  are  found  for  each  surface. 
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Figure  8-6.  Test  Data  and  Back  Predictions  of  v  versus  u  ; 

Quasi-Static  Loading;  0  =10  psi,  §°  Surfaced 

Associative  and  Nonassociative  Case 
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Figure  8-7.  Test  Data  and  Back  Predictions  of  v  versus  u 
Quasi-Static  Loading;  CT  =20  psi;  5°  Surface 
Associative  and  Nonassociative  Case 
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Figure  8-8.  Test  Data  and  Back  Predictions  or  v  versus  ur; 

Quasi-Static  Loading;  o  =50  psi;  9°  Surface; 
Associative  and  Nonassociative  Case 


185 


Figure  8-4  presents  plots  of  shear  stress  versus  relative  tan¬ 
gential  displacement  for  the  9-degree  surface  under  different  levels  of 
normal  stress.  Use  of  the  associative  or  nonassociative  flow  rule  does 
not  significantly  alter  the  back  prediction  of  these  curves;  hence, 
only  back  predictions  utilizing  the  associative  flow  rule  are  presented 
here.  The  model  captures  the  effect  of  different  levels  of  normal 
stress  and  hardening  behavior  quite  well.  Back  predicted  curves  appear 
to  be  in  reasonable  agreement  with  test  data.  The  only  significant 
discrepancy  is  at  a  normal  stress  of  50  psi  (344.5  KPa) . 

Figures  8-5  to  8-8  are  plots  of  relative  tangential  displace¬ 
ment,  ur>  versus  relative  vertical  displacement,  v^.  Each  figure 
depicts  results  from  different  levels  of  normal  stress  including  back 
predictions  using  associative  and  nonassociative  flow  rules  and  the  test 
data.  In  all  cases,  the  back  predictions  which  employ  the  associative 
flow  rule  significantly  overpredict  dilation.  The  back  predictions 
which  utilize  the  nonassociative  flow  rule  are  in  better  agreement  with 
test  data.  However,  some  error  between  experiment  and  theory  still 
exists.  The  model  appears  to  be  much  more  sensitive  to  the  effect  of 
changes  in  levels  of  normal  stress  on  the  prediction  of  relative 
vertical  displacements  during  shear  than  experimental  results  indicate. 
At  normal  stresses  of  5  and  10  psi  (34.5  and  68.9  KPa),  the  back  pre¬ 
dictions  utilizing  the  nonassociative  flow  rule  still  overpredict 
dilation.  At  20  psi  (137.8  KPa),  the  back  prediction  is  good,  and  at 
50  psi  (344.5  KPa),  not  quite  enough  dilation  is  predicted. 
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Average  Parameters.  In  this  case,  growth  function  constants,  a 
and  b,  and  the  parameter  <  are  averaged  over  the  range  of  surface  rough¬ 
nesses  considered.  The  ultimate  parameter  y  is  still  determined 
separately  for  each  roughness.  Therefore,  in  this  discussion,  it  is  im¬ 
plicitly  assumed  that  y  is  significantly  affected  by  the  macroasperity 
characteristics  of  the  surface,  while  hardening  and  nonassociativeness 
are  more  affected  by  characteristics  of  the  microasperities. 

Figures  8-9  to  8-11  present  results  of  back  predictions  uti¬ 
lizing  the  average  parameters  presented  in  Table  8-1.  Only  back 
predictions  utilizing  a  nonassociative  flow  rule  are  considered.  Re¬ 
sults  are  presented  for  normal  stresses  of  5  and  50  psi  (34.5  and 
344.5  KPa) .  Comparing  the  results  utilizing  average  parameters  with  the 
results  of  the  previous  section,  it  appears  that  there  are  some  differ¬ 
ences,  but  the  differences  are  not  striking  for  the  9-degree  surface. 


The  case  of  quasi-static  loading  with  unloading  and  load  re¬ 
versal  was  investigated  with  test  results  from  the  9-degree  surface.  As 
described  in  Chapter  7,  this  case  is  modeled  by  considering  two  distinct 
phases  of  displacement,  the  forward  and  reverse  passes.  During  the  for¬ 
ward  pass,  the  joint  is  dilating  as  "up  asperity  travel"  takes  place 
(Fig.  6-2(b)) .  Parameters  for  this  phase  are  identical  to  those  discussed 
in  the  preceding  section  on  quasi-static  loading.  During  the  reverse 
pass,  the  joint  is  compressing  as  "down  asperity  travel"  takes  place 
(Fig.  6-2(b).  Parameters  for  this  pass  are  found  in  a  manner  similar  to 
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that  for  the  forward  pass,  previously  described.  The  following  section 
describes  the  determination  of  model  parameters  from  the  reverse  pass 
data.  A  summary  of  all  the  parameters  utilized  in  this  section  is  pre¬ 
sented  in  Table  8-3. 

Determination  of  Parameters 

Figure  8-12  is  a  plot  of  ultimate  shear  stress  versus  normal 
stress  obtained  from  the  reverse  pass  of  quasi-static  testing  with  the 
9-degree  surface.  The  square  of  the  "best  fit"  line  through  these 
points  yields  a  value  of  y  for  the  reverse  pass  of  0.45. 

As  previously  discussed,  it  is  assumed  that  the  phase  change 
parameter,  n,  is  related  to  the  characteristics  of  the  microasperity 
pattern  as  opposed  to  the  macroasperities.  Assuming  the  microrsperities 
are  similar  during  travel  up  or  down  the  asperity  sides,  the  same  value 
of  n  -  2.5  is  employed  for  both  the  forward  and  reverse  passes. 

With  the  values  of  y  and  n,  the  growth  function  parameters  may 
be  determined.  Figure  8-13  is  a  plot  of  a  versus  £  from  data  obtained 
during  the  reverse  pass.  Data  from  tests  performed  under  four  levels  of 
normal  stress  5,  10,  20,  and  50  psi  (34.5,  68.9,  137.8,  and  344.5  KPa) 
is  presented.  Using  this  data,  growth  function  constants  a  and  b  are 
0.024  and  -0.149,  respectively.  The  solid  line  in  Fig.  8-13  is  a  plot 
of  the  determined  growth  function. 

Calculations  for  the  nonassociative  parameter  <  during  the  re¬ 
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verse  pass  are  presented  in  Appendix  III. 
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Back  Predictions  of  Quasi-Static  Loading 
With  Unloading  and  Load  Reversal 

Using  the  parameters  discussed  in  the  preceding  section,  back 
predictions  of  both  forward  and  reverse  passes  under  four  levels  of 
normal  stress  on  the  9-degree  surface  are  computed  and  compared  with 
test  data.  Both  associative  and  nonassociative  flow  rules  are  employed. 

Algorithm.  The  algorithm  is  similar  in  concept  to  the  algorithm 
described  in  the  section  on  quasi-static  loading.  The  difference  is 
that  the  back  prediction  continues  after  completion  of  the  forward  pass. 
Upon  completion  of  the  forward  pass,  unloading  takes  place  prior  to  load 
reversal,  as  described  in  Chapter  7.  Unloading  is  considered  to  be 
elastic. 

After  elastic  unloading,  parameters  for  the  reverse  pass  are  in¬ 
put,  and  the  computations  resume  at  Step  2  of  the  algorithm  described  in 
the  section  on  quasi-static  loading.  Using  this  procedure,  parameters 
which  involve  the  history  of  deformation  are  reinitialized  at  the  onset 
of  reverse  loading.  However,  the  history  of  deformation  is  not  erased 
since  parameters  for  the  reverse  pass  were  determined  from  reverse  pass 
test  results  obtained  on  a  surface  that  previously  underwent  deforma¬ 
tions  of  the  forward  pass. 

Results ■  Figures  8-14  to  8-18  present  results  of  back  predic¬ 
tions  of  quasi-static  shear  tests  with  unloading  and  load  reversal.  Back 
predictions  and  test  data  are  presented  considering  shear  with  the 
9-degree  surface  under  different  levels  of  normal  stress  including  5,  10, 
20,  and  50  psi  (34.5,  68.9,  137.8,  and  344.5  KPa) .  Both  associative  and 
nonassociative  flow  rules  are  considered. 
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Figure  8-14.  Test  Data  and  Back  Predictions  of  T  versus  u  ; 
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Figure  8-16.  Test  Data  and  Back  Predictions  of  v  versus  u 
Quasi-Static  Tests  With  Unloading  and  Load  r 
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Quasi-Static  Tests  With  Unloading  and  Load  r 
Reversal;  9°  Surface;  0^  =  20  psi;  Associative 
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Figure  8-18.  Test  Data  and  Back  Predictions  of  v  versus  u 
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Quasi-Static  Tests  With  Unloading  and  Load 
Reversal;  9°  Surface;  a  =50  psi;  Associative 
and  Nonassociative  Cases 


Figure  8-14  is  a  plot  of  shear  stress  versus  ur  with  results 
from  all  levels  of  normal  stress  presented.  The  back  predictions  appear 
satisfactory  with  some  significant  error  evident  at  normal  stress  level 
50  psi  (344.5  KPa) . 

Figures  8-15  to  8-18  are  plots  of  uf  versus  v^  considering  un¬ 
loading  and  load  reversal.  Each  figure  presents  results  from  testing 
under  one  level  of  normal  stress  on  the  9-degree  surface,  including  back 
predictions  utilizing  associative  and  nonassociative  flow  rules  and  ob¬ 
servations.  Back  predictions  of  the  reverse  pass  utilizing  the 
associative  flow  rule  fail  to  indicate  compression  of  the  joint,  and 
predict  dilation,  which  is  not  a  physical  reality.  Improved  results  are 
obtained  with  the  introduction  of  the  nonassociative  flow  rule  although 
significant  errors  still  prevail,  particularly  at  normal  stress  levels 
of  5  and  50  psi  (34.5  and  344.5  KPa).  At  the  5  psi  (34.5  KPa)  normal 
stress  level,  too  much  dilation  is  indicated  at  the  conclusion  of  the 
reverse  pass.  At  the  50  psi  (344.5  KPa)  normal  stress  level,  too  much 
compression  is  indicated.  This  is  similar  to  the  errors  discussed  with 
respect  to  u^  versus  v^  curves,  back  predicted  utilizing  the  nonassocia¬ 
tive  flow  rule,  for  the  case  of  quasi-static  loading.  The  error  at  the 
end  of  the  reverse  pass  is  higher  than  that  at  the  end  of  the  forward 
pass  because  the  errors  of  the  forward  and  reverse  passes  are 
compounded. 


Cyclic  Shear 

Cyclic  shear  is  investigated  with  test  results  from  the  9-degree 


surface  under  normal  stress  levels  of  5,  10  and  20  psi  (34.5,  68.9,  and 
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137.8  KPa)  with  a  maximum  amplitude  of  cyclic  displacement,  uj*,  of 
0.1  in.  (2.54  mm.).  As  described  in  Chapter  6,  this  case  is  modeled  by 
considering  four  distinct  phases  of  displacement  defined  by  direction 
[refer  to  Fig.  6-6 (b)].  By  inspection  of  cyclic  shear  test  results,  it 
appears  that  the  behavior  during  Phases  1  and  3  and  2  and  4  is  similar. 
The  direction  of  travel  is  different,  but  the  tendency  to  compress  or 
dilate  is  identical,  and  the  hardening  behavior  and  mobilized  shear 
stress  at  the  end  of  the  phases  are  very  close.  Since  the  parameters  de¬ 
termined  for  Phases  1  and  3  and  4  and  2  should  be  the  same,  only  two  sets 
of  parameters  need  be  determined  to  describe  the  four  phases  of  displace¬ 
ment  during  a  given  cycle. 

From  the  test  results ,  the  mobilized  shear  stress  at  the  end  of 
Phases  1  and  3  increases  with  number  of  cycles.  This  cyclic  hardening 
shall  be  incorporated  into  the  plasticity  model  by  making  the  ultimate 
parameter,  y,  a  function  of  number  of  cycles.  In  this  manner,  y  includes 
the  effects  of  the  history  of  the  deformation  on  the  behavior  from  cycle 
to  cycle,  while  during  a  given  cycle  the  growth  function,  a,  incorporates 
the  effect  of  history  of  deformation. 

The  manner  in  which  y  is  made  a  function  of  number  of  cycles  as 
well  as  the  determination  of  the  other  parameters  utilized  in  the  model 
is  discussed  here.  An  algorithm  for  predicting  results  of  cyclic  tests 
is  also  presented,  and  a  comparison  of  observed  and  back  predicted  re¬ 


sults  is  made. 
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Determination  of  Parameters 

Within  each  phase  of  a  given  displacement  cycle,  the  same  five 
parameters  discussed  in  the  section  on  quasi-static  loading  are  utilized. 
In  a  sense,  the  isotropic  hardening  model  developed  for  the  case  of 
quasi-static  loading  is  applied  independently  to  each  phase  of  displace¬ 
ment.  In  the  determination  of  parameters,  consideration  is  given  to  the 
history  of  deformation  so  that  the  isotropic  hardening  model  may  incor¬ 
porate  the  effect  of  cyclic  shear.  Table  8-4  presents  a  summary  of  the 
parameters  used  for  the  back  prediction  of  cyclic  shear  tests. 

Ultimate  Parameter  y.  Figure  8-19  presents  plots  of  ultimate 
shear  stress  versus  normal  stress  for  Phases  1  and  3  at  various  displace¬ 
ment  cycles.  The  solid  lines  in  the  figure  are  ultimate  surfaces  fitted 
through  the  data  points  via  a  linear  regression  analysis.  As  before, 
the  square  of  the  slope  of  these  ultimate  surfaces  yields  the  value  of  y. 
The  manner  in  which  y  increases  with  number  of  cycles  is  apparent.  A 
plot  of  y  versus  number  of  cycles,  N,  is  depicted  in  Fig.  8-20.  It 
appears  that  an  expression  of  hyperbolic  form  would  fit  this  data.  An 
expression  for  y  during  Phases  1  and  3  is  determined  as 

y  =  N/(1.67  +  0.7/N)  (8-1) 

The  solid  line  in  Fig.  8-20  is  a  plot  of  this  function  indicating  the 
closeness  of  fit  to  experimental  points. 

The  peak  shear  stress  during  Phases  2  and  4  of  the  cyclic  dis¬ 
placement  is  not  consistently  apparent  from  the  test  results.  It  is 
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Table  8-4. 


Parameters  for  Back  Prediction  of  Cyclic  Test 


9°  Surface 

#  a  =  20  psi 

n 

Urm  -  0.1  in. 


Phase  1 

Cycle  1 

Phase  1  (Njfcl) 
and  3 

Phase 

2  +  4 

Y 

N/  (1 . 67  +  0.71  N) 

N/(1.67  +  .71  N) 

0.228 

n 

2.5 

2.5 

2.5 

a 

0.00455 

N/(200  +  50  N) 

0.0096 

b 

-0.364 

-0.90 

-0.206 

< 

0.60 

0.60 

1.20 

N  is  number  of  cycles. 


I 


I 
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Figure  8-19.  Ultimate  Envelopes  from  Cyclic  Test  Results  at 
Different  Cycles  of  Displacement  (N)  for 
Phases  1  and  3;  9°  Surface 


assumed  that  negligible  surface  wear  occurs  during  "down  asperity  level" 
and  that  the  maximum  mobilized  shear  stress  does  not  vary  with  number  of 
cycles.  This  assumption  is  also  made  by  Plesha  (1986).  Since  the  ul- 

i 

timate  shear  stress  is  computed  directly  from  the  peak  shear  stress,  it 
also  does  not  vary  with  number  of  cycles.  The  value  of  y  is  determined 
from  Fig.  8-21  which  is  a  plot  of  ultimate  shear  stress  versus  a  from 
data  obtained  during  Phases  2  and  4  of  Cycle  1.  The  ultimate  surface 
passed  through  these  data  points  yields  a  value  for  y  of  0.228. 

Phase  Change  Parameter  n.  The  same  value  of  n  =  2.5,  which  was 
employed  for  the  two  previously  discussed  cases,  will  be  used  for  the 
case  of  cyclic  shear. 

Growth  Function  Constants.  Some  special  considerations  in  de¬ 
termining  the  proper  growth  function  for  different  phases  during  a 
particular  cycle  of  displacement  need  to  be  discussed.  During  deforma¬ 
tion  of  a  particular  phase,  hardening  is  described  through  the  use  of 
the  growth  function.  The  question  arises,  "What  is  the  value 
of  the  growth  function  at  the  beginning  of  a  particular 
phase?" 

For  the  first  phase  of  Cycle  1,  the  answer  is  more  obvious  than 
at  other  points  in  the  deformation  history.  Since  no  previous  shear 
displacement  has  occurred,  the  initial  yield  surface  should  be  that 
corresponding  to  the  state  of  stress  after  application  of  the  normal 
stress  but  prior  to  the  application  of  shear  to  the  joint  plane. 

At  the  beginning  of  Phase  2,  previous  deformation  has  occurred 


during  Phase  1  which  must  somehow  be  accounted  for  in  the  model.  During 
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Phase  2,  unloading  is  assumed  to  be  elastic.  Therefore,  the  state  of 
stress  at  the  onset  of  plastic  deformations  is  once  again  normal  stress 
in  the  absence  of  shear  along  the  joint  plane.  Thus  the  initial  yield 
surface  for  Phase  2  should  reflect  this  state  of  stress.  Previous  de¬ 
formation  of  the  joint  surface,  which  occurred  during  Phase  1,  is 
considered  since  determination  of  parameters  for  Phase  2  is  based  on  de¬ 
formations  which  occur  to  the  surface  after  being  disturbed  during 
Phase  1.  The  same  arguments  presented  in  this  paragraph  hold  for 
Phase  4  but  with  respect  to  Phase  3. 

During  Phase  1  (after  Cycle  1)  and  Phase  3,  the  situation  is  a 
bit  different.  No  :r.  oading  in  shear  occurs  at  the  beginning  of  these 
phases,  and  in-t'^ily  the  state  of  stress  is  that  which  exists  at  the 
end  of  Phases  2  and  4,  respectively.  Therefore,  the  value  of  the  growth 
function  at  the  beginning  of  Phase  1  (after  Cycle  1)  and  Phase  3  must  be 
such  that  the  initial  yield  surfaces  reflect  this  state  of  stress. 

A  program  was  developed  which  computes  the  growth  function 
parameters  for  each  phase  of  displacement  at  various  cycles.  Required 
input  for  each  cycle  of  interest  consists  of  the  values  of  y,  n,  the 
elastic  shear  stiffness,  and  data  points  consisting  of  shear  stress, 
normal  stress,  and  relative  vertical  and  tangential  displacements  from 
each  phase  of  the  displacement  cycle.  The  program  computes  components 
of  plastic  deformation,  the  value  of  the  growth  function,  a,  and  the 
trajectory  of  plastic  strain  for  each  data  point  of  a  phase.  The  program 
then  creates  a  data  file  for  plotting  a  versus  E,  and  performs  a  linear 
regression  analysis  of  logarithmic  transformed  data  points  to  determine 
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growth  function  constants  for  each  phase  of  displacement.  The  program 
is  entitled  "GROCYC. FOR. "  Table  8-5  shows  growth  function  constants  a 
and  b  determined  by  the  program  from  cyclic  shear  test  data  for  the 
9-degree  surface  with  a  normal  stress  of  20  psi  (137.8  KPa) . 

Growth  function  parameters  are  determined  based  on  the  values  of 
the  ultimate  parameters  used  for  a  given  cycle  of  displacement  and  may 
be  functions  of  number  of  cycles  themselves.  Growth  function  parameters 
from  Phases  1  and  3  are  a  function  of  number  of  cycles  while  those  from 
Phases  2  and  4  are  considered  constant.  For  Phases  2  and  4,  average 
values  of  the  growth  function  parameters  determined  from  different  cycles 
are  employed.  Expressions  and  values  for  growth  function  parameters  em¬ 
ployed  in  the  back  predictions  may  be  seen  in  Table  8-4.  Only  data  from 
the  test  with  a  normal  stress  of  20  psi  was  used  to  determine  growth 
function  constants,  and  these  were  used  for  back  prediction  of  all  tests 
considered  herein. 

Figures  8-22  to  8-28  are  plots  of  a  versus  E,  for  various  phases 
and  cycles  of  displacement.  Experimental  data  points  and  plots  of  the 
growth  function  using  both  constants  determined  only  from  data  pertaining 
to  the  cycle  of  interest,  and  those  determined  with  expressions  given  in 
Table  8-4  are  presented.  From  this  information,  the  proximity  of  the 
employed  growth  function  to  observation  may  be  ascertained. 

Nonassociative  Parameter  k.  Computations  for  the  values  of  the 
nonassociative  parameter  <  for  different  phases  and  cycles  of  displace¬ 
ment  are  presented  in  Appendix  III.  Since  the  plots  of  ur  versus  vr  do 
not  vary  from  cycle  to  cycle,  the  computations  do  not  reflect  significant 
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Table  8-5.  Growth  Function  Constants  Determined  for 
Cycles  of  Displacement;  9°  Surface; 
a  =  20  psi 


Cycle 

Phase  1  & 

a 

3 

b 

Phase  2  &  4 

a  b 

1 

0.0000 

-4.22 

0.0096  -0.206 

2 

0.0096 

-0.851 

5 

0.0116 

-1.011 

10 

0.0174 

-0.914 

20 

0.0157 

-1.028 

100 


0.0211 


-0.892 


Figure  8-23.  a  versus  £;  Phase  3;  Cycle  1 
Cyclic  Test:  9°  Surface; 
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variation  of  <  with  N.  Some  variation  may  be  predicted  between  cycles 
for  Phases  1  and  3  due  to  variation  in  Y,  but  this  does  not  appear  to  be 
significant.  There  is,  however,  significant  difference  in  the  values  of 
K  between  Phases  1  and  3  and  Phases  2  and  4.  Therefore,  only  two  values 
of  <  are  required  for  all  cycles  of  displacement,  one  for  Phases  1  and  3 
and  one  for  Phases  2  and  4. 

Back  Predictions  of  Cyclic  Shear  Tests 

Using  the  parameters  discussed  in  the  preceding  section,  back 
predictions  of  cyclic  shear  tests  with  the  9-degree  surface  are  per¬ 
formed.  The  computations  are  carried  out  for  100  cycles  of  displacement, 
and  compared  with  test  data.  Both  associative  and  nonassociative  flow 
rules  are  employed. 

Algorithm.  This  section  describes  the  algorithm  used  to  deter¬ 
mine  increments  of  shear  stress  and  relative  vertical  displacement  based 
on  applied  increments  of  relative  tangential  displacement  under  the  con¬ 
dition  of  constant  normal  stress.  The  algorithm  describes  the 
determination  of  the  initial  yield  surface  for  a  given  phase;  i.e.,  the 
initial  value  for  the  growth  function,  the  formulation  of  the  elastic- 
plastic  constitutive  matrix,  and  determination  of  the  unknown  quantities. 

At  the  beginning  of  each  phase  of  a  cycle,  parameters  are  deter¬ 
mined  from  appropriate  expressions.  Using  these  parameters,  the 
following  algorithm  is  employed  during  each  phase  of  the  deformation 
history: 

1.  Calculate  the  value  of  the  growth  function  at  the 


beginning  of  the  phase. 
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a.  Initially  for  Phase  1  of  Cycle  1,  the 
growth  function  is  determined  so  that 
the  initial  yield  function  is  equal  to 
zero.  The  shear  stress  is  zero  and  the 
normal  stress  is  the  applied  normal 
stress . 

b.  For  Phases  2  and  4,  after  elastic  unloading 
and  just  prior  to  load  reversal,  the  initial 
value  for  the  growth  function  is  computed  as 
in  a. 

c.  At  the  beginning  of  Phase  1  after  Cycle  1, 
and  Phase  3,  the  value  of  the  growth  function 
is  determined  based  on  the  history  of  de¬ 
formation  at  the  end  of  Phases  4  and  2, 
respectively.  However,  if  this  value  is 
greater  than  that  required  to  satisfy  the 
condition  that  the  current  yield  surface 
equals  zero,  the  growth  function  is  modified 
so  that  this  condition  is  satisfied  given  the 
current  state  of  stress. 

d.  The  growth  function  at  any  point  during  the 
deformation  of  any  phase  of  any  cycle  other 
than  the  initial  condition  is  computed  based 
on  the  history  of  deformation  at  the  end  of 


the  previous  Increment. 
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2.  Set  the  value  of  a  equal  to  the  value  of  a  just 
after  the  application  of  normal  stress,  but  prior 
to  the  initiation  of  shear.  At  this  point,  £  is 
equal  to  £v  and  r^  is  unity. 

3.  Using  these  values  of  a,  a  ,  r^  and  the  parameters 
for  the  particular  phase  and  cycle,  compute  a^,  the 
derivatives  of  the  yield  and  potential  functions, 
and  form  the  K6*5  matrix. 

4.  Apply  an  increment  of  relative  tangential  displace¬ 
ment,  Du  *  0.0001  in. 

5.  Solve  for  increment  of  shear  stress,  dt,  and  rela¬ 
tive  vertical  displacement,  dv. 

6.  Compute  elastic  and  plastic  components  of  displacement 
using  the  associative  or  nonassociative  flow  rule. 

7.  Obtain  updated  values  for  £,  £v»  r^  and  a.  If  this  is 
the  first  increment  of  shear,  check  to  see  that  F  *  0, 
and  if  F  is  negative,  adjust  the  elastic  shear  stiff¬ 
ness,  as  described  in  the  algorithm  for  quasi-static 
loading  such  that  F  >  0. 

8.  If  the  phase  is  not  over,  go  to  (3).  If  the  phase 
is  over,  go  to  (1).  If  this  is  the  end  of  Phase  4, 

Cycle  100,  stop. 

This  algorithm  was  coded  into  a  program  entitled  CYCDIS,  which 


was  then  used  for  back  predicting  results  from  cyclic  shear  tests. 
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Results .  Figures  8-29  to  8-33  present  results  of  back  predic¬ 
tions  of  a  cyclic  shear  test  performed  on  the  9-degree  surface  with 
an  *  20  psi  (138  KPa) .  Back  predictions  are  for  the  case  of  a  nonasso- 
ciative  flow  rule  only.  Results  with  a  ■  5  and  10  psi  (34.5  and 
69  kPa)  may  be  found  in  Appendix  II. 

Four  figures  of  shear  stress  versus  ur  are  presented  (Figures 
8-29  to  8-32) .  The  first  of  these  depicts  back  predictions  of  cycles 
ranging  from  1  to  100,  and  the  remaining  compare  back  predictions  to 
test  data  for  Cycles  1,  10  and  100.  The  transition  between  phases  is 
readily  apparent  in  these  back  predictions.  Comparison  of  back  predic¬ 
tion  to  observation  is  considered  satisfactory  although  there  is  some 
error  present  for  Cycle  1.  This  is  due  to  the  function  which  describes 
the  variation  of  y  with  number  of  cycles  [Eq.  (8-1)].  Inspection  of 
Fig.  8-20,  which  compares  this  function  to  values  of  y  found  for  each 
cycle,  indicates  that  the  error  between  the  data  point  and  the  "fitted" 
curve  is  greater  for  Cycle  1  than  for  Cycle  100.  It  is  conservative  to 
underpredict  the  shear  strength  of  the  first  cycle;  therefore,  the  rela¬ 
tionship  which  describes  the  increase  in  y  with  number  of  cycles  is 
considered  satisfactory. 

Results  of  vf  versus  u^,  Fig.  8-33,  present  back  predictions 
for  cycles  ranging  from  1  to  100,  and  compare  these  to  test  data.  Since 
observations  do  not  vary  significantly  from  cycle  to  cycle,  only  observa¬ 
tions  from  Cycle  1  are  presented.  Although  back  predictions  of  Cycle  1 
are  acceptable,  considerable  errors  begin  to  predominate  beyond  the 
tenth  cycle  of  displacement.  It  seems  that  at  higher  cycles  the  vertical 


—0.23  -0.20  -0.13  -0.10  -0.03  0.00  0.03  0.10  0.15  0.20  C 

u  /*  ( IN) 


Figure  8-30.  Test  Data  and  Back  Predictions;  T  versus  u 
Cyclic  Test;  9°  Surface;  0  =20  psi; 
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displacements  during  Phases  1  and  3  are  not  predicted  well  so  that  at 
the  conclusion  of  the  cycle  a  significant  net  error  persists.  Since  ob¬ 
servations  indicate  no  perceptible  increase  or  decrease  in  joint 
thickness  at  the  end  of  a  cycle,  these  errors  were  corrected  at  the  end 
of  a  cycle  such  that  errors  would  not  accumulate  during  back  prediction. 


CHAPTER  9 


SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

Summary 

A  constitutive  model  to  describe  the  behavior  of  individual  rock 
joints  during  shear  has  been  proposed.  The  model  is  based  on  the  theory 
of  incremental  plasticity  and  is  a  specialization  of  the  Generalized 
Plasticity  Model  developed  for  geologic  solid  materials.  In  its  present 
form,  the  joint  model  can  consider  initial  state  of  stress,  states  of 
shear  and  normal  stress,  plastic  hardening,  volume  changes  at  joints, 
monotonic  shear  loading,  unloading  and  load  reversal,  and  cyclic  shear 
loading. 

A  comprehensive  program  of  shear  testing  was  performed  on  ideal¬ 
ized  joint  specimens.  In  order  to  perform  these  tests,  certain 
modifications  to  the  CYMDOF  device  were  necessary.  These  modifications 
allowed  the  rock  joint  to  compress,  dilate,  or  rotate  in  the  plane  of 
shear.  Quasi-static  and  cyclic  shear  tests  were  performed  under  variable 
conditions  of  joint  geometry,  initial  state  of  stress,  and  cyclic 
displacement. 

Inductance  coils  were  employed  for  the  measurement  of  vertical 
displacement  at  the  center  of  the  joint  plane.  To  the  author's  knowl¬ 
edge,  this  is  the  first  time  these  devices  have  been  employed  for  this 
purpose.  Special  calibration  and  verification  of  the  performance  of 
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these  instruments  was  necessary.  The  effects  of  coil  misalignment 
during  shear  and  the  presence  of  the  metal  sample  holders  were  con¬ 
sidered  in  calibrating  the  instruments  for  this  application. 

Parameters  for  the  plasticity  model  were  determined  from  the 
test  data  and  back  predictions  of  the  test  data  were  performed.  In 
order  to  determine  parameters  for  cyclic  tests,  the  displacement  cycle 
was  considered  to  consist  of  different  phases  of  displacement.  It  was 
demonstrated  that  during  back  predictions  this  amounts  to  collapse  or 
expansion  and  rotation  of  the  yield  surface  at  appropriate  instants  of 
time  during  the  deformation  process.  Cyclic  shear  test  results  show  an 
increase  in  mobilized  peak  shear  stress  with  number  of  cycles.  The  in¬ 
crease  in  mobilized  shear  strength  with  number  of  cycles  was  described 
by  making  the  ultimate  parameter,  Y,  a  function  of  number  of  cycles. 
Programs  were  written  for  determination  of  parameters  and  back 
predictions. 

Conclusions 

In  conclusion,  the  plasticity  model  described  in  Chapter  7  is 
capable  of  predicting  the  behavior  of  rock  joints  during  monotonic  or 
cyclic  shear  loading.  Displacement  induced  anisotropy  and  cyclic 
hardening  behavior  are  incorporated  in  the  model.  Back  predictions  pre¬ 
sented  in  Chapter  8  indicate  that  shear  response  of  the  joint  is 
predicted  well  for  both  monotonic  and  cyclic  shear.  To  reasonably  pre¬ 
dict  the  vertical  behavior,  a  nonassociative  flow  rule  must  be  employed 
or  dilation  will  be  overpredicted.  Predictions  of  vertical  response 
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indicate  that  the  model  is  much  more  sensitive  to  changes  in  normal 
stress  than  test  results  indicate.  Modifications  are  necessary  to  im¬ 
prove  these  predictions. 

The  inductance  coils  utilized  for  measurement  of  dilation/ 
compression  provide  a  very  useful  technique  for  taking  measurements 
within  the  interface.  The  instruments  must  be  calibrated  for  misalign¬ 
ment  and  the  presence  of  metal  sample  holders,  as  described  in 
Chapter  4. 


« 


Recommendations 

The  current  research  has  considered  the  case  of  cyclic  shear  on 
dry,  unfilled,  structured  joint  surfaces.  Although  this  is  a  signifi¬ 
cant  contribution,  many  improvements  and  modifications  in  the  model  are 
required  before  the  model  can  be  applied  to  a  wide  range  of  problems. 

It  is  anticipated  that  these  improvements  and  modifications  will  build 
upon  the  current  model,  utilizing  the  hierarchical  approach,  as  outlined 
by  Desai  et  al.  (1986).  To  accomplish  this,  the  following  improvements 
are  recommended: 

1.  The  model  needs  to  be  modified  to  describe  the  effect  of 
normal  stress  better.  In  the  present  model,  the  growth  function  was  de¬ 
termined  mainly  on  the  basis  of  deformation  which  occurs  during  shear. 

A  split  form  of  the  growth  function,  similar  to  that  proposed  by  Faruque 
(1983) ,  in  which  a  portion  is  dependent  on  normal  stress  may  predict 
the  response  more  accurately  at  all  levels  of  normal  stress. 
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2.  The  model  should  be  extended  to  the  case  of  strain  softening. 
It  may  be  possible  to  describe  strain  softening  by  using  the  model  pro¬ 
posed  by  Frantziskonis  and  Desai  (1986)  in  which  softening  is  described 
by  considering  some  part  of  the  material  (joint)  to  be  damaged.  The  be¬ 
havior  of  the  undamaged  zone  would  still  be  described  with  the  isotropic 
hardening  model  described  in  this  report. 

3.  Tests  should  be  performed  on  natural  rock  joints  using  the 
CYMDOF  device  to  determine  parameters  for  the  plasticity  model.  The 
parameters  should  be  compared  to  those  obtained  from  testing  of  idealized 
joints.  Back  predictions  should  be  performed  to  show  that  the  model  is 
still  easily  applied  to  the  case  of  naturally  occurring  rock  joints. 

4.  The  effects  of  filler  materials  in  rock  joints  should  be  con¬ 
sidered.  Naturally  occurring  rock  joints  may  be  filled  with  sand  and 
clay  which  leach  into  the  joint  space. 

5.  The  effect  of  water  in  the  joint  plane  should  be  considered. 
Natural  joints  may  be  saturated  with  water.  If  free  drainage  is  not 
allowed,  pressure  variations  can  occur  during  shear  due  to  the  tendency 
to  compress  or  dilate.  Changes  in  water  pressure  then  contribute  to 
changes  in  effective  stress  at  the  joint  plane. 

6.  The  model  should  be  implemented  in  a  finite  element  program. 

7.  Boundary  value  problems  in  which  the  presence  of  joints  is  a 
factor  should  be  studied.  Measurements  could  be  made  in  the  field  or 
laboratory.  Results  from  analysis  using  methods  which  incorporate  the 
proposed  joint  model  could  be  compared  with  measured  results  for  verifi¬ 
cation  purposes. 
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APPENDIX  I 


DERIVATION  OF  K®P  MATRIX 


Define  the  following: 


{a}  -  {°n} 


{e}  -  elastic  strain 


{e}  =*  {  r}  *  -r  {e}p  -  plastic  strain 

u_  t 


K  0 
n 


[K1  -  [„  R  ] 


t  -  joint  thickness  taken  as  unity 


The  stress-elastic  strain  relationship  may  be  written  as 


{da}  =  [K]  [def 


Decompose  the  total  strain  into  elastic  and  plastic  parts 


(I-U 


{de}  =  {de}e  +  {de}p 


d-2) 


Rewrite  Eq.  (1-1)  as 


{da}  =  [K]  ({de}  -  {de}p) 


(1-3) 


Associative  Case 


Assuming  associated  plasticity,  the  flow  rule  may  be  written  as 


-  *  t£> 


(1-4) 
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Write  the  consistency  condition,  assuming  a  is  a  function  of  £, 


dF  -  {|£}  {da}  +  ||  d£  =0 


(1-5) 


T  1/2 

dC  *  ({de}p  (de}p) 


and  Eq.  (1-4) 


dC  -  Ay. 


(1-6) 


where 


3F  T  3F 

<!§» 


Utilizing  Eqs.  (1-3),  (1-4),  (1-5)  and  (1-6),  obtain  an  expression  for 
A  as  follows: 


{da}  -  [K]  ({de}  -  A  {||}) 

dF  -  tf}1  [K]  ((de)  -  x|I)  +  ||x  YF  -  o 


<£>  <■*> 
(1I}T  [K1  fl£l  -  II  Y 

^3a^  [KJ  13crf  3C  YF 


(1-7) 


By  substituting  Eqs.  (1-4)  and  (1-7)  into  Eq.  (1-3),  the  following  ex¬ 
pression  is  obtained: 


{1I}T  rK1  {11}  _  H  Y 

l3a'  LKJ  Ha'  3£  YF 


-)  {de} 


Equation  (1-8)  may  be  rewritten  as 


where 


{da}  =  [K]ep  {de} 


M  <|f>  t||>  IK) 

<H>T  «>  <i>  - 1  yf 


(1-10 


Nonassociative  Case 

Assuming  nonassociated  plasticity,  the  flow  rule  may  be  written 


{de}p  =  A  {|S} 


(1-11 


and  Eq.  (1-5)  becomes 


dF  =  {-^}  [K]  ({de} 


(1-12 


Defining 


Yr  -  <>  <§» 


and  utilizing  Eqs.  (1-3),  (I— 11),  (1-12)  and  (1-6),  an  expression  for  A 


is  obtained  as 


APPENDIX  II 


BACK  PREDICTION  FIGURES 

Appendix  II  contains  figures  for  determination  of  parameters 
and  results  including  test  data  and  back  predictions  which  were  not  dis¬ 
cussed  in  the  main  body  of  the  text.  These  include  results  from 
quasi-static  testing  with  flat,  5-degree  and  7-degree  surfaces,  as  well 
as  cyclic  testing  with  the  9-degree  surface  at  normal  stresses  of  5  and 
10  psi  (34.5  and  69  KPa) . 
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Figure  II- 1.  Ultimate  Envelope;  Quasi-Static  Test  Series, 
Forward  Pass;  Flat  Surface 


Figure  I I— 1 1 .  Test  Data  and  Back  Predictions  of  T  versus  u 
Quasi-Static  Loading;  7°  Surface  1 
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Back  Predictions 


0.1590 


u  (in.) 
r 

Figure  11-15.  Test  Data  and  Back  Predictions  of  v  versus  u  ; 

Quasi-Static  Loading;  a  =10  psi;  r  r 

5°  Surface;  Associativenand  Nonassociative  Case 
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0.1990 


Test  Data  S 
Back  Predictions 


Figure  11-16.  Test  Data  and  Back  Predictions  of  vr  versus  ur; 
Quasi-Static  Loading;  0  =20  psi; 

5°  Surface;  Associative  and  Nonassociative  Case 
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Figure  11-20.  Test  Data  and  Back  Predictions  of  vr  versus  u^.; 
Quasi-Static  Loading;  0  =20  psi; 

7°  Surface:  Associativeand  Nonassociative  Case 
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Back  Predictions 


Figure  11-25.  Test  Data  and  Back  Predictions  of  v  versus  u 
Quasi-Static  Loading;  Average  Parameters; 

Flat  Surface;  Nonassociative  Case 


(in. 


Figure  11-27.  Test  Data  and  Back  Predictions  of  v  versus  u 
Quasi-Static  Loading;  Average  Parameters; 

5°  Surface;  a  =50  psi;  Nonassociative  Case 


Figure  11-28.  Test  Data  and  Back  Predictions  of  v  versus  u 
Quasi-Static  Loading;  Average  Parameters; 

7°  Surface;  a  =5  psi;  Nonassociative  Case 


0.1940 


Figure  11-29.  Test  Data  and  Back  Predictions  of  v  versus  u  ; 

Quasi-Static  Loading;  Average  Parameters; 

7°  Surface;  a  =50  psi;  Nonassociative  Case 
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Figure  11-35.  Test  Data  and  Back  Predictions;  T  versus  u  ; 

Cyclic  Test;  9°  Surface;  0  =10  psi;  Cycl 
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Test  Data  and  Back  Predictions;  v  versus 
Cyclic  Test;  9°  Surface;  cr  =5  psi 
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APPENDIX  III 


Derivation  of  these  expressions  can  be  found  in  Chapter  7  and  all  parame¬ 
ters  have  been  previously  defined. 

Tables  III-l  and  III-3  show  the  determination  of  <  from  test  data 
and  the  parameters  used  in  the  calculations.  The  parameter  <  is  deter¬ 
mined  using  data  from  the  following: 

1.  The  forward  pass  of  the  quasi-static  test  series 
(Table  III-l). 

2.  The  reverse  pass  of  quasi-static  testing  with  the 
9-degree  surface  (Table  III-2) . 

3.  Cyclic  testing  at  different  cycles  and  phases  of 
displacement  (Table  III-3) . 
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Table  III-l .  Determination  of  <;  Quasi-Static  Test 
Series,  Forward  Pass 


Surface 

Y 

a 

(psx) 

T 

(psi) 

(V 

du 

u 

aQ 

ai 

r 

V 

< 

Flat 

0.364 

8.10 

4.88 

0.0 

0.1023 

0. 1279 

0.0000 

0.7999 

5° 

0.42 

5.0 

3.24 

-0.077 

0.1324 

0.1887 

0.0593 

0.746 

0.42 

10.0 

6.48 

-0.065 

0.0956 

0.1334 

0.0515 

0.755 

0.42 

20.0 

12.96 

-0.066 

0.0674 

0.0943 

0.0502 

0.753 

0.42 

50.0 

34.40 

-0.076 

0.0419 

0.0551 

0.0551 

0.744 

7° 

0.783 

5.0 

4.42 

-0.080 

0.2547 

0.3502 

0.0872 

0.797 

0.783 

10.0 

8.85 

-0.064 

0.1839 

0.2476 

0.0860 

0.812 

0.783 

20.0 

17.69 

-0.057 

0.1310 

0.1751 

0.1013 

0.832 

0.783 

50.0 

44.24 

-0.077 

0.0809 

0.1107 

0.0946 

0.806 

9° 

0.81 

5.0 

4.50 

-0.110 

0.2544 

0.3622 

0. 1022 

0.782 

0.81 

10.0 

9.00 

-0.0995 

0.1823 

0.2561 

0.1050 

0.795 

0.81 

20.0 

18.00 

-0.090 

0.1303 

0.1811 

0.1208 

0.818 

0.81 

50.0 

45.00 

-0.086 

0.0839 

0.11455 

0.1182 

0.820 

n  -  2.5  for  all  tests. 
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Table  III-2.  Determination  of  K ;  Reverse  Pass; 

9°  Surface 

I 


Surface 

Y 

a 

n 

psi 

T 

u 

psi 

dv.  P 

u 

aQ 

ai 

r 

V 

< 

9° 

0.45 

5.0 

-3.35 

-0.0979 

0.1634 

0.2012 

0.1100 

0.91 

0.45 

50.0 

-33.54 

-0.1052 

0.0589 

0.0636 

0.1200 

1.05 

I 


< 


< 


I 


« 


Table  III-3.  Determination  of  c;  Cyclic  Test; 

9°  Surface;  a  =20  psi 
n 


Phase 

Y 

a 

n. 

psi 

T 

u . 
psi 

(— )P 
w 

u 

“Q 

ai 

r 

V 

< 

Cycle 

_1 

1 

0.42 

24.0 

15.60 

-0.221 

0.0450 

0.0857 

0.1126 

0.59 

2  and  4 

0.228 

20.0 

±9.54 

-0.0507 

0.0450 

0.0510 

0.0904 

0.970 

3 

0.42 

24.0 

-15.60 

0.202 

0.047 

0.0857 

0.1259 

0.625 

Cycle 

100 

1  and  3 

1.376 

20.75 

24.34 

-0.309 

0.1780 

0.3021 

0.1453 

0.689 

n  =  2.5  for  all  cycles  and  phases  of  displacement 
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